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SECTION I 
INTRODUCTION 


One principal task of assessing aircraft vulnerability 
to blast loadings is to identify appropriate damage modes and 
to use this information for extension to predicting levels of 
failure for given loadings. For example, an aircraft ina 
tactical air or ground situation can be modelled for first 
analysis as a free free beam or beam on multiple elastic sup- 
ports when considering large scale structural damage. Con- 


sideration of a further localized assessment of damage to a 


part of the fuselage or control surface can require addi- 


tional spatial coordinates to define the structural member. 


Thus a complete identification of the potential damage toler- 
ance levels in aircraft requires appropriate modelling at 
selected structural component levels. A general representation 
of this problem is shown in Figure 1 which classifies loading 
types in a broad sense and identifies the area of primary con- 
cern in this report. 

Very specifically it is the purpose of the following re- 
port to assess those theories most applicable to predicting 
both elastic and plastic dynamic response of parked aircraft, 
modelled as stiffened metal cylindrical shells, to lateral 


blast loads. The structural response so generated can then be 
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used to more extensively study the local effects which ulti- 
mately govern failure of a structural configuration, such as, 
end fixation, rivet connections and overall panel separation 
and perforation. 

In reviewing the literature for analytical predictors for 
such loadings an extensive search utilizing DOD, DDC, NASA, and 
GRA filings was surveyed. As a result of this literature re- 
view, the following predictive type models of Table I emerged 


and were considered for modification. 


TABLE I. LIST OF CYLINDRICAL SHELL ANALYSES 
WITH BLAST LOADINGS 


Classification of 
Investigator Model References 


Semi-Analytical 


Analytical Force 
Equilibrium 


Analytical Modal 
Analysis 


Schuman (1-5)* has tested over six hundred monocoque shell 
configurations subjected to lateral blast pressures at varying 


*Numbers in parentheses () are references and numbers in brac- 
kets [] represent equation numbers. 
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stand off distances and quantified his work in tabular form. 


These tests represent the most exhaustive experimental studies 
found by the current investigators. The data obtained by 
Schuman has been used by Greenspon to check his theory which 

in this report has been labelled as semi-analytical. The 

above classification was made to emphasize that the analytical 
development requires some information based upon experimental 
data to quantify failure modes and/or loadings. The analytical 
code DEPROP developed by Mente and Lee (18) represents an ambi- 
tious analytical representation directed towards predicting the 
time dependent response of the structure and thus obtain a 
history of the response leading up to ultimate structural 
failure. 

Each of the above approaches is useful in predicting cer- 
tain features of experimentally observed shell response which 
is shown in the accompanying Figures 2-4. The particular type 
of response predicted is inherent in the basic assumptions with- 
in the scope of each of the theories presented and these are 
summarized in Table II. 

The major objective of this study was to investigate 
effects of axial stiffening of cylindrical shell subjected 
to transverse blast loads. It was decided the most direct 
approach would be one of modifying an existing model(s) to 


include axial stiffening. In that the DEPROP code of Mente 
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and Lee (17,18) was on file it was selected to be modified 

and serve as the detailed method for calculating localized 
stresses, strains, and deflections. The method of Greenspon 
(6-15) was chosen as a 'first-cut' design type predictor. 

Both methods are based on an energy approach and were used with 
some reasonable accuracy in predicting damage in unstiffened 
cylindrical shells. The selection of only these two models to 
be modified in no way reflects upon the usefulness of the 
other models, but was simply considered the most expedient 
engineering method. 

The analytical development for the response of the cylin- 
drical shell with axial stiffening for the two separate methods 
will be given in the following Sections II and III. In that 
the modification must parallel the original basic method only 


the stiffened or modified version will be developed. The 


basic unstiffened case will be included in the stiffened case 
and a solution for the unstiffened case may be obtained by 


assuming zero number of stiffeners. 
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SECTION II 
EXTENSION OF GREENSPON SHELL THEORY TO INCLUDE STIFFENERS 


2.1 Introduction 

In this section, the semi-analytical model of Greenspon 
(15), previously described in Section I, has been extended to 
include the use of stiffeners in blast loaded shell analysis. 
The selection of this analytical approach for further con- 
sideration and development is based upon (a) its adaptability 
to facilitate an engineering approach for inclusion of stiff- 
eners, (b) the models realistic characterization of the princi- 
pal modes of deformation occurring (collapse and/or circunm- 
ferential buckling) and (c) the relative ease of characterizing 
failure from design curves. 

It should be remarked that the inclusion of stiffeners in 
any existing analytical development should necessitate the 
introduction of anisotropic constitutive equations into the 
formal shell governing equations which in turn would lead to 
complicated coupled equations involving the shells displacement 
coordinates for solution. The approach adopted here, as men- 
tioned, has been based upon an engineering approach as intro- 
duced in reference (15). This technique incorporates the 
influence of stiffness anisotropy in an uncoupled manner into 


the energy equations which is readily adaptable to Greenspon's 
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analytical model. This development, which incorporates details 
of the unstiffened shell analysis is given in the following 


paragraphs. 


2.2 Analytical Development 

This development is based on the work of Greenspon (6-15) 
for the cylindrical shell as shown in Figure 5. The develop- 
ment parallels the case of unstiffened cylindrical shells and 
the modifications to include stiffening effects are shown with 


broken underlines in the respective terms. It is noted that 


full credit for the basic unstiffened method is given and is 

documented in References 6 to 15. The choice of words in the 

following text may in many cases closely parallel those as 
given in the references and when given are selected for clarity 


and for lack of more descriptive method of presentation. 


2.2.1 Potential Energy 
The work of deformation, V per unit volume of an elastic- 


plastic body can be written (13,14), 


V - [ ode, + Ke, (1) 


where K, is the bulk modulus and 


‘ ; 

Ci = Ge -GiGz + Gy" + 3%yF = | 

ot ee 2 t 2 : | 
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Direction of Loading 


Stiffener 


Figure 5. Coordinate System for Greenspon Analysis. 


Direction of Loading 


Stiffener 


Figure 6. Coordinate System for Mente Analysis. 
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@ = ey + Ey + & , [3] 
* g, “3k, ) Ey re, “ZR, ) Yas, = x - 23 Zz) C4a,b,c] 
- ou Oury* ee ; 

E,= yan +($3) , &= + s -2+ £(2% 2 [Sa,b] 
= oe tou O Ww 

u se taby + Oe 2, [Sc] 


: C6a,b,c] 


For very large deformations under intense lateral load- 
ing the midsurface strain involving w should probably be greater 
than the linear terms involving u and v. Assuming that u and 


v and their derivatives are much smaller than w and its deriva- 


tives, 
- L dur) _ 4 dtu 
x= 2 (an) ~4 due e | 
=a Pe 2 | 
&y= 2 Fev a ~2. pag a | 
+ dee gar. O*ur 
ony = Oe ade Si ow C7] 


Further, the material response is restricted to one that 
obeys an elastic-linear hardening law as shown in Figure 7. 


For an incompressible material, that is, @ = 0, the stresses 
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ee 


e, é; 


Figure 7. Stress-Strain Curve. 


can be written in terms of strains, following Reference (15), 


for a Poisson ratio of 0.5 for both the elastic and plastic 


case. 
Ox ss: 52, = (Ey + hie [8a] 
Oy = = = (6, + S#), Ceb] 
Cy = sa oo Cec] 
in which the stress-strain law can be written as, 
Se = E[!-o@)J, C9] 
where 


He: ) = A(i- 3), 


C10] 
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é 
For the elastic case WG.) = Ad-=-\ 20, 6:26, ,6:<@, 
and 
é 
for the plastic case 22, ) = AGze) » Ese 4 Gre. 
Substituting Equation [10] into Equation [1] for the 
case of ¢ = 0, the work done by the internal forces of the 


stiffened cylindrical shell (pg. 5, Ref. 13) is 


V JSS : Ee, [i-AC -2)\de, adadnd3 
DSSS Be-t- -r(i- 2) ]de: dn. aa 


C11] 


= Var + Met 
where Aj indicates integration over the cross sectional area 
of the ith stiffener. The first term on the left side of 
Equation [11] is associated with the work of the shell Vgp 
while the second term represent the contribution of the stiff- 
ener. The dotted underlined terms refer to the stiffener 
contributions. Considering Figure 1 and Equations (3] and [10] 
Vsph may be obtained in the form: 


2 sy ee (i-A) + EA é,6, | adodn Az 


[12] 


~EA&,TAaLt. 
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Rewriting e; of Equation [26] as 


2 tot ss. 
Q-BIF J E = Ey + Ey Ey + Ey + a Uy 
and substituting into Equation [12] the work done by the 


internal forces of the shell becomes 


4 $ * 
Vc), EC-08r + Ete aff ade -EAe,WoLt. [13] 


The strain term — when expanded using Equation [5a,b,c] 


becomes 


Fo =a%o) +3 6lre) + 3x%lee) 


where 


alr/e) = 724) +} 1 (Que 
+ alae). es 


(4,6) “(Beta (See) + de 


and 


* du 


¥(#,0 = “(5x ae 


Substituting Equations [14] and [15] into Equation [13] the 
potential energy for the shell without stiffeners becomes, 


£ 2r 


Ve ‘JJ {2e(-) re +e Jacdedwy 
poaee B+ y)fa+vz4 83° 
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[16] 
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The work done by the internal forces in the stiffener 
represented by the second term on the right hand side of 


ec ima [11] may be rewritten as 


Vee S{SS | ES(1-r) +E Ae, @:|deda - a LA [17] 


me 


The stiffener is assumed to resist only axial and radial 


forces, therefore, for a stiffener, Equation [7a] is expressed 


as 
= Sitewy os 
6, = €, = 4 (se) tt 
or [18] 
7 = xe dur} Dur 4 gt /D'ury* 
G. = gigs! Colom sae (aoa 
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Using Equations [17] and [18] the stiffener work term becomes 


ve" } (\ ah Let “ya. + Meee yaa.] C19] 


+ EAS | 
* se | 


Equation [11] gives the potential energy of the stiffened 
shell and it was obtained from the stress-strain relations [8] 
assuming a Poisson's ratio of 0.5 in the plastic region as 
well as the elastic region. Using the following stress-strain 


relations for plane stress, 


Oz = — (E% + VEy) } os = (éy + DEx), 


i ele tee 
Cxy ~ 2d) 


the combined potential energy expression becomes 
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IE) FGe)” + og ple Me a) 
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iii nai — 


It may be shown that Equation [21] reduces to Equation [11] 


for v = 0.5. Equation [21] can be further rewritten as 


V Py Sot. ae, gee es  Rys Soy [23] 
Lie A oral ZAE (aah +V, Js Te, dPe at 
where 


T -[ [ade'de a} JB an'ae 
ONG va) (2) (2) 6 (BE) A: + (2 )(2)  edath ay’ 


C24) 
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and A is the underlined portion of Equation [21]. For buckling 


the spatial function f becomes 


5 -ke 
F(%/0) =pimi7n’ EG Cot NO ) epee ee 
-t 7-8) 
F(nje) = = pun 7H EO Con (27-6) » TLOL 27, 


where n is given by Reynolds (20) as 


457 C27] 
“td 5 
(SE) Vv 


and D is diameter of the shell. 
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i Here n is not a function of the load intensity or the number 


of stiffeners. For stiffened shells these parameters should 
be included. 

1 : The above results represent a rather complete and real- 
istic first order approach to obtaining the energy absorbed 
within a shell structure and for estimating damage data of 
blast loaded stiffened shell structures. As damage criteria 
such non-dimensional parameters as “2, = and t/2a can be 


used for calculation purposes, when the functional form of w 


expressed as Wof(x',6), is specified. (Note: x' = x/L.) 


2.2.2 External Energy Imparted to Shell 
For short duration pulses an approximate empirical 
engineering approach to generating the blast load intensity 
appears useful for predicting shell deformation in those cases 
where the impulse may be known and the loading less certain. 
Consider the energy flux/unit area of a fixed charge weight 
at a certain stand off distance from the cylinder to be given 


approximately as in Reference (21) in the form, 


W 
Ey ae 


In the above expression, W is the charge weight, R the stand- 


ee ee Tey On ee 


off distance, c a constant, and Es the energy flux. The total 


energy imparted to the cylincrical shell is Et = EgA(2taL) 


Saitaatisiieneeneiatahent annie rete er a ere ees 
2 PT. Se ee 
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where a is the cylinder radius and L the length. Thus, the 
total energy available for damage can be expressed as 

E, ~e# (erat). 
Equating the above expression to the potential energy of the 
shell during deformation as developed in Section 2.2.1 pro- 
vides a means for estimating the maximum deflection once the 
solution of Equation [23] is determined. 

A slight modification of the above development can be 
made if some assumption concerning the type of forcing func- 
tion is made. For example, if an exponential pressure 
variation is assumed, then, we can write for a peak pressure 


P, and decay constant a; P(t) = Poet/a, The impulse/unit 
ao 
area can be easily found from, I = f P(todt. In addition, 


if an expression from Reference (16) is introduced, 


Ry nk Aes iY Pe] dt , where t = time, 


then E¢ can be expressed as E¢ = P28/2p6e5 where 5 represents 
the density of the medium and co the velocity of sound in the 

medium. Once again equating the total energy generated by the 
explosive process to the energy absorbed during the deforma- 


tion of the shell, that is, E,~V, leads to 


V= pL (7a. L) ) (27) 
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where V is given in the preceeding development. The above 


result appears most interesting in that the essence of th 


P-I (Pressure-Impulse) or so-called iso-damage curve is re- 
cognizable in the form of the equation, and graphical displays 
of this quantity can be readily plotted. An extensive amount 


of work in this area has been discussed in Reference (16). 


2.2.3 Radial Deflection of the Shell 

The preceeding discussion relates to the mechanics of 
establishing the energy stored within the shell during the 
deformation process and the external energy imparted to the 
shell generated by an impulse or pressure pulse loading. In 
order to establish any meaningful damage or failure criterion 
some statements regarding the form or criteria for establish- 
ing the form of the displacement functions are necessary. 

As mentioned in Section 2.2.1, the in-plane displacement 
components are generally neglected since they are considered 
small relative to the radial displacement coordinate. For 
determining the latter quantity several approaches are 
available. 

One direct approach for finding w(t) when w(A) and the 
applied external pressure loading are known is to use Hamil- 
ton's Principle for dynamic systems. Recall that this Prin- 


ciple can be written as 
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et 
SJ. (T =) dt a ts ae for time t, [28] 


with U = V - W, where V represents the potential energy and 
W the work done by the external forces. The Kinetic energy 


T can be written as 
e 2 
Tb) pW, [29] 


where yp(A) represents the mass/unit-area of the structure, 

dA the element of surface area of the structure, and 

wo Wo(t)£(A) represents the radial deflection in terms of 
time and spatial dependency. The corresponding quantity U 
represents the difference between the potential energy and 

the work done by the external forces in the (u-v-w) directions. 
The quantity V can be further expanded in terms of a power 


series of the form, 
— > = 2 = s 
V=A+Bar+Cwr + Dus +---- [30] 


and the work done by the external forces as, 

J PAA) FA LA, [31] 
with P(A,t) representing the time varying pressure applied 
to the external surface of the structure. Substituting the 
above in Hamilton's equation results in an equation of the 


Euler form, that is, we can recast Hamilton's Principle as a 
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problem in the Calculus of Variations of the form 


. b 
L - J FC%,Y,4') ax , [32] 


a 


and obtain Euler's equation, 


dF _o oF) so. [33] 


AY. AX VOY 
In the present example, y' = Wo x = t, and y = w(t). 
Upon substitution in the above Euler's equation results in 


) the following nonlinear differential equation for Wo» 
ee 2 - -_ _ 2 
liz ) HAF A dA + (B+ 2luzs + 3D uz +.) 030) 
= [ Pla,t) FA) La. 
with initial conditions, wo(0) = 0, wo(0) = 0, so long as 
the radial deflection continuously increases without unload- 
ing taking place. 


As an alternative to solving the above equation the de- 


flection can be considered to be of the form w = wof(x',@) in 


the post impact range with x' = x/L and L being the shell | 
length. Then for a collapse type behavior of a shell loaded ) 


laterally on the side facing due to a blast load we can write, : 


ar(4j@) = awe /a* —[De)'(/-2y')* , oxn'ck, 
(3s]* 
urxie) = 0 on0/@ (BY (pen) ) -$ex'co, 


& 


*D5 shown in Figure 2. 
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a ae 


while correspondingly for circumferential buckiing the pattern 


is, 
UW (x;Q) = We an Tn'E  corne » OLE LTT, 
; , -€(e7-8) [36] 
uy (%/ 8) =Ws pm 7T¥ EC Cotn (27-8) , LO Ler. 


The above functional forms for w can then be substituted into 
the energy expression for the deforming shell, that is, V and 
equated to the external energy in order to determine Wo as an 


explicit function of V. 


2.3 Numerical Example 

In order to demonstrate the theoretical developments pre- 
sented, calculations for three different L/D (Length/Diameter) 
ratios (1, 2, 4) and varying D/t (Diameter/thickness) ratios 
between 100 to 1200, for each L/D, have been made. The two 
specific buckling energies V,,I, as defined in the text of this 
report have been calculated for both the unstiffened and stiff- 
ened shells and are presented in the accompanying Tables III - 
V. Each of these tables represents data generated for eight 
stiffeners. The influence of increasing the number of stiff- 
eners on the specific buckling energy within the shell is in- 
cluded in Table VI. The last case of Table VI, with 105 stiff- 
eners, represents the maximum packing density of stiffeners 
which could be suitably retained within the outer shell geometry 


and was obtained by taking the width of the stiffeners used in 
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re 


the present calculations and dividing this number into the 
inner circumference of the shell. The specific input data 
and assumptions used in the computations are included below 


in Figure 8. 


0.06, 


ae 


0, #07 


Ww 


Figure 8. Stiffener Dimensions and Assumed Stress- 
Strain Curve. 


The effect of the stiffener is seen by comparing the 


— Ar cee ee SR Ie a 
sachs ksh ie intial nln tinal eda 


values of the VSH1 and VST1 or ISH1 and IST1 of Tables III - 
VI. In all cases the shell terms (VSH1,ISH1) are almost 

two orders of magnitude greater than the stiffener terms 
(VST1,IST1) which points out one of the interesting conse- 
quences of the results obtained, i.e., the shell itself is ; 
much stronger than is the influence of adding reinforcing 


stiffeners. This result occurs even though added stiffness 


$1 


in the boundary attachments of the stiffeners has been included 
in the theory. This is further demonstrated by the fact that 
for the case of maximum stiffener packing the specific buck- 
ling energies are less than those for the case of a shell having 
twice the original thickness. In order to show trends for 
various parameter changes several cases of the non-stiffened 
shells are shown in Figures 9 - 11. Perhaps a more exact re- 
presentation of the stiffener effects, that is, the use of 
anisotropic constitutive equations and the associated solution 
of the resultant coupled governing equations would lead to a 
more significant contribution of the stiffeners in the speci- 
fic buckling energy terms. 

One important result as determined from observation of 
several of the post buckled reinforced shells tested and 
described in Reference (19) is that the addition of stiffeners 


can influence the wave number of the buckled shell. This can 


be summarized by use of the following table. 


TABLE VII. STIFFENER BUCKLING INFLUENCE 
| 


with wave number equal 
to number of stiffeners 
in the shell 


number associated with either number 

of stiffeners or calculated by Reynold's 
eriteria, whichever yields the higher 
number. 
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Figure 195. Unstiffened Cylindrical Shell Parameter Study. . 
vs w/a. Results of Table V. 
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The results obtained thus indicate that while the stiff- 
ener influence in the specific buckling energy is weak at 
best, the stiffeners may play an important role in any failure 
eriteria because of the manner in which they influence the 
circumferential buckling mode. Such a result would be extremely 
important if for example a failure criteria adopted by Shuman, 
based upon empirical data collected through exhaustive testing, 
were used, Reference (4). That is, specifically setting fail- 
ure as some fraction of the half amplitude deflection relative 
to the radius or diameter of the shell specimen tested. In 
Shuman's (4) work this criteria has been set at 5 to 10% of the 


shell radius. oa 
Original 
Shape 
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uo 


Figure 12. Buckled Modes Shapes. 


Figure 12 shows the inner relation between the number of cir- 


cumferential waves and half amplitude deflection wo. Since 
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the change of perimeter during the loading period is negligible 


the number of circumferential waves as well as wo will increase 


with the increase of such pressure or peak impulse whatever 
criteria is selected as defining loading. Such dependency of 
the wave number n on the loading and geometrical parameters 
is not included in either Shuman (4), or Reynolds (20) failure 


criteria. 


SECTION III 
DYNAMIC ELASTIC-PLASTIC RESPONSE OF FINITE-LENGTH 
STIFFENED CYLINDRICAL PANELS TO BLAST LOADS 
3.1 Introduction 

In this section the analytical model (DEPROP) of Mente 
(17, 18) previously described in Section I has been extended 
to include axial stiffeners. The selection of this model for 
modification is based upon its availability and neasonable 
Success in predicting maximum deflections of unstiffened 
cylindrical shells subjected to intermediate blast loads as 
reported in Reference (19). 

For this particular study the stiffeners have been in- 
cluded only for a single layer with the stiffeners having the 
Same material properties as the panel material. The portions 
of the DEPROP code (18) associated with multilayer and honey- 
comb panels have not been modified in this study. The code 
has been modified to include stiffening for a single layer 
with both static and dynamic response for either elastic or 
elastic-plastic stress strain relations. 

The stiffeners are assumed to be placed on the center 
line of the panel thickness. The selection of the stiffener 
size is determined in the following manner. Assume the 
stiffener resists forces only by bending about its transverse 


axis and by tension or compression in the axial direction. 
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The stiffener is assumed not to buckle. The actual stiffener 


and the modelled stiffener are shown in Figure 13. 


A= X-SECT.AREA | m4 
3 
Telyy 
+39 x, 
ts Wi | 
a) actual stiffener b) modified stiffener 


Figure 13. Stiffener Dimensions. 


The coordinate system and deflection are shown in Figure 
6 of Section I. In many cases the same words or phrases as 
used in Reference (18) are used in this report. Full credit 
for the basic model of the unstiffened panel is extended to 
the authors of Reference (18) and any duplication of sections 


of Reference (18) in this report are intended for clarity. 


3.2 Analytical Development 
The governing equations of motion for the stiffened 
cylindrical panel are obtained from the principle of virtual 


work for a dynamic structural system which is given by (22) 


FUSS ap sdyav +$ SSSaise, av 


st -3r - SFsraape=o (37) 


where 


035 stress component 

€ij strain component 

s number of stiffeners 

t time 

T Kinetic energy 

x external force 

V volume 

A area 

ra deflection 

u,V5W deflection in x,y,z direction and 


components of general deflection uj 


and 
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The dotted underlined terms represent the contributions of 


the stiffeners. If it is assumed that T = Tuy; U2; U3), then 
OT ° 

= <= u. 38 

ST = sR, $4, [38] 


with Equation [38] and using integration by parts, 


CTdt = - 46 or) Suid, 
(39] 


The cylindrical coordinates (x,®,z) and the axial, tan- 
gential and radial displacement components (u,v,w) are shown 
in Figure 6 on the coordinate surface which is located at the 
median surface of the panel. It is assumed that p(x,6,t) 
acts on the coordinate surface of the stiffened cylindrical 
panel. As the panel surface deforms, the elemental pressure 
force vector also changes as the element surface area of the 
deformed panel changes. Thus, the vector scalar product of 


the force and virtual displacement is expressed as 
FSF = PH, €)(MnSTy +My 8% +, S14)- [40] 


By linearizing nx, Ny» Nz and recasting 6r,, Sry, Srz 
in terms of the virtual displacements é6u, 6v, dw, the virtual 


work done by the external forces is given by 
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SJF -SFdA = §f P(%,t) (Ni Su+ N, So +N, Sur)LA [41] 


where 


The subscripts on the displacement components denote spatial 
derivatives and ® denotes the initial radial imperfection in 
the panel. Assuming 645 = €;5(u,v.w), ey = Ex (usw), by the 
chain rule, 


5é., = Bos Su + 2b Far +28 Sw, 


~ _ 06 O€é 
SE, = aif Su + Soe bur, 


and substituting these terms and Equations [41, 42] into Equa- 


tion [37] the equation of motion becomes 


” RiEACe i) aes + SMe, Ses dy 
+) $95 Uy Derdy - SS p(%,6,t) Nu] Su 


oe) “SSSeu4 Cae & dv - Sauna da] Sar 
pre Sele) + SSS 05 Pea 


+) S6Se 2 av “ ay) e wat) Nur LA | or 
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The displacement components are assumed in the following 


truncated series form 


MN 
uot) = ) S Uncle) Ser) KO)» 


ar (X,0,4) = y Ahan (€) dw) d.(2), 
wna) ¥ ¥ ume) £0) FO), 


UL (4) ef ute) LO) #0) 
wat) t ie whl) Blot) (0), 


m= 31 


gee ee 


where $m(x) and $,(@) are functions that satisfy the geometric 


boundary conditions of the stiffened cylindrical panel. The 


as 


initial radial imperfection of the stiffened cylindrical panel 


is represented by 


m oN re 
W(H,0) = yy » Ann BH) ©) [46] 


mal rz 


ur (W, @,) = a3 Ann 4 &) $G) 
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where Any are prescribed values based upon assumed deviations 

from the ideal shape of the stiffened cylindrical panel. 
Determining the virtual displacements from Equations 

C45] and substituting into Equation [44] the 3MN equations 


of motion may be determined and are given below. 
5 
1 @T)+ )d (ar ) + SVSo: Ed 
at DLLign =, sa.) a Oden C47a] 


+ F 5950s 38, av eke =O 


4(BE.) + SJSq Sa av - [fem Se da = 0 toms 


5 > 
4G )+ aRE) + [[Sou pi 
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Equations [47] may be simplified by introducing the 


following notations, 


! 
a 
= 
& 

| 


Qn 


oNe =— | C48] 


Onn 


Qs = ° Nu Our ‘ 


mA 


The kinetic energy of a cylindrical panel is given as 
£ £8, 
ae neti) (WU + rr* 4+ Ur*) dnde. [49] 
ovo 


The strain-displacement relations used in this analysis 


are based on the following assumptions (18). 
(i) the strains are small compared with unity 
(ii) the thickness of the shell is small compared 
with the radius and 
(iii) the Kirchhoff-Love hypothesis that straight 
fibers which are normal to the undeformed co- 
ordinate surface remain straight and normal to the 


deformed coordinate surface, and are not elongated. 
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The total strain consists of membrane and bending components 
expressed by the form 


E (m6, zt) =E(4e, 4) + 2K C4,0,4)- [so] 


The membrane elongation and shear strains (e€xx, €99, €,9) on 
the coordinate surface are expressed in terms of the dis- 


Placement components and their spatial derivatives, 


Eyy = Uy tt (uy +l, +z *) + Wy Wy » 

[S1a] 
> ! ee ey we 
, = Uy +f (uy tUy ) + Wy Wy 


C51b] 
m f | 
Ey Mp tale + MG (Warm) + Za (e-wr) 
[51e] 
to Us, + x (Uy We + UZ UZ) - 
a ee Ue Me o a 


For moderate rotations, the change of curvature quantities 


(Kxx» Kgg, Kyg) are given by 


K un = Wry (1+ 22-4 e uy)» 


[52a] 


Sif l if os 
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{52¢] 


The behavior of the stiffened cylindrical panel material 
is treated as elastic-plastic and the deformation theory of 
Plasticity is used. In the deformation theory of plasticity 
the total strain is a function of the state of stress and 
consists of a recoverable elastic component and a nonrecover- 
able plastic component. It is assumed that the material is 
incompressible. Furthermore, it is assumed that the material's 
biaxial state of stress is represented by the bilinear repre- 
sentation shown in Figure 14 in which the strain hardening is 
defined by the slope Ey. This stress-strain representation 
employs the effective stress (a) and effective strain (€) con- 
cept, in which the secant modulus (E,) indicated in Figure 14% 


is defined by 
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Figure 14. Effective Stress-Strain Bilinear Representation. 


ee om = E. £-€. 
Ec. oe ra = mth, (ores) ) [53a] 


and 09, €, are the yield stress and strain, respectively, from 


the material's uniaxial bilinear representation. The effective 
stress and strain, expressed as o = f(oj4) and € = g(ej5) are 


functions of the total stress and strain components, respec- 


tively. For the stiffeners Equation [15a] reduces to 
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The Hencky stress-strain relations for deformation theory 


(23) are used in the plastic region and are given in the follow- 


ing form: 


t 
) (%; > 3 On $5) 


where E is the modulus of elasticity, v is Poisson's ratio 


and 6j5 is the Kronecker delta. For stiffeners, Equation [53] 


reduces to 


% 


E — 
= 4 
E. [54b] 
For use in Equation [47], the stress-strain relations in Equa- 


tions [52] are inverted into the form ojj = f(eij) for the case 


of plane stress (i.e., O22 = Ogz = Oy, = 0) and are given by 


= HOME +Y eee] (Osh = 12), cosa 
s 


ease tr 
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En 7s ee a En = Eco. 


CBee [55b] 


The initiation of yielding for a biaxial state of stress is 


based on the Mises-Hencky yield criterion and given as 


2 at 
f= G+, -GG. +362) [se] 


where o is the equivalent or effective stress and 


Gi, “Uy J 6, = So ) Oe = Oye 


A kinematic hardening model is employed in conjunction with 
the Mises-Hencky yield surface which accounts for the Bauch- 
inger effect when yielding occurs again due to the strain 
reversals during unloading. It is assumed that during plastic 
deformation the yield surface translates as a rigid body in 
stress space. This type of kinematic hardening model, as used 
in this report, is shown in Figure 15. The change in total 3 
stress components from position (i) to position (f) as in- 


dicated in Figure 15 are defined by aj; and, the corresponding 


~ 


change in the total strain components are defined by 8jj so that 
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Figure 15. Kinematic Hardening Model 
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where 
r = the number of elastic unloadings from the yielded 


condition (r = 1,2,...) and initially 


@i,=O > Aye = BY = 0, 


=—-_— 


(i) indicates initiation of yielding or reyielding, 
(f) indicates the final position prior to unloading. 
The yield criterion for the translated yield surface is based 


on the effective stress given by 
G = [Ct -at) -G -ax)(Gq- aX) 


L 
+ (Gz ay + 3 Chie a) C58] 


é = ESS, mae 2 IIE. BP g * Ge Fe) ] | 
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The general stress-strain relations based on the form of 


Equation [55] is given by 
n E nx a ~ ~ 
Oy =O, + aol WE, Bi) +0 Fee” 8%) Sip] ss 
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and for the stiffeners 
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The boundary conditions on the stiffeners are assumed to 
be the same as those imposed on the panel which are described 
in Reference (18) and established as input data. 

The governing equations of motion (47) for elastic-plastic 
deformation are developed further by performing the indicated 


spatial integrations. For convenience the dimensionless 


quantities 
Se ab yee 2 £ Pe 
mee. ) Vas ) Us OL J W = a ) L. TOW } g 2, 
2 PAL . U9 _ && 
R=¥ } RE =KQ ») o= x C 6 = Oo. ) On Ss ’ 


are introduced into the formulations. In Equation [49] ®, must 
be given in radians, however in the input data @, is given in 


degrees. Care must be exercised when using 0, by itself, such 
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as in J above and in the second term of Equation [59a] below, 


it must be converted to radians for computational purposes. 
With this notation and the spatial integration for the kinetic 
energy in Equation [49] performed analytically, the governing 


equations of motions for clamped opposite edges are given as 


2ek' LB. +g ) Vd (a) ¢ Us $, (9.) 
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The spatial integrations described by Equations [61] are 


to be accomplished numerically. For integration through the 
thickness of the panel in the z direction, it is convenient 


to separate the integrand into parts which either are or are 
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not explicitly dependent on the z variable, that is, those 
variables involving membrane strains and bending strains. 
The total strain quantities ej; consist of the membrane and 


bending components given by 
Ey = Exe + 3 Ves C63] 


Therefore, the integrand can be given by 


where 
O& E, 
£™ = Oy <7 . Foe See + Ore dere ) Cexa] 
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mn man n 
for cylinders and 
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+ a Ge a [ésa] 
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for stiffeners. Although Equations (64], [65] are given in 
terms Umn, Similar expressions can be obtained by replacing 
Unn With Van and Win; Up, with Wp, respectively. 
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The Legendre-Gauss quadrature formula (18) was chosen for 
the numerical integration in the z direction where L is the 
number of points selected through the thickness of the panel. 
For integration in the y and 8 directions Simpson's quadrature 
formula was used. The number of spatial points selected in 
the y and 8 directions are given by M and N, respectively, 
where M and N must be odd numbers. By performing the indicated 


numerical integrations, Equations [61] take the following form: 
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where for the w equation 
by = Rp = j/2 for c-c, c-s opposite boundaries 
By = Re = WS2 for s-s opposite boundaries 


and for u and v equations 


hate = Ze, 


and —; are the zeros of the Legendre polynomial PCE). 
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When symmetry is present in both the y and 8 direction, only 
one quarter of the panel need be considered (18). 

The second order differential given by Equations [66] 
are solved numerically in time and the time wise steps are 


based on the central difference equation 


Nous. thal) Rha. Nae [67] 


where X represent the normalized time dependent displacement 
coefficient Umn>s Vmn> Wmn» At is the time increment and k 


denotes the current time step. The At required for calculations 
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with stiffeners included is independent of @ and for n greater 
than m the At calculated by the program for a given panel thick- 
ness will suffice for the case with stiffeners. If instability 
problems arise the time step should be reduced by approximately 
the ratio of panel thickness to stringer height. 

The reaction forces are calculated in the same manner as 
described in Reference (18) and the stiffener reaction forces 
are simply added to the shell reaction forces at the stiffener 


locations. 


3.3 Panel Loading Program Options 

Four transient pressure loading options (identical to those 
of DEPROP (18)) are incorporated into the DEPROSP program to de- 
scribe the pressure function p(x,8,t) which is designated as 
positive in the negative z-direction. The first loading option 
is an analytical representation of the transient pressure on a 
flat panel generated from the detonation of a projectile fired 
into an adjacent fluid medium. The second option provides a 
method for a more arbitrary pressure loading through spatial 
and temporal discretization of the pressure fielc. The third 
option contains a simple uniform pressure distribution over the 
panel with a combination of exponential and triangular time de- 
cay behind a sharp-edged blast front. The fourth option is also 
for a uniform pressure loading, but with an arbitrary temporal 
discretization. These loading options are discussed briefly in 
the following paragraphs. 
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3.3.1 Loading Option 1 

This option describes a transient pressure loading 
p(x,y,t) generated from test data obtained on a flat plate 
subjected to the pressure loading from the detonation of a 
projectile fired into an adjacent fluid medium. As shown in 
Figure 16 the trajectory of the projectile is defined by the 
obliquity angle ¢ and the detonation position is a distance 
Z from the flat panel. This pressure model is based on the 
assumptions that the shock wave shape is spherical, shock 
velocity is constant at 58,800 inch/second and the pressure 
pulse is triangular with a duration of tg. If it is assumed 
that time (sec) initiates when the blast wave reaches the cen- 
ter of the panel, the time of arrival of the blast wave at an 


arbitrary point (x,y) on the panel is given by 


t, = RA [68] 


where R is the distance from the detonation point given by 


(x? + y? + gt;?? ‘ 


The normal pressure p(x,y,t) at this arbitrary point on the 


panel is described as: 
@C%, 4.4) =O , *£<&, 
P(t) 0 , tortrty C69] 
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O(H +) =Gn era) (1+) , tke St Hy 


81 


FE IF ES 


Projectile 


Trajectory Detonation Point 


Fluid Medium 


Flat Panel 


Figure 16. Loading Geometry. 


where 


Z 7°28 


= 484.95 (t + 
Pm ‘t * Seace 


(psi) 


cos a = Z/R 


3.3.2 Loading Option 2 

A non-uniform load can be applied to either a curved or 
flat panel by specifying a discrete pressure-time history for 
an array of selected points covering the surface of the panel. 
The spatial array of points need not coincide with the integra- 
tion grid point (determined by M, N), but must be a regular 
grid in the sense that all points remain in rows and columns 
in the x-@ plane, although the spacing between rows (and 
columns) need not be constant. The spatial grid should also 
span the entire portion of the panel analyzed or the program 
will be forced to linearly extrapolate pressures toward the 
edges. 


The timewise variation is specified at a set of evenly 


spaced times - the spacing, At, is the same for all points. 


However, the time history of each point does not begin until 
time corresponding to a unique delay time has been reached. 
This delay time corresponds to the time of shock arrival and is 
specified on input for each grid point. It is important that 
the first point in the array to be engulfed have a delay time 
equal to zero. 

Pressures at intermediate times and interior spatial 
points are determined by linear interpolation. No attempt 
has been made to estimate shock arrival at interior points, 


thus, the shock wave will tend to be smeared unless a great 


number of grid points are used. For times beyond the last 
time allowed for in the loading, a pressure equal to the last 


value specified at that grid point is used. 


3.3.3 Loading Option 3 
The third load option assumes a uniform distribution over 
the surface of the panel, with simultaneous engulfment. A 
single pressure-time history describes the entire loading se- 
quence. The pressure loading is an analytical representation 
of a combination of triangular and exponential decay, as in- 


dicated in Figure 17. The pressures in the three regions in- 


dicated in Figure 17 are given in analytical form as follows: 


3 t 
pr(t) = py, (1 - = ) Chee" 
et 
prr(t) = po (1 - = ya. * (t" <°t © €) [70] 
Prrr(t) = 0 Ct > tg) : 


It should be noted that the second function is used for time . 
greater than or equal to t'; hence, by specifying t' = 0 the 
special loading cases indicated in Table VIII (step, triangular, 
impulse, exponential) can easily be generated, where I is the 


impulse and At is the integration time interval. 
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Figure 17. Analytical Pressure Time History. 
(See Table VIII) 

3.3.4 Loading Option 4 

Like the previous option, loading option 4 is appropriate 
for uniformly applied loads without consideration of engulfment. 
With this option, discrete values of pressures are specified 
at a set of times, beginning at zero. For other times, linear 
interpolation is used, except after the last time in the table, 


when a pressure equal to the last value is assumed. 
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TABLE VIII. LOADING PARAMETERS FOR SPECIAL CASES 
(See Figure 17) 


tep Triangular 
Function Deca Impulse 
0 0 0 


Po = 


xponential 
Deca 


1.0 x 1049 


3.4 Numerical Example and Discussion 


In order to check the theory included in the program and 

study the effect of stiffening three cases were run using the 

DEPROSP program. A schematic of the cylindrical shell and the 

stiffener locations are shown in Figure 18. The kind of material 

and pressure loading was the same for all three cases and are 

given in Figure 18. The same basic 180° cylindrical shell axially 
stiffened at the 45°, 90° and 135° positions was used for runs 

with 1) zero stiffeners, 2) stiffeners with height equal to shell 
thickness and 3) a stiffener height of .75 in (1.91 cm). The 


stiffener width for cases 2 and 3 was .121 in (.307 em) which for 
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the latter case gives an approximate bending stiffness of a 
~75x.75x.063 in (1.91x1.91x.16 cm) equal leg angle. The pressure 
loading, simulating a plane wave impinging on a cylinder (24), 
was assumed constant in the longitudinal direction and variations 
with 6 and time are shown in Figure 18. The odd mode numbers n 
between 1 and 13 inclusive were used to define the deflection 
circumferentially and the fundamental mode cf m = 1 was used for 
the axial direction. 

The results of these three cases in general show that the 
effect of the stiffener is minimal when compared to the shell 
without stiffeners. This is in agreement with the results of 
Section 2.2.3 which showed the same general results for the 
modified Greenspon solution of a 360° cylindrical shell. 

For the DEPROSP calculation of the cylindrical panel shown 
in Figure 18, the difference between calculated deflections for 
the stiffened and unstiffened case is so small that if presented 
in graphical form they would be indistinguishable. However, the 
results of the stiffened case are shown in Figure 19 along with 
deformed shape of a similiar cylinder tested under the loading 
conditions specified in Figure 18. The calculated deflection 
shape of Figure 19 shows no local influence of the stiffeners, 
whereas the local influence of the stiffeners is quite discernible 


for the experimental curve. 
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The influence of stiffeners on the stress and strain dis- 
tributions of the unstiffened panel studied (Figure 18) is also 
minimal with only slight reduction in stress and strain. Re- 
distribution of stresses and strains occur but the overall 


changes are small. 
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Simply supported 
Sides and ends 


1 


2.0 in 


Material: 6061-T6 Aluenun 


Yield stress 0 psi 
Elastic modulus 9. 8x 10 psi 
Work hardening slope 60,000 psi 
Density 0. i 1b/in$ 
Poisson's ratio 0.3 


Pressure Loadin 
The distributed load In psi and time in millisec are listed 
below for the given angle. 


0° 22°25 45 : 


280 - .0 260 - .07 170 - .0 106 -. ae 

175 - .14 154 - .21 104 - .23 59 - .26 29 - .27 

102 - .28 86 - .35 60 ~ .37 31 - .40 12 - .41 
§3 - .42 42 - .49 30 - .S1 14 - .S4 5 - .55 
21 - .66 16 - .73 12 - .75 5 - .75 1 - .79 
0 - .70 0 - .77 0- .79 C - .82 0 - .83 


Figure 18. Schematic of Shell, Material Properties and 
Pressure Loading for Numerical Example. 
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Figure 19. 
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Experimental shape 
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SECTION IV 
CONCLUDING REMARKS i 


The emphasis of this study was specifically directed to- 
ward investigating the large plastic deformations of stiffened 
cylindrical shells subjected to blast loadings, that is, as 
associated with vulnerability studies. Many of the previous 
analyses and studies have not been considered in this report 
for reasons stated in previous sections. However, Jones (26) 


has recently completed a literature review of the dynamic 


Plastic response of structures with many studies referenced 

in his review. The references contained in this study and 
those of Reference (26) represent a rather complete biblio- 
graphy of both analytical and experimental plastic response of 
cylindrical shells. 

The primary results of this study indicate that the effect 
of axially stiffening a cylindrical shell using stiffeners 
typical of those in aerospace applications is very small. This 
conclusion is based upon an analytical approach which incor- 


porates the additional stiffness of the stiffeners into the 


energy equations directly by simply adding the reinforcing 

element to the potential and kinetic energy terms of the basic 
shell equation. This type of analysis thus eliminates bending- 
membrane coupling which would exist if anisotropic constitutive 


equations were introduced. 
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This conclusion has been verified by using two different 
methods for studying responses of the blast loaded shells and 
agrees with experimental tests conducted by the USAF Armament 
Lab, Eglin AFB, Florida. Both of the methods used have been 
incorporated into computer algorithms which allow an investigator 
to determine failure modes of blast loaded shells by either an 
engineering approach or a more sophisticated mathematical 
approach. 

These data can then be used as first order approximations 
to examining failure regions at critical shell sections and 


attachment points. 
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APPENDIX A 
PROGRAM FOR GREENSPON SHELL THEORY WITH STIFFENERS 


A-1 Introduction 

The purpose of this program is to evaluate the total 
potential energy of shell configurations with/without stiff- 
eners (stringers). The potential energy of the stiffened shell 
is divided into four parts, there being Vish and Tish for the 
shell and V,,, and I,,;4 for the stringers. These are given by 


the following equations: 


os 
Vek = Aw) A (v1 8) dedy’, (71] 


6.5 ATT , 
[un 4), § [Reve +S 6H e]dode’, [723 


eee resent 


~2: : ; * ie Ds * ie 
B, =@ Fcosne)am(te) » @, = : [75] 
an (i-q) , a0 


in which le is the independent variable. 
The procedure for evaluating the above energy terms is 
included in the user's manual, Section A.2.1, with a listing 


of the program variables used in Section A.2.3. 


A.2 User's Manual 
A.2.1 The formal procedure used to evaluate the two double 
integral types defined by Vyg, and TI,,, in equations [71] and 
[72] is given by the following steps: 

(a) First evaluate the inner integrals keeping x‘' constant, 
that is, 


Yo 


§. Alea) is Dd Hi A(*;@;) = F,C¢), [76] 


j=! 
¥ = = 
5 [a(v,e) +4 @ (w,0) | dedy’ 
£77] 


No 
= du fa@ai+se@iad] = Rw) - 


Flo 
(b) Next, evaluate the outer integrals in a similar manner. 


We thus have 


) 


5 Ny! 
Vek = at F (¥') dy He Fy) wei 
° AZ 
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where Ye and Ny' are the number of steps choosen in @ and ye’ 


coordinates respectively. Hi and Hj are given by the combina- 


tion of Simpson's rule and Newton's 3/8 rule, as contained in 
IBM System 1360 Scientific Subroutine Package (360 A-CM-03X) 
Version II (H20-0205-2), page 88. 


A.2.2 Program Description 
A.2.2.1 Usage 

The program consists of a main program and three sub- 
routines. In main program, two major steps are necessary: 
(1) to evaluate V,,} and I,,, by using the numerical integra- 
tion scheme previously described; (2) to evaluate Vise and 
Tist+ If the input number of stringers is zero, then the 
second step is eliminated (see, for example the flow chart 
description of Figure 20). All input and output is included 
in the main program and all subroutines are also called from 


the main program. 


A.2.2.2 Subroutine Required 


DEL(DELTA,ALBT) evaluates the integrals given by equations 


71] and [72] at the nodal points. 
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WA=WAMIN 
K=0 


DO WHILE 
WA 
WAMAX 


Cee ee a 


Figure 20. Flow Diagram for Greenspon Theory with Stiffeners. 
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DEL 
(DELTA, ALBD 


_ QSF 
(Visn»Tisn) 


SH1(K) =Vish 
ISH1(K) “Tis 


WA= 
WA+DWA 


bs A 


QSF(DX,FY,FZ,NX) evaluates the integrals using a single 
integration. 
scientific subroutine QSF which has been referred to previously. 


ESTRI(VST,EIST,NST,ARAT,ZAL) computes equations [73] and 


C74]. 


A.2.2.3 Program Variables 


A 
AK 
ALBT 


DWA 


ETSH1L 
Ei ST1l 
FY 


Re 


This subroutine is modified from the IBM supplied 


D/2 

K Out of roundness parameter, Equation 26 
Matrix of &(X;,Y;) +96 (%,Y;) 

in equation [70] where i = 1, ..., NX and 
fF = Ls «sey NY 

N 

Area of stringer 


Diameter of Shell 


Matrix of A(Xj,Y;) in equation [71] where 
4 dey ae TX ered Ff oe dp teeny 

D/T 

Step size of WA 

Step size of X 

Step size of Y 

Vector of Ty>, 

Vector of Tist 


Input vector of QSF 
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Output value of QSF 
Length of shell 
Integer defined by (WMAX-WMIN)/DWA + 1 
Number of stringers 
Number of increments of X 
Number of increments of Y 
L/D 
Thickness of shell 
Poisson's ratio 
Vector of Vy>, 
Vector of Vist 
WO/A 
Vector of WA 
Upper bound of WA 
Lower bound of WA 
Radial displacement of Shell 
x! 
Y Radial displacement of shell 
ZBAR Centroid of area 


A.2.2.4 Dimension Requirements 


The required dimension of the variables in the main pro- 


gram are, DELTA(NX,NY), ALBT(NX,NY), FY(NX), WA1(NP), VST1(NP), 
EISH1(NP) and EISTI(NP) where the parameters NX, NY and NP are 


definec in the description of variables used. 


A.2.2.5 Input Formats 


CARD TYPE FORMAT CONTENTS 


1 315 NX, NY, NST 

2 8F10.5 WAMIN, WAMAX, DWA, U, AK, T, AREA, 
ZBAR 

3 2F10.5 SLD, DT 


(CARD TYPE 3 may be repeated as may be 
necessary ) 


Final card Blank 


A.2.2.6 Output 

(1) If NST # 0, the vectors WA1(K), VSH1(K), VST1(K), 
EISH1(K) and EIST1(K) are printed in columns. 

(2) If NST = 0, the vectors WA1(K), VSH1(K), EISH1(K) 


are printed. 


A.2.2.7 Summary of Users Requirements 
(1) Determine values for NX, NY and NST. If NST = 0, 
then the values of T, AREA and ZBAR are neglected. 


(2) Adjust the DIMENSION statements in main program and 


subroutines. 


A.2.3 Test Problem 

The following parameters were used to test the program: 

NX = 26, NY = 26, NST = 8, WAMIN = .0125, WAMAX = .8, 
DWA = .05, U = .5, AK = .25, T = .094%5, AREA = .063 and ZBAR = .2. 
The listing and output for this problem are contained in the 


following section. 
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A.3 Program Listing and Example Output 


Cc Cc 
Cc Ll - LENGTH OF SHELL Cc 
Cc Oo - OLAMETER OF SHELL Cc 
i Cc A - 072 Cc 
‘ Cc T — TRICKNESS OF SHELL Cc 
j Cc OT ~ OsT Cc 
; Cc sto - LsDd Cc 
Cc x - xe Cc 
c Y — RADIAL DISPLACEMENT OF SHELL Cc 
4 Cc wo - RAOIAL DISPLACEMENT OF SHELL c 
| Cc WA - wO/sA c 
i Cc WMAX = UPPER SOUND OF WA Cc 
i Cc wMIN —- LOWER BOUND OF WA c 
Cc OwA - STEP SIZE OF WA Cc 
Cc AREA — AREA OF STRINGER c 
Cc ZBAR — CENTROID OF AREA Cc 
Cc U - PCISSON®S RATIO Cc 
c AK - K Cc 
Cc AN - N Cc 
Cc NST - NUMBER CF STRINGERS c 
Cc NXeNY- NUMBER OF INCREMENTS OF X AND Y al 
Cc 
Cc 
Cc 


IMPLICIT REAL*8 (A-He0-Z) 
DIMENSICN DEL TA(26 026) oF Y( 26) 5 


w 1 (50) -eVST1¢50) 
OIMENSICN ppp hae thay tapes ta at 
J 


~z= 


COMMON UsPI sAKeANecAL oTAe WAcDX 

COMMON 8A2.WAS » WAS cAL2 sAL3 cD AL4S i 

READ LeNXeNY¥eANST | 
1 FORMATC( 31S) } 

READ 3eAMINes WAMAX 0D WA oU eAKeAREAe Ts ZBAR 
3 FORMAT(8F10¢e5S) 

P1=32141593 


INPUT SLC»OT 


140 READ SeStOe0T 
5 FORMAT(C32F10e5) 

IF COTeLE2O) GG TO 150 
AL=1e/SLD/26 
AL 2=AL*#AL 
AL3=AL2% AL 
AL4=AL37AL 
TA=1e/0T 
TA2Z=TASTA 
AN=1657%2.* AL *®DSOQRT(10615/22/AL/DSORT(TA)D-1e) 
N=AN 
TF ( AN—Ne GEo 0 5 JN=NF1 
AN=N 
WA=WAMIAN 
DxX=Le/(NX—-1) 
OY=2e*PI/S(NY-1) 
OxX=O0X/2e6 
DY=DY/26e 
K=0 

130 K=K+1 
WAT (K)=RA 
WA2Z=WAPWA 
WAZ=WAZ AWA 
WwA4=WA3* WA 


EVALUATE DELTA AND ALPHA#BETA/3 AT NODAL POINTS ANO STORE IN 
DELTACI+J) AND ALBT(IeJ) RESPECTIVELY 


CALL DEL(DELTAsALBT) 


SS 


aan 


anna 
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Cc 
Cc EVALUATE ViSH 
c 


OO 100 [=1+sNx 
OO 110 J=1eNY 
110 FYCJI=DELTACI.J) 
} CALL QSF (DY sFY sFZsNY) 
. 100 DELTAC1I»1 )=FZ 
: 00 120 I=1LeNX 
i {20 FY(C E)=SDELTA(T 21) 
CALL OSF (DX sFYsFZeNX) 
VSH1(K)=42*FZ 


co me ht “heb On nh 


Speman ide inten 


EVALUATE 11SH 


00 190 I=1+Nx 
DO 200 J=ieNY 
200 FY( JI=ALBT(I2S) 
CALL QSF (DY »FYeF ZeNY ) 
190 ALBT(1el )=FZ 
OO 210 I=1eNx 
210 FYCII=ALBT( 101) 
CALL QSF (DX oF YeFZeNX) 
EISHI(K )=42 *FZ 
IFQNSTeL E29) GO TO 170 


EVALUATE viST AND I1ST 


ARAT=AREA/( T#¥T#OT/2. ) 
ZAL=ZBAR/(SLD*DT *T) 
CALL ESTRIC VST s-EIST es NSTsARAT»ZAL)D 
} vST1¢(K )=VST 
EISTi(K I=EIST 
170 WA=RA4tDnA 
IF ( WAeLE eWAMAX) GO TO 130 
NP=K 
PRINT 3C 
30 FORMAT(1H1) 
PRINT 162SLDeDTsAKsANeU 
10 FORMAT (4/72 o1 7X0 PL/D=* sEL5S eS eS ol TX s *D/T=H GELS: Sef sli 7Xo*K=* E1565 
Cc oFeolTKXe *N=* cE15 059/01 7Xo *NU=" E1505) 
IFCNSTeLE2QO) GG TO 180 
PRINT 12eARATeZALeNST 
=. FORMAT (1 7X0 *AREASA/ST=* s9E15 05001 7Xs * ZBARSA/ST=" 20E1S Sete 
17Xe*NUMBER OF STRINGERS=° 515) 
<eRINT 2¢ 
| 25 FORE eee er eethiciel ett ee i 
¢c ®ISH1* eo l3xXe*ISTi ts H 
0O 160 I=1»sNP \ 
1 CI) oVSTICIVSETSHIC TI) »9EISTICI) : 
e5)) 


nan 


nan 


160 PRINT 15,WA1(1)sVSH 

1S FORMAT (15X%e5(2xX2E15 
GO TO 140 

1680 PRIAWT 40 

40 FORMAT (/ o28X0°WOSA® 51 3X0 *VSH1® o13Xe ISHII s/) 
00 220 I=1.NF 

220 PRINT 3£2WA1( 1) e VSHI (I) SEISHICI)D 

3S FORMATC15Xe 3( 2X0E15e5)) 
GO TO 140 

150 PRINT 30 
STOP 
END 
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110 


120 
130 
100 


SUBROUTINE DELCOELTAeALBT) 

IMPLICIT REAL#S (A-HeO~Z) 

OTMENSICN DEL TAC 26226) sALBT( 26226) 

COMMON UsPI sAK seANoAL oT Ae WAsDXeDY eNXeNY 

COMMON BA2e WA2 eWAS PAL2Z cALZeALS STA2Zs T 

AS INH(X ) =OLOG (X*#DSORT(C 1 +X¥X)) 

X X=—DX 

OO 100 I=1eNX 

XR=AX+OX 

Y Y=-DY 

DO 100 yv=1eNY 

YVY=YV*DY 

yY=vYY 

X=XX 

IF CVVeGTePI )Y¥=2-ePI-YY 

SPXZDOSIAN(PI®X) 

CPX=DCOS (PI*#X) 

EK Y=DEXP (-AK®Y) 

SNY=DSINCAN#Y ) 

CNY=DCOS(AN*Y) 

F=SPX*EK Y* CNY 

F X=PLI SCR X*E KY*CNY 

FY=—AN®SPX*EK Y#SNY—-AK*SPX*#EKY*ECNY 
FXXs=PLAPI* SP XE KYSCAY 

F YV=2e *AN*®AK* SP X*E KY SSNY+ ( AK #AK—AN#®AN) *SPXBEKY*CNY 
FXY=—AN*PI #CPX*EKY €SNY-AK#PI *CPXREKYECNY 

IF CVVeLEePIIGO TO 1190 

FY=-FY 

FXY=—FXY 

F 2=FeF 

FRX2=FXFX 

FY2=FY*FY 

F X4=F X24F X2 

FY4=F Y24FY2 

ALPRA=WA4S FAL4 FF X4/4 0 CWAS FALZ*¥FXZRFY 2/20 -U*WAS FAL 24F SF XD 
CWA SSF Y4/4e— BASZ*FUFY 24 WA2S*F2 

GAMMA=— WAS#AL4 CTASFE X2*FXX—WA SETAE YY €F Y242 0 FWAZRTAS FY YRF 
C-US WASP AL2* TA FFX 24FY Y-U WAZ FAL 2% TAR FXXRFE VY 2425 CU*SWAZBAL2ZETA 
CRE REXX—4 0 (LeU) Ze  BASFAL ZETA KF XY SF REF Y 

BET A=WA2Z FAL 44 TAZ RE XX FF XX 42 oe FUP WASKTAZKAL 2HF XX*FYY4 WAL 
CTA2Z SF VYRF VV O20% (19 —-U DE WAZRAL2AKTASRFXYVRFXY 

DEL 1=( 2 « #BETAtGAMMA) *DSQRT (ALPHAt*GAMMAtTBETA)/42a/BETA 
DEL3=(—-2 oe #BETAtGAMMA )*DS ORT ( AL PHA-GAMMAtBETA)/42 /BETA 
ABG=4.¢ * ALPHA* BE TA~GAMMA* GAMMA 

ABGLM=e 1E-12 

TF ( ABGeGTe-ABGLMe AND ec ABGeLToABGLM) GO TO 120 

DEL 22 (4,4 * AL PHA*BETA— GAMMA® *2 )/8,/BETA/DSORT (BETA) ®ASINH 
ClC2e*BE TA*+GAMMA D/DSQRT (46 *® ALPHA®BE TA—GAMMA ® 2 ) } 

DEL 4=( 4. * AL PHA*BE T A~GAMMA# #2 )/8e/B8ET A/DSQRT (BETA) AS INH 
CCC  2e8BE TAtGANMA )/DSOART( 42 *ALPHA*BE TA- GAMMA #2 D ) 

GO TO 130 

DEL 2=0 

DEL4=0.6 

DEL TAC 1 e JI=DEL1I ¢DEL2-DEL 3-DEL4 

ALBT(I » J} =ALPHASBETA/3 ¢ 

CONTINUE 

RETURN 

END 
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SUBROUTINE QSF(HsYsZ »sNOIM) 
IMPLICIT REAL*8 (A-HeO-Z) 


OIMENSICN Y¥(26) 
COMMON UsPT sAKs ANec AL oT An WA 2D Xe DV eNXe NY 
COMMON WA2Z ce WAI WAS 2ALZe ALI 2 ALS eo TAZ0T 


HT= 03333333 *h 
TF CNDIM—-S)7 07 ol 


NOIM IS GREATER THAN 5e PREPARATIONS OF INTEGRATION LOOP 
SUMI=Y(2)4#+Y(2) 

SUM1=SUML+SUM1 

SUMI=HT* (V (1) #¢SUM1 47 (3)) 

AUX1I=V¥(4)+4Y(4) 

AUX 1=AUX1 FAUX1 

AUX1=SUM1 4HT* (CY (3) ¢AUX1+Y(5)) 

AUX2=HT# (¥(1) 430 A75#(V(2 DF Y¥(5) ) 420625" (Y (3) 4¥(4) 4069) 
SUM2=Y(5 )#+¥(5) 

SUM22SU¥24SUM2 

SUM2=AVU X2-HT *(Y¥(4)4+SUM24Y(6)) 

AUX=¥(3)t+Y(3) 

AUX=AUXt# AUX 

IFCNDIM-6 152502 


INTEGRATION LCGP 
DO 4 I[=7eNDIMe2 
SUM1I=AUX1 
SUM2=AUX2 
AUX1=Y( I-1) 
AUXIL=AUX1+4A 
AUX1=SUM1 4H 
IFC I-NOIM)3 
AUX2=Y( 1) FY 
AUX 2=AUX24AU 
AUX2=SUPM244T#(Y CI=- L)+AUX2evC 1e1)) 
CONTINUE 

Z=AUX2 

PETURN 

Z=AUX1 

RETURN 

END 


-1) 


*Y¥CI- 

UxX1 

T#e(¥ (1-2) ¢AUXI1FV(I)) 
20606 

(1) 

X2 
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SUBROUTINE ESTRIC VST sEISTeNSTcARAs ZL) 
IMPLICIT REAL*8 (A-HeO-Z) 

COMMON U oP I eAK sANoAL oTAe WA20X DY oNXo NY 
COMMON BA2e WAI eo WAS eAL2sAL3vAL4 go TA2ZeT 
vsT¥=0 

EItst=0. 

NY=nNST 

00 10 I=1.N1 

J=I-1 

Y=2e8PI*#J/N1 

TEC V¥eGTePI) Y=2.*PI-Y 

PHI=DEXP (-AK*#Y)*#DCOS (AN*®Y) 

PHY 2=PHI*PHI 

PHI 3=PHI*PHI2 

PHI42PHI2*PHI 2 

PI2=PI#*¥FI 

PIZ3=PI*PI2 

PY4=P12ePI2 
VST1i=107320508*WA2Z/2 oe FAL Z2¥PI2*PHI2/2¢*ARA 
VST 2=1e 732050 &* WA*AL #20*PI*PHI SZL#ARA 
vsT=vsT+t+vsTi¢evst2 

UU1L=1 e-UFU 

EISTI=VUU 1 #e 25 FWAS RAL 243, /8 2 PISA *PHIGSARA 
EIST2Z=UU1 *WAZ*AL3*22. /30*PI S*PHIS*ZL*ARA 

E IST3=UU1 FWA2ZFAL2€ eS #PI4*€PHI 2*( ZL+TAKAL ) F¥2KARA 
EIST=EISTHEISTI*FEIST2tEISTS 

CONTINUE 

RETURN 

END 
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APPENDIX B 
PROGRAM FOR MENTE DEPROP CODE WITH STIFFENERS (DEPROSP) 


B.1 Introduction 

The modified DEPROP code will be called DEPROSP, acronym 
for Dynamic Elastic Plastic Response of Stiffened Panels. 
The DEPROSP program has all the capabilities of DEPROP plus 
the capability of handling axially stiffened flat and cylin- 
drical panels for a single layer of isotropic material. The 
stiffeners (stringers) must be of the same material as the shell 
and have the same boundary conditions. All the provisions of 
the DEPROP program as pertaining to either static or dynamic 
loadings remains the same. The words stiffener and stringer 
as used interchangably here and in the text have the same 
meaning, i.e., the axial stiffening element added to panel. 


In Section B.2 the subroutines are listed and described, 


flow diagrams are shown, and the major program variables are 


listed. Program input data is given in Section B.3 and program 
output and error messages are described in Section B.4. A 


complete program listing is given in Section B.5. 


B.2 Description of Subroutines, Flow Diagrams, and Program 
Variables 
B.2.1 Subroutines and Flow Diagrams 

Table IX lists the 19 routines and all commen blocks 
which make up DEPROSP. The decimal length of each common is 


also indicated. 
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Flow diagrems of the majcr routines are presented in 
Figures 21 to 26, while brief descriptions of the purpose of 
all the routines are given below. Also included are lists 
of the routines which are referenced by other routines and 


vice-versa. 


DEPROSP 


Main program. Reads preliminary input data and controls pro- 
gram flow. Calls PINIT, PROP. 


BOLT 


Sets up W mode shapes for boundary conditions selected. 
Called by DSET3. 


DERV2 

Computes strains, displacements, and accelerations in the 
main integration loop. 

Calls LIST1, LIST2, PRESS, REIT, SIGMA, SOLVE. 

Called by PROP. 

DSETL 


Reads DEPROSP input data and calculates constants. 
Called by PROP. 


DSET2 

Calculates constants used in DEPROSP. 

Calls LEGEND, DTSTEP. 

Called by PROP. 

DSET3 

Calculates additional constants and writes out a description 
of input data. 


Calls BOLT. 
Called by PROP. 
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RRS 


DTSTEP 
Computes an integration time step small enough to avoid 
numerical instabilities. 

Called by DSET2. 

HIM 


Numerical timewise integration routine. 
Called by PROP. 


LEGEND 

Sets up constants for Gaussian integration through the thick- 
ness for an elastic-plastic solution. 

Called by DSET2. 

LIST1 


Output routine for the elastic-only option. 
Called by PROP, DERV2. 


LIST2 


Output routine for the elastic-plastic option. 
Called by PROP, DERV2. ‘ 


Reads in pressure loading time history. Sets up loading 
functions. 
Called by PROP. 


PINIT | 


PRESS | 


Calculates pressure loading at a given time. 
Called by DERV2. 


PROP 


Executes the static and dynamic panel solutions. 
Calls DERV2, DSET1, DSET2, DSET3, HIM, LIST1, LIST2, RELAXP, SEC. 


Called by DEPROSP. i] 
REIT | 
Computes reaction forces at the boundaries and corners. 
CALLED BY DERV2. 
| 
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RELAXP 


Solves simultaneous nonlinear equations representing preblast 
conditions using a relaxation procedure. 
Calls SOLVE. 
Called by PROP. 
SEC 
Finds elapsed CP time. 
Calls system routine SECOND. 
Called by PROP. 
SIGMA 
Computes stresses in the main integration loop for the elastic- 
Plastic option. Is not needed for the elastic-only option. 
Called by DERV2. 
SOLVE 
Solves a set of simultaneous linear algebraic equations. 
Called by DERV2, RELAXP. 
B.2.2 Major Program Variables 

The major program variables are defined in this subsection. 
These variables are listed alphabetically with a brief descrip- 
tion devoted to each one. An asterisk preceding a variable 
name indicates that the variable is input as run data. The 
dimension of a variable is given parenthetically after the 
variable name, where a numerical dimension indicates the fixed 
amount of storage required for the variable. There is no need 
to change the dimension of such a variable. However, other 
variables have variable dimensions. As an example, one of the 


first dimensioned variables listed is BETR(NBAR). The dimen- 


sion is a variable, NBAR, which represents the number of spatial 
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integration points used in the beta-direction in the panel. 
This dimension must be the largest number of points the user 
intends to employ. In the DEPROSP program, this dimension is 
23. If additional integration points are required, the dimen- 
sion of NBAR must be increased, thus increasing the dimensions 
for all variables with the dimension, NBAR. The current 
dimensions provided for the variables in the DEPROSP routines 
are given in Table X. 

Almost all of the variables which may require dimension 
changes as indicated above are contained in the COMMON blocks. 
There are a few exceptions and, in such cases, the subroutine 
in which the variable is dimensioned is indicated in the list 
of variables or in Table XI. If the dimensions are changed, 
certain additional changes in the program may be required. 
These changes are also indicated in Table XI. 

Many of the variables found in the program result from 
their use in the larger NOVA-2 code (Reference 25). In most 
eases, these variables have little or no use in DEPROSP and 


are so indicated. 
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DEPROSP 


(SETUP) (STATIC BRANCH) 


BOLT 
DSET1 
DSET2 
DSET3 
STSTEP 
LEGEND 


(DYNAMIC 
LOOP) 


LIST1 
LIST2 
SOLVE 


Figure 21. Major Program Flow (DEPROSP). 
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Figure 22. Program DEPROSP Flow Diagram. 
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Figure 23. Subroutine PROP Flow Diagram. 
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Figure 24, Subroutine DERV2 Flow Diagram. 
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Figure 25. Subroutine SIGMA Flow Diagram, Part 1. 
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Figure 26. Subroutine SIGMA Flow Diagram, Part 2. 


TABLE X. DIMENSIONS OF DEPROSP VARIABLES 


LBAR 
LBARST 
MB 
MBAR 


MBAR®NBAR*LBAR+ 


MG 

MGMB2 

NBAR 

NBNBTS 

NKP 

NL 

NPX 

NPY 

NSTR 

NSTR*NBAR 

NSTR*NBAR*LBARST 
This constraint is only significant for an elastic-plastic 
run (NDERV=2). Three possible combinations using maximum 
dimensions are: (17x17x6), (19x19x5), (21x21x4). 
The total number of modes selected (MGMB) from the total 
possible (MB*MG) cannot exceec 49. 


The total number of spatial integration pointsallowed 
(NGNBT) from the total possible (MBAR*NBAR) is 361. 
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B.2.3 List of Program Variables 


A Radius of the cylindrical panel, a, in. Is 
set equal to 1.0 for flat panel 

*AA Pressure loading parameter, a, dimensionless 
(Load Option 3) 

\ 

ALTT (LMAX ) Change in stress component, agg, psi 

ALXT (LMAX ) Change in stress component, tg. psi 

ALXX (LMAX ) Change in stress component, . psi 

ALXXST (LMAXST ) Change in stringer stress component, 
xx» PS 

*ANN Pressure loading parameter, n, dimensionless 


(Load Option 3) 
AZ Constant equal to a/ty, l/sec (Load Option 3) 


AAU (MGMB , MGMB ) Matrix of U,, coefficients, making up part 


of the set of simultaneous linear algebraic 
equations used in stringer calculations, 
1/sec2 


AAW (MGMB ,MGMB) Matrix of W,, coefficients, making up part 
of the set of simultaneous linear algebraic 
equations used in stringer calculations, 
1/sec2 


BBU (MGMB ) Right hand side vector of computed constants 
making up part of the set of simultaneous 
linear algebraig equations for stringer 
calculation of Umn, 1/sec 


BBW(MGMB ) Right hand side vector of computed constants i 
making up part of the set of simultaneous 
linear algebraig equations for stringer 
calculation of Wn, 1/sec? 


BE1(LMAX) Change in strain component, on in/in. 
BE2 (LMAX ) Change in strain component, Fess in/in. | 
| 
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BE3 (LMAX ) 


BE4 (LMAXST) 


BETR(NBAR) 


BTL(NL) 


*BSTR 
BXL(NL) 


BXLST(NL) 


CCRIT(NL) 
CC1(MG) 
CC2(MB) 
CC5(MG) 
CC& (MB) 
CINST(3) 
CK(6) 


CM11, CM12, 


CM22, CM33 


CM11ST 


CN10 


% 
Change in strain component, 8,g, in/in. 


Change in stringer strain component, 
XX in/in. 


Integration positions in the beta-direction, 
rad 


Constants used in elastic option, EK/(1-vkv§), 


psi 


Width of stringer, bs» in. 


psi 


Constants used in elastic option for stringers, 


EX, psi 

Not used 

Constant, m, dimensionless 
Constant, n, dimensionless 
Constant, mtl1, dimensionless 
Constant, ntl, dimensionless 
Not used 


Constants, 1/ky >» 1/kg; for the equations of 
motion, dimensionless 


Stiffness constants Cjj used in elastic, 
multi-layer integration through the thick- 
ness, lb/in. 


Stiffness constant Cj] used in elastic, 
single layer integration through the 
stringer thickness, lb/in. 


Constant, 2L?R for elastic option, 2L2 for 
elastic-plastic option, dimensionless 
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Constants used in elastic option, ES/C-vEve), 


| 
| 


CN1l 


CN6 


CN7 


CN8 
CNS 
COSB(MB*NBAR) 


COSG (MG*MBAR) 


COS2B(MB*NBAR) 


COS2G (MG*MBAR ) 


CRIT(3) 
DC 
DELT 


*DELTIM 

DELX 

DELX (3*MGMB) 
*DET(NPY ,NFY) 


DM11, DM12, 
DM22, DM33 


Constant, 2L*R for elastic option, L?/6R2 
for elastic-plastic option, dimensionless 


Constant for elastic-plastic option, E/(1-v?), 
psi 


Shear modulus for crate ini option, 
E/2(1+v), psi 


Constant, L?, dimensionless 
Constant, L2/2R, dimensionless 


Cosine functions of the beta angles, 
cos ((nt+1)85), dimensionless 


Cosine functions of the gamma angles, 
cos ((mtl)yjz), dimensionless 


Cosine functions of the beta angles, 
eos (n 85); dimensionless 


Cosine functions of the gamma angles, 
cos (m yj), dimensionless 


Not used 
Not used 


Grid point spacing in B-direction, Ay, in., 
or aA®8, in-rad 


Integration time interval, sec 

Grid point spacing in y-direction, Ax, in. 
Working array in RELAXP, dimensionless 
Delay time to when spatial point is first 
engulfed by pressure wave, sec (Load 


Option 2) 


Stiffness constants Djj used in elastic, 


multilayer integration through the thickness, © 


1lb/in. 
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DM118ST 


DPRT 


DPRT1 
*DTIM 


DWB(NGNBT ) 


DWG (NGNBT ) 


DWO(NGNBT ) 


DX1 (MBAR ) 


DY1(NBAR) 
EC 
EL 


#EM(NL) 


ENX (2*NBAR-2) 


ENT (2*#MBAR-2) 


*EP 


EPB(LMAX ) 


Stiffness constant Dj] used in elastic, 


Single layer integration through the 
stringer thickness, lb/in. 


Running time, in multiples of DPRT1, used to 
flag next printout, sec 


Time interval between printouts, sec 


Time interval between specified pressure 
data, sec (Load Option 2) 


Values for imperfection-related partial 
derivative Wes dimensionless 


Values for imperfection-related partial 
derivatives Wy, dimensionless 


Values for imperfection-related displacement 
a, dimensionless 


Fractional distance in y-direction, locating 
grid point between two pressure-mesh points, 
dimensionless (Load Option 2) 

Same as DX1l, only in §-direction 


Not used 


Modulus of elasticity for elastic-plastic 
option, £, pei 


Modulus of elasticity for each layer for 
elastic-plastic option, E", psi 


Tangential reaction force per unit length 
along boundary, N,,, 1b/in. 


Tangential reaction force per unit length 
along boundary, Ng, lb/in. 


Strain hardening slope for elastic-plastic 
option, Ey, psi 


Temporary storage used for either eorg 
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EPBO(LMAX ) 


EPBOST (LMAXST ) 


EPBSTR(LMAXST ) 


EPO 


*EPSIF 


ERR(3*MGMB ) 


*ET (NL) 


ETT 
ETT1(LMAX ) 


*EX (NL) 


EXT 
EXT1(LMAX ) 


EXX 
EXXSTR 
EXXST1(LMAXST) 


EXX1 (LMAX ) 


Equivalent stress, squared, when response 
is still elastic for the elastic-plastic 
option, 62, 1b2/in4. 


Equivalent stringer stress, squared, when 
response is still elastic for the elastic- 


plastic option, a2, 1b2/in4 


Temporary storage used for stringer var- 
iables € or & 


Equivalent yield strain corresponding to 


SIGO for elastic-plastic option, Eo? in/in. 


Ultimate (fracture) strain for elastic- 
plastic option, es, in/in; not used 


Allowable error in displacement coefficients 
in the static solution 


Modulus of elasticity in the theta-direction 
for elastic option, Ek, psi 


Temporary value of strain, ERE > in/in. 


Strain component at the time of last yield- 


Modulus of elasticity in x-direction for 
elastic option, Ex, psi 


Temporary value of strain, en? in/in. 


Strain component at the time of last yield- 
ing, €,9, in/in. 


Temporary value of strain, eo o4 in/in. 
: ; - ier Pe 
Temporary value of stringer strain, €,,, in/in. 


Stringer strain component at the time of 
last yielding, €yy, in/in. 


Strain component at the time of last yield- 
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*FG(MB,MB) 


FM11, FM12 
FM22, FM33 


FP1(G*MBAR) 


FP2(MG*MBAR) 


FP3(MG*MBAR) 


FP4(MG,2) 


FP5 (MB*NBAR) 


FP6 (MB*NBAR) 


FP7 (MB*NBAR) 


FP8(MB,2) 


GAM(MBAR) 
GAMMA (41) 
GC 

GX (LBAR ) 


GXSTR(LBARST) 


Modal displacement coefficients for the 
initial imperfections, in. 


Stiffness constants Fyjj used in elastic, 
multilayer integration through the thick- 
ness, lb 


Displecement function, on (x), dimensionless 


m3 ¢"(x), dimensionless 
ay mm 


2 i 
~F on (x), dimensionless 
Y 
9394 
m = ; 3 
at y = 0, 1, dimensionless 
oy? 


Displacement function, on (@), dimensionless 


a one 

PAC ANS » dimensionless 
op 

a7on (8) 


» dimensionless 
38? 


a> oW(e) 
oe at 


= 0, t, dimensionless 


Integration positions in the y-direction, rad 

Not used i 
Not used } 
Zeroes of the Legendre polynomial for the 
elastic-plastic Gaussian integration through 

the thickness, &;, dimensionless 
Zeroes of the Legendre polynomial for the i 


elastic-plastic Gaussian integration through ' 
the stringer thickness, &,;, dimensionless ; 
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®GXT(NL) 
H 


HBAR 


HGO (LBAR) 


HGOSTR(LBARST) 


*HM (NL) 


*HSTR 


ICOMP 
ICOUNT 


IFIRST 


IMASTR 


INOUT 
INZ(2) 


INZSTR(2) 


IP(3*MGMB) 


Shear modulus, G9, psi 


Thickness of cross section for elastic- 
plastic option, in. 


Distance from the inner panel surface to 
the coordinate surface, H, in. 


Weighting factors for the elastic-plastic 
Gaussian integration through the thickness, 
Hj, dimensionless 


Weighting factors for the elastic-plastic 
Gaussian integration through the stringer 
thickness, Hsj, dimensionless 


Distance from the inner panel surface to 
the outer surface of layer, h, in. 


Thickness of the stringer cross-section 
for the elastic-plastic option, in. 


Not used 


Counter, initially set at 37777000000000000000B. © 


and incremented for each integration step in 
program 


Code designating first pass through subroutine 
SIGMA, dimensionless 


Code which indicates if integration point 
currently under consideration has a stringer 
associated with it 


Not used 

Code designating the appropriate layer 

number corresponding to the two panel sur- 
faces, dimensicnless 

Code designating the appropriate layer number 
corresponding to the two stringer surfaces, 
dimensionless 


Working array in RELAXP, dimensionless 


ST 


IXI (MBAR) 


JFIRST 


JL 


JLT (NPY ,NPX) 


JSTRFT 


JYJ(NBAR) 


KALT 
KB 
KC 
KDAM 
*KDS 
KERR 


KOK 
*KPG(NKP ) 


*KPB(NKP) 


KSUMA (NGNBT ) 


KSUMAS (NSTR*NBAR) 


Integer locating grid points in y-direction 
relative to pressure-mesh (Load Option 2) 


Code which indicates if the panel has 
yielded for the elastic-plastic option 


Lower index on timewise interpolation 
table (Load Option 4) 


Lower index on timewise interpolation 
table (Load Option 2) 


Code which indicates if stringer has 
yielded for the elastic-plastic option 


Integer locating grid points in 8-direction 
relative to pressure-mesh (Load Opticon 2) 


Not used 
Not used 
Not used 
Not used 
Response option code 


Dynamic response error code - 0, no error; 
1, error 


Not used 


Mesh point number (y), when paired with 
KPB, calling for printout 


Mesh point number (8), when paired with 
KPG, calling for printout 


Number of z points which have not yielded 
at each spatial station; used only in elastic- 
plastic option 


Number of stringer z points which have not 


yielded at each spatial station; used only 
in elastic-plastic option 
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*KTYPE Code designating panel type 
1, single-layer metal 
2, honeycomb panel (3 layers) 
5, multilayer panel (elastic response 
only) 


KY (LMAX ) Code in elastic-plastic response, indicat- 
ing number of times an integration point 
has yielded, unloaded, etc. 


KYSTR(LMAXST ) Code in elastic-plastic response, indicat- 
ing number of times a stringer integration 
point has yielded, unloaded, etc. 


KZ Code deciding whether the output routine 
should print 


*LBAR Number of integration points through the 
panel thickness; is assumed to be one for 
elastic option 


*LBARST Number of integration points through the 
stringer thickness; is assumed one for 
elastic option 


LMAX Maximum number of integration points being 
used; equal to MBAR*NBAR*LBAR 


LMAXST Maximum number of integration points involv- 
ing stringers which are being used; equal 
to NSTR*LBARST*NBAR 


*®LOCSTR(NSTR) Integration positions in beta direction 
matching stringer positions 


*MB Number of beta modes to be incorporated 
into the solution 


*®MBAR Number of gamma integration points to be 
used; for a symmetric boundary condition, 
only approximate half as many points are 
required 


*MBSTR Number of stringers to be used; for a symmetric. 


boundary condition in the beta direction, 
only approximate half as many are required 
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*MG 


*MGM (MG ) 
MGMB 


MGMB2 
*MOUT (MG*MB) 


MUSE(MB ,MG) 


NB 


*NBAR 


*NBN (MB) 
*NBND 
NBT 
NCALL 


NCASE 
*NCASES 
NCHPT 
*NDBUG 


Number of gamma modes to be incorporated 
into the solution 


Gamma mode numbers 


Constant, equal to the total number of 
modal combinations used, MG*MB-NNOUT 


Constant, 2*MGMB 


Gamma modes not to be included, in combina- 
tion with the corresponding NOUT mode 


Code designating which modal combinations 
are to be used 


Beta-point number corresponding to center- 
line of panel 


Number of beta integration points to be 
used; for a symmetric boundary condition, 
only approximately half as many points are 
required as otherwise 
Beta moce numbers 
Boundary condition code 
Total number of beta points monitored 
Code describing program phase: 

0, find dynamic response 

1, find static solution 

2, read data, set up constants 
Case number currently being run 
Number of cases to be run 
Not used 
Output debugging control: 

0, no debug cutput (normal option) 


1, most debug output 
2, all debug output 
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*NDERV 
NELP 


NG 
NGNB 
NGNBT 
NGT 


*®NKP 


*NL 
*NLOAD 


NLZ(2*NL) 


NMASS 
*NNOUT 


*NOUT (MG*MB) 


*NPLT 
*NPX 


Response code: l-elastic (multilayer, 
orthotropic), 2-elastic-plastic (single 
layer, isotropic) 


Response code for the elastic-plastic option: 
l-keep solution elastic, 2-allow solution to 
go elastic-plastic 


Gamma-point number corresponding to center 
line of panel 


The spatial integration point number corre- 
sponding to the center of the panel 


Total number of spatial integration points 
used 


Total number of gamma points monitored 
Number of spatial points for which printout 
of strains, stresses, reactive forces, dis- 
placements and pressure is required 
Number of layers 
Pressure load option code: 

1, nonuniform, functions 

2, nonuniform, discrete 

3, uniform, functions 

4, uniform, discrete 


Layer number corresponding to each layer's 
upper and lower surfaces 


Not used 


Number of modal combinations (<MG*MB) to 
be eliminated from the solution 


Beta modes not to be included, in combina- 
tion with the corresponding MOUT mode 


Panel type code: OQ-flat, l-curved 
Number of mesh points in y-directions for 


which pressure data is provided (Load 
Option 2) 
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*®NPY 


®NSYMB 


*NSYMG 
NTECO 
*NTIME 


NTRIAL 
NU 


NUSE(NBAR, MBAR) 


NY2 

NZP 
OTTO 
OTT1 
P(NGNBT ) 
PB(40) 
PDAM 
*PHI 


Number of mesh points in §-direction for 
which pressure data is provided (Load 
Option 2) 

Not used 

Number of stringers actually needed in 
calculations; MESTR or (MBSTRt1)/2, 
truncated 


Symmetry code in the beta-direction: 
O-symmetric, l-not symmetric 


Same as NSYMB, except in the gamma-direction 
Not used 


Number cf times for which pressure data is 
previded (Load Options 2 and 4) 


Not used 


Code indicating whether loading is spatially 
uniform or not: O-not uniform, l-uniform 


Use code for the spatial integration stations: 
O-not used, l-not used fer printout only, 
2-used for integration only, 3-both 

Constant, equal to 3*MGMB 

Not used 

Constant equal to 1/t,, 1/sec (Load Option 3) 
Constant equal to 1/tj, 1/sec (Load Option 3) 
Pressure at each spatial point, psi 

Not used 

Not used 

Angle projectile trajectcry makes with the 


normal to the panel (z-axis), $, degrees 
(Loae Option 1) 
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PI 
PIMA (MBAR) 


PINAC(NBAR) 


PPP 


®PP] 
*PPO 
* PRINT 


PRES (3*#MGMB ) 
*PRT(NTIME,NPY,NPX) 


PRTT (NPY ,NPX) 


*PT(20) 


PX (3*MGMB) 

Ql 

Q2 
RA(NBAR , MBAR) 
RFR 

RHO 


Constant, equal to Tm 


Constants associated with the equation of 
motion and Simpson's Rule; Gamma-direction, 


in2/1b-sec2 


Same as PIMA, cnly in the beta-direction, 
in2/1b-sec2 


Calculated pressure on panel (uniform dis- 
tributions only), psi 


Pressure Pj, psi (Load Option 3) 
Pressure Pg, psi (Load Option 3) 


Output frequency - integration steps per . 
printout 


Working array in RELAXP, dimensionless 


Pressure specified on panel vs time, psi 
(Load Option 2) 


Temporary storage for pressure on panel 
after interpolation on time, psi (Load 
Option 2) 


Table of uniform pressures specified on 
panel, psi (Load Option 4) 


Working array in RELAXP, dimensionless 
Constant A (Load Option 1) 

Constant B (Load Option 1) 

Range from detonation point, R, in. 
Not used 


Density of panel, 1b-sec?/in* 
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*RHOM(NL) 
PR(4) 

RRES (2*MGMB) 
RTRIAL(5) 
SAC (NL) 


SAT (NL) 


*SIGO 
SIGO02 
SIGTT1 (LMAX) 


SIGX (3*MGMB) 
SIGXT1(LMAX ) 


SIGXX1(LMAX ) 


SIGX1S (LMAXST ) 


SINB(MB*NBAR) 
SING (MG*MBAR) 


SIN2B(MB*NBAR) 
SIN2G (MG*#MBAR ) 


STRCN1 


STRCN2 


Density of layer, lb-sec?/in4 


Reaction force at corners of panel 

Working array in RELAXP, dimensionless 

Not used 

Compressive yield stresses for metal material, 
compressive ultimate stress for plastic material, 
psi; not used 

Tensile yield stress for metal materials, 

tensile ultimate stress for plastic material, 
psi; not used 

Yield stress for elastic-plastic option, psi 


Constant, equal to SIGO squared, 1b2/in" 


Stress component at time of last yielding, 
Ogg, psi 


Working array in RELAXP, dimensionless 


Stress component at time of last yielding, 
Oxg>s psi 


Stress component at time of last yielding, 
Cxx» psi 


Stringer stress component at time of last 
yielding, o,,,, psi 


Sines of beta functions, sin((n+1) 85) 
Sines of gamma functions, sin((m+l)yj) 
Sines of beta functions, sin((n)B3) 


Sines cf gamma functions, sin((m)y;) 


Constant for stringer option, bei80 
dimensionless ato 
bshs180 


Constant for stringer option, 


dimensionless athéy 


SMAX 

STT (LMAX ) 

SXT (CLMAX ) 

SXX (LMAX) 
SXXSTR(LMAXST ) 
S1A(NGNET ) 


S2A(NGNBT ) 


S3A(NGNBT ) 


S4A(NGNBT ) 


SSAC(NGNBT ) 


S6A(NGNBT ) 


S7A(NSTR*NBAR) 


S8A(NSTR*NBAR) 


+ THETAO 


*THNU (NL) 


*TITLE(20) 
TMAX 
*TNU 


Not used 

Stress component, Ugg, psi 

Stress component, Ox¢> psi 

Stress component, o,.5 psi . 
Stringer stress component, Og,,, psi 


Stress component, o™ , psi 
xX 


Stress component, Be? psi 
Stress component, an psi 


Stress component, oP >» psd 


Stress component, o2., psi 


Stress component, o-., psi 


Stri bit . 
tringer stress component, Ssxyx? PSs 


Stringer stress component, oP >» psi 


SXxXx 


Total angle subtended by cylindrical 
panel, 6,, deg, or width of flat panel, 
Dg Lic 


Poisson's ratio in the theta-direction, 
Vo> dimensionless 


Descriptive title of case 


Not used 


Poisson's ratio for elastic-plastic option 
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a = 


*TPRINE 
*TSTOP 
*TT(20) 
TTNU (LMAX) 
TTNUST (LMAXST ) 
U(NGNBT) 
| UB(NGNBT) 
| UG (NGNET) 
UU (MB ,MG) 
U1(MB MG) 


V(NGNBT ) 
VE(NGNBT ) 
VGC(NGNBT) 

VRX (2*NBAR- 2) 


VRT(2*MBAR-2) 
vs 


VV (MB,MG) 
VX0(3*MCMB) 
V1(MB,MG) 


W(NGNBT) 
WE(NGNBT ) 


Time t' (Load Option 3), sec 

Integration stop time, sec 

Time table (Load Option 4), sec 

Value of vg, dimensionless 

Value of v, for stringers, dimensionless 
Value of U, dimensionless 

Value of Ug, dimensionless 

Value of Uy » dimensionless 

Displacement coefficient, Up;,, dimensionless 


Displacement coefficients, Upn, representing 
the static conditions, dimensionless 


Value of V, dimensionless 
Value of Ve> dimensionless 
Value of Vy» dimensionless 


Normal reactive force per unit length along 
boundary, Vy, lb/in. 


Normal reactive force per unit length along 
boundary, Vg, lb/in. 


Shock velocity, equal to 5.88x107 in/sec 
(Load Option 1) 


Displacement coefficients, V,,,,, dimensionless 
Initial velocity coefficients 


Displacement coefficients, Vpn, representing 
the static conditions, dimensionless 


Value of W, dimensionless 


Value of Wg, dimensionless 
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WBB(NGNET ) 
WBBB(2*MBAR-2) 
WG (NGNET) 
WGB(NGNBT) 
WGBB(2#NBAR-2) 
WGG (NGNBT ) 
WGGB( 2*MBAR-2) 
WGGG (2*NBAR-2) 


WW(MB,MG) 
W1(MB,MG) 


XB(NBAR) 


XG (MBAR) 


XJ 


Value Weg» dimensionless 
Value Weegg> dimensionless 
Value Wy dimensionless 
Value Wy ge dimensionless 
Value Wye? dimensionless 
Value Wey» dimensionless 


2 — 
Value Wye? dimensionless 


Value of Wyy? dimensionless 


Displacement coefficients, Wnn, dimensionless 


Displacement coefficients, Wmn, representing 
the static conditions, dimensionless 


Integration positions in the beta-direction 
inches for flat panel, degrees for curved 
panel 


Integration positions in the gamma-direction, 
inches 


Constant, J. equal to 180/85 (dimensionless) 
for curved panel and t/b (inches) for flat 
panel 


Constant, 
Constant, 
Constant, 
Constant, 2J/L 


Temporary value of strain, See? in/in. 


Temporary value of strain, yi in/in. 


XKXX 
XKXXST 


XL 


*XLP 
XLP1 
XLP2 
XLP3 
XL1 
XL2 
XL3 
XL4 
XL5 
XL7 
*XP(NPX) 


XRES (3*MGMB ) 
XX (3*MBMG ) 
*®XXNU (NL) 


XX1(3*MGMB) 
X1A(NGNBT) 


X2A(NGNBT) 


X3A(NGNET) 


Temporary value of strain, ie in/in. 
Temporary value of stringer strain, 
eD , in/in. 
Sxx 
Constant, &£/m, for flat panel (inches), 
&/na for curved panel (dimensionless) 
Length of panel, &, inches 
Constant, 2L7R 
Constant, 2LR 
Constant, 1/2L2R 
Constant, 1/L 
Constant, L2 
Constant, 2/L 
Constant, 2L2 
Constant, 2L¢R 
Constant, 1/2 


X-position of pressure-mesh (Load Option 2), 
in. 


Working array in RELAXP, dimensionless 
Array composed of UU, VV, and WW, dimensionless 


Poisson's ratic in the x-direction, vx, 
dimensionless 


Working array in RELAXP, dimensionless 


Strain component, Core? in/in. 


m 


66° in/in. 


Strain component, € 


Strain component, ot in/in. 
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X4A(NGNBT) 


XSA(NGNBT ) 


X6A(NGNBT) 


X7A(NSTR*NBAR) 


X8A (NSTR*NBAR) 


*YP(NPX) 


YY (3*MGMB) 


ZA(2) 


ZASTR(2) 


ZB(2) 
ZBSTR(2) 
ZC (2*NL) 
ZCSTR(2) 
*ZEE 


ZF (LBAR ) 


ZFSTR(LGARST ) 


ZG (LBAR ) 


b 


Strain component, ea in/in. 


Strain component, e2,, in/in. 


66? 


Strain component, eb in/in. 


xe” 


Ti 


Stringer strain component, foxx? in/in. 


Stringer strain component, eP Si in/in. 


Y-position of pressure-mesh (Load Option 2), 
in. or deg 


Array composed of acceleration coefficients, 


Uns Vans tems Lfeee* 


ZB normalized with a, ZB/a, inches for flat 
plate, dimensionless for a curved panel 


ZBSTR normalized with a, ZBSTR/a, inches for 
flat plate, dimensionless for a curved panel 


I+ 


layer thickness/2, in. 


+ 


HSTR/2, in. 
+ layer thickness, in. 


+ HSTR/2, in. 


Distance from panel to detonation, Z, in. 
(Load Option 1) 


ZH normalized with a, ZH/a, inches fora 
flat plate, dimensionless for a curved panel 


ZHSTR normalized with a, ZHSTR/a, inches 
for a flat plate, dimensionless for a curved 
panel 


‘ 2 2 ‘ . 
Gaussian station squared, G39 dimensionless 
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ZGSTR(LBARST ) 


ZH(LBAR) 


ZHSTR(LBARST) 


2Z1(9) 


Gaussian stringer station squared, a 
dimensionless 


z coordinates corresponding to the Gaussian 
integration points through the panel thick- 
ness for the elastic-plastic option, in.: 


z coordinates corresponding to the Gaussian 
integration points through the stringer 
thickness for the elastic-plastic option, in. 


Not used 


B.2.4 Program Operation 


The DEPROSP program is written in FORTRAN IV and consists 


of 19 user supplied routines on approximately 3300 cards. The 


code was developed on the Control Data Corporation (CDC) 6600 


computer under the NOSBE1 operating system. 


In order to minimize the amount of central memory core re- 


quired to execute the program, the user should eliminate at 


least one of two options prior to loading. 


is either the elastic static solution capability, where subrou- 


tines RELAXP, LIST1 and SOLVE are required, or the elastic- 


plastic option where subroutines LIST2 and SIGMA are required. 


For stringer calculations, the subroutine SOLVE is also needed 


in the elastic-plastic option. These options are outlined in 


Tables XIIand XIII, and the corresponding program core requirements 


are given. Elimination of the unnecessary subroutines can be 


accomplished by removing them from the deck (or file) completely, 
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This choice of options 
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replacing them with dummy routines, or using an SLOAD instruc- 
tion to selectively load the required routines. The use of 
blank common enables the program to load and execute at the 
same core level. Input and output are equated with logical 
files TAPES and TAPE6, respectively, and there are no addi- 
tional tape or disk file requirements for operating the code. 
The FIN compiler has been used to compile the code under 
"OPT=1" and "R=2" options. Compilation requires approximately 
80,000 words and 24 seconds CP time. 
Computation time will vary considerably, depending on 
the panel, the complexity of the model, whether or not the 
solution goes inelastic and the response time requested. 
A very rough approximation for an inelastic response is 6 x 107" 
CP seconds per mode, per integration grid point, per time step 


of response, for a non-stringer case. 


TABLE XII. CORE REQUIREMENTS FOR MAJOR PROGRAM OPTIONS 

ore Required) 
Program Response | Input Parameters| Subroutines |to Load and 
Option NDER Eliminated _|Execute* 


2 1 RELAXP, SOLVE] 151 Kg 
LIST2, SIGMA 
LIST2, SIGMA | 226 Kg 
RELAXP, SOLVE,} 254% Kg 
LIST1 
oe ad etn is} 
3 2 LIST1 330 Kg 


*Code compiled under FIN option 1 


Elastic-Plastic 
(Dynamic) 


Elastic-Plastic 
with Stringers 
(Dynamic) 


Static, followed 
by Elastic- 
Plastic (Dynamic) 


TABLE XIII. PANEL RESPONSE OPTIONS 


Panel Type 
(KTYPE) 


Response Options 
(NDERV ) 


*The program forms an equivalent single 
layer for KTYPE = 3, NDERV = 2 
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B.3 Program Input Data 

Specific instructions are provided in this section to 
enable the user to provide the necessary card input for the 
DEPROSP program. These instructions are identical to those 
of the DEPROP program (18) except for the cases where stiffeners 
are included. Hopefully the program could be run without the 
use of Reference (18) but for greater detail the user is re- 
ferred to that reference. The stiffeners are included only 
for the single layer and the three related response options 
of static only, dynamic only, and static followed by dynamic. 
The stiffener is assumed to have the same mechanical properties, 
the same boundary conditions, and same stress-strain relations 
as the single layer for either elastic or elastic-plastic 
response. 

The input data are specified in groups, where each group 
begins on a separate card. More than one card may be required 
for a group, however. The variable type and format correspond- 
ing to each data group is given in parentheses in the input 
instruction and is always in fields of 12. For convenience, 
floating point numbers can be left justified in the field as long 
as the exponent is right justified. Also, zero values can be 
replaced by a blank field. Columns 73 through 80 are not used 
for data and can be used for card identification or other 


purposes. 
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All input parameters, where appropriate, should be compared 
with the maximum dimensions provided for in the program, as 
delineated in Table X. This is very important since the program 
does not attempt to check the input for such violations. 

| Group 1 provides for the execution of several jobs (or 
cases) in the same run. All subsequent data groups (2 to 23) 
must be repeated for each case. 

The panel type, response option, and debug options are 
specified in Group 3. It is important to check Table XI to 
insure the correct sub routines are included to the response 
option selected. It is suggested that the first debug option 
(NDBUG=0) be selected. 

Group 4 contains the number of modes to be used in the 
solution and the number of integration points to be used. The 
accuracy of the solution is based on the degree of convergence 
of stress and strain quantities. These quantities converge 
less rapidly than the radial displacement. Also, cases involv- 
ing a clamped edge condition will converge less rapidly than 
simply-supported cases. Since both computer time and accuracy 


increase with more modes and points, a trade-off usually be- 


comes necessary. It has been found that about 25 modes give 
an acceptable general solution for most panels. But more modes 


are required for better accuracy in determining edge strain and 


| 
| 
| 
| 
| 
| 
| 


reaction force quantities for clamped panels. 
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The actual mode numbers are specified in Groups 5 and 6. 
The maximum value that the mode numbers can assume in the pro- 
gram is 19. When symmetry is taken in either direction (Group 
7), or if the pressure loading is symmetrically oriented, only 
the odd numbered modes (1, 3, 5, ...) are required in that 


direction. 


Spatially, the desired number of integration points (MBAR 
and NBAR) for a full panel should be approximately two times 
the maximum mode number used in that direction, plus three. 
However, when NBN or MGM is large, this condition may not be 
satisfied for nonsymmetrical panels, since MBAR and NBAR are 
dimensioned at 23 in the program (see Table X). For symmetric 
solutions, MBAR (or NBAR) need only be approximately one-half 
the value for a full panel since only one-half (or one-quarter) 


of the panel is actually analyzed in the solution. For a non- 


symmetric condition, MBAR (or NBAR) must be an odd number. For 
an elastic-plastic solution, a minimum of four integration points 
through the thickness is recommended, and a maximum of six is 
provided in the program. | 
Information relating to the use of stringers is specified 
in Groups 8 through 10. In order to be included in the calcula- 
tions, the stringer locations must coincide with beta integra- 
tion points selected in Group 13. It is suggested that the number 
of integration points through the thickness be at least that for 


the panel, with a maximum of ten allowed in the program. 
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In Group 11, the user is given the option of a purely 
elastic solution, or an elastic-plastic solution. The elastic- 
plastic option will tend to be slower and require more computer 
memory. It should be noted that honeycomb panels are reduced 


to an equivalent single-layered panel for elastic-plastic 


response. Again, Tables XII and XIII should be consulted to insure 


the program is compatible with the response option selected. 

Groups 12 and 13 provide a mechanism for selecting a maximum 
of 49 modal combinations from a 13 by 13 combination array 
(MG=MB=13). Thus, the more significant modal combinations for 
an optimal solution with respect to accuracy and computer time 
can be selected and the other combinations eliminated. A 
general rule of thumb is to eliminate the higher frequency modes 
which are usually associated with modal combinations having the 
larger MG+MB values. An example of this would be the selection 
of MG=MB=7 for a symmetric problem, but eliminating 24 combina- 
tions as indicated inFigure 27. The relative importance of each 
modal combination can be evaluated by examining the response 
output and comparing the magnitude of the displacement 
coefficients. 

Groups 14 and 15 are responsible for selecting the points 
in the integration grid for which printout of strains, stresses, 


displacements, reactive forces (or boundaries), and pressures is 
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Figure 27. Example of Modal Selectiv.. 


required. Strains and stresses are computed at the inner and 
outer surfaces of the panel Each point in the grid is designated 
by a pair of integers, the first integer referring to the gamma- 
position, the second to the beta-position. Actual positions are 


found from 


* & (I-1) I = Lgjsceght 
ala: (M-1) (symmetric in x-direction) 
- (I-1) I = ‘oe | 
M-1 (full in x-direction) 


and similar expressions for y (or @). For example, the corner 

point in a symmetric panel would be numbered (1,1); the center, | 
(MBAR,NBAR). It is suggested that these points include the 

stringer locations so that the effect of the stringer on the 


strains and stresses can be displayed. 
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Group 23 contains the modal components, 6,,, for the 


initial radial imperfections. The analyst must compute the 


i 
: 
i 
3 
i 


6mn's from measured data using the integration technique applied 
to Fourier series coefficients. Generally, such data will not 
be available, and zero values should be specified for the §pn's. 
The capability of considering initial imperfections also enables 
the analyst to determine the sensitivity of panel response to 
initial imperfections. 

Group 24 provides the integration time increment, the 
response stop time, and the printout interval. If the user 
specifies a zero time increment, the program computes an appro- 
priate At which, in most cases, will give a stable sclution. Be- 
cause it is approximate, the analyst may want to make comparable 

r runs using different At's. In general, an elastic solution which 


is numerically stable will be accurate. Hence, the optimum At 


is the largest which remains stable. For an elastic-plastic 
solution, however, the accuracy of the solution may deteriorate 
Slightly as the point at which the solution diverges is approached. 
For an elastic-plastic solution with stringers, a smaller At may 
be required. Once a time increment is selected, it should be 
valid for moderate changes in response level. 

Although the stop time can vary a great deal, the total 


number of integration steps required to capture peak response 
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will be roughly between 500 to 1500. One exception to this 
may be a curved panel experiencing snap-through buckling, in 
which case considerably larger response times may be required. 
A printout frequency of once every 20 steps is usually adequate 
for monitoring the response time history. 
Groups 25 to 36 provide for the appropriate pressure load 
on the panel. The user should refer to Figurel6 orl? for a 
definition of certain input parameters. 
Group 1: (112) NCASES 
Number of cases to be run (NCASES) 
Group 2: (20A4) TITLE 
Identifying title per run, date, etc. Free 
field. . (TITLE) 
Group 3: (3112) KTYPE, KDS, NDBUG 
Code designating panel type (KTYPE) 
1, Single-layer panel 


3, Honeycomb panel (3 layers) 

5, Multilayer panel (elastic response only) 
Response option code (KDS) 

1, Static only 

2, Dynamic only 

3, Static, followed by dynamic 
Debug option (NDBUG) a. | 


0, No debug output 


152 


et deren mme ae eS 


See aca eae 


1, Most cebug output 
2, All debug output 
Group 4: (5112) MG, M83, MBAR, NBAR, LBAR 
Number of gamma modes to be used (MG) 
Number of beta modes to be used (MB) 
Number of gamma integration points actually 
used over the portion of the panel analyzed. 
Must be an odd number for full panel (see 
Group 7) (MBAR) 
Number of beta integration points actually used 
over the portion of the panel analyzed. Must 
be an odd number for full panel (see Group 7) 
(NBAR) 
Number of z integration points used through the 
thickness. (LBAR) [Not needed for NDERV=1 
(see Group 11)] 
Group 5: (6112) (MGM(I), I=1, MG) 
Gamma mode numbers, m 
Group 6: (6112) (NBN(I), I=1, MB) 
Beta mode numbers, n 


Group 7: (2112) NSYMG, NSYMB 


Symmetry code in gamma direction (NSYMG): 


1/2) | 


0, Symmetry assumed (O<y< 
1, No symmetry (0 < y < 7m) 
153 
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Symmetry code in beta direction (NSYMB): 
B < sf2) 


[A 


0, Symmetry assumed (0 


JA 


1, No symmetry (0 < B < m) 
Group 8: (2112) MBSTk, LBARST 
Number of stringers to be considered (MBSTR) 
Number of z integration points used through the . 
stringer thickness (LBARST) [Not needed for 
NDERV=1 (see Group 11)] 
If MBSTR=0, skip to Group 11. 
Group 9: (2F12.1) BSTR, HSTR 
Width of stringers, in. (BSTR) 
Thickness of stringers, in. (HSTR) 
Group 10: (6112) (LOCSTR(I), I=1, MBSTR), if NYSMB=1 
or (LOCSTR(I), I=1 , (MBSTR+1)/2), if NYSMB=0 
Position of stringers. Must coincide with beta 
integration points actually used (LOCSTR) 
Group 11: (3112) NPLT, NBND, NDERV 
Panel type (NPLT): 
0, Flat panel 
1, Cylindrical panel 
Boundary condition code (NBND) 
y-direction 8-direction 
1, Clamped-clamped; Clamped-clamped 


2, Simple-simple; Simple-simple 
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3, Clamped-clamped; 
4, Simple-simple; 
5, Clamped-simple; 
6, Clamped-clamped; 
7, Clampedc-simple; 
8, Simple-simple; 


9, Clamped-simple; 


Note: Whenever a clamped-simple condition is selected, the full 


panel is analyzed in that direction, and NSYMG, NSYMB, 


Simple-simple 
Clamped-clamped 
Clamped-clamped 
Clamped-simple 
Simple-simple 
Clamped-simple 


Clamped-simple 


MBAR and NBAR should reflect this. 


Response option (NDERV): 
1, Elastic only 
2, Elastic-plastic 


Group 12: (112) NNOUT 


Number of modal combinations to be eliminated 


from solution (NNOUT). 
(G < NNOUT < MG*MB) 
If NNOUT=0, skip to Group 14. 
Group 13: (2112) MOUT(I), NOUT(I) 
Gamma mode. (MOUT(I)) 
Beta mode. (NOUT(I)) 


Repeat Group 13 for I=1, NNOUT. The cards in Group 13 may be 


arranged in any order. 
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Group 14%: (112) NKP 

Number of spatial points at which printout of 

stresses, strains, displacements, reactive forces 

and pressures are requested. If NKP=0, all of 
the above information will be suppressed. 

(NKP) 

If NKP=0, skip to Group 16. 

Group 15: (2112) KPG(I), KPB(I) 

Integration point in gamma-direction at which 
printout is requested. Points are ordered 
1-MBAR, beginning at y=0, and evenly spaced 
from there. (KPG(I)) 

Integration point in beta-direction at which 
printout is requested. Points are ordered 
1-NBAR, beginning at 6=0, and evenly spaced 
from there. (KPB(I)) 

Note: These two indices are taken as pairs where each pair 
designates a particular spatial point. The pairs may 
be specified in any order. 

Repeat Group 15 for I=1, NKP. 

Group 16: (112) NL 

Number of layers. (NL) 
(NL must be 1 for KTYPE=1, and 3 for KTYPE=3) 
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Group 17: (3F12.1) XLP, THETAO, A 
Full length of panel, 2, in. (XLP) 
Full width of flat panel, b (short 
direction), in. (NPLT=0) 
or (THETAO) 
Full subtended angle of cylindrical 
panel, 85> deg. (NPLT=1) 
Radius of cylindrical panel, in. (A) 
(Not needed for NPLT=0) 
If NDERV=2, skip to Group 21. 
Group 18: (2F12.1) HM(I), RHOM(I) 
Distance (h) from inner panel surface to the 
outer surface of layer I, in. (HM(I)) 
Mass density of layer I, lb-sec?/in*. (RHOM(I)) 
Group 19: (5F12.1) EX(I), ET(I), XXNU(I), THNU(I), GXT(I) 
Modulus of elasticity in the x-direction, psi. 
(EX(I)) 
Modulus of elasticity in the theta-direction, 
psi. CETC(I)) 
Poisson's ratio in the x-direction. (XXNU(I)) 
Poisson's ratio in the theta-direction. (THNU{(I)) 


Shear modulus, psi. (GXT(I)) 
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(2F12.1) SAT(I), SAC(I) 
Tensile yield stress for metal panels; tensile 
ultimate stress for plastic panels, psi. 


(SAT(TI)) 


Absolute value of compressive yield stress for 


metal panels; absolute value of compressive 
ultimate stress for plastic panels, psi. 


(SAC(I)) 


Repeat Groups 18-20 for I=l, NL. 
Skip to Group 23. 


Group 21: 


(3F12.1) HM(I), RHOM(I), EM(TI) 

Distance (h) from inner shell surface in the 
outer surface of layer I, in. (HM(I)) 

Mass density of layer I, lb/sec?/in*. (RHOM(I)) 


Modulus of elasticity, psi. (EM(I)) 


Repeat Group 21 for I=l, NL. 


Group 22: 


Group 23: 


(4F12.1) TNU, SIGO, EP, EPSIF 

Poisson's ratio. (TNU) 

Yield stress for a metal panel, psi. (SIGO) 
Strain hardening modulus (Et), psi. (EP) 
Ultimate strain, in/in. (EPSIF) (not necessary) 
(6F12.1) ((FG(N,M), N=1,MB), M=1,MG) 

Modal displacement coefficients for initial 


radial imperfections, in. (FG (N,M)) 


Group 24: (3F12.1) DELTIM, TSTOP, PRINT 
Integration time increment, sec. If DELTIM=0.0, 
the program determines the time increment re- 
quired for stability. (DELTIM) 
Integration stop time, sec. (TSTOP) 
Print frequency (integration steps per printout). 
If PRINT=0.0, printout of intermediate data 
will be suppressed. (PRINT) 
If KDS=2, skip Group 25. 
Group 25: (F12.1) PS 
Uniform static pressure load, psi. Can be either 
positive or negative value. 
' If KDS=1, skip Groups 26-36. 
Group 26: (112) NLOAD 
Dynamic load option 
1, Special Eglin analytical function over 
space and time. (See Figure 16) 


2, Discrete point by point, time by time dis- 


tribution 
| 3, Spatially uniform, with an analytical 
function for time histery. (See Figure 17) 
4, Spatially uniform, with a discrete time 


history. 
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If NLOAD=2, skip to Group 28. 
If NLOAD=3, skip to Group 34. 


If NLOAD=4, skip to Group 35. 
Group 27: (2F12.1) ZEE, PHI 


Distance of detonation from panel, Z, in. (ZEE) 
Angle projectile trajectory makes with the normal 


to the panel (z-exis), $¢, degrees. (PHI) 


- 


Skip Groups 28-36. 
Group 28: (3112) NPX, NPY, NTIME 


Number of spatial points in the gamma-direction 
at which pressures are to be specified. (NPX) 
(Must be at least 2) 

Number of spatial points in the beta-direction at 
which pressures are to be specified. (NPY) 
(Must be at least 2) 

Number of times specified in the pressure-time 
history. (NTIME) (2 < NTIME < 6) 

Group 29: (F12.1) DTIM 

Time interval between samplings (DTIM). The time 
history for each point has the same time inter- 
val, but distinct delay times (Group 32) 

Note: Be sure to allow for first point to be engulfed at time=0. 


Program will extropolate data past last time in table. 
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Group 30: (6F12.1) (XP(I), Isl, NPX) 
x-positions at which time histories are specified, 
in. (XP? 
CGFIZ.1) CYPCI), T=1, NPY) 


y-positions at which time histories are specified, 


in. (YP) 
Group 32: (6F12.1) (DET(J,1I), J=1, NPY) 

Delay time for pressure wave to reach grid point, 
sec (DET) (One point must have delay time of 
zero) 

Repeat Group 32 for I=1, NPxX. 
Group 33: (6F12.1) (PRT(K,J,I), K=1, NTIME) 


Pressure for each time and grid point, psi (PRT). 


Repeat Group 33 for J=1, NPY. 
Repeat Group 33 again, for I=1, NPX. 


Skip Groups 34-36. 
Group 34: (6112) PP1, PPO, TTO, TPRIME, AA, ANN 


Pressure, Pj, psi (FP1) 
Pressure, Pg, psi (PPO) 

Time, tg, sec (TTO) 

Time, t', sec (TPRIME) 
Parameter a, dimensionless (AA) 


Parameter n, dimensionless (ANN) 


Note: See Figure 17 for definitions. 


Skip Groups 35 and 36. 


A A LLL I DRA RETR BRITE ANE mt 


Group 35: (112) NTIME 


Note: 


Group 36: (2F12.1) TT(I), PT(I) 


Note: 


Repeat Group 36 for I=1l, NTIME. 
Repeat Groups 2 to 36 for each additional case, as specified 
in Group 1. 


B.4 


The output for DEPROSPis directed to the printer. Although 
program output is largely self-explanatory, the normal output is 


described in detail in Table XIV. Where possible, the correspond- 


Program Output 


Number of points to be specified in point-by- 
point load description. (NTIME) 
(2 < NTIME < 20) 
Be sure to include time=0 and also an end time which 
exceeds TSTOP. Otherwise, the last value in the table 


will be used. 


Time, sec. (TI(I)) 
Pressure, psi. (PT(I)) 


One time and one pressure per card. 


ing program variable is given parenthetically. 


Certain errors, if detected by the program during execution, : 


are brought to the user's attention by means of a printed error 


message. Table XVprovides a list of such messages, along with an 
indication of the routine associated with the message and the sub- 


sequent action the program takes. In most cases, the program will 


cycle back to attempt the next case if it cannot continue with 


the current one. 
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TABLE XIV. DEPROSP STRUCTURAL RESPONSE QUTPUT 


Time-History Output 


Time from shock arrival, sec (TIME) 


AR IP Ni ihe OER IRS 2 init Nok Rp NN 


Normalized axial, tangential, and radial displacement modal 
coefficients for all modes, with the beta mode index varying i 
most rapidly ((UU(J,1), VV(J,I), WW(J,1I), J=1,MB),I=1,MG) ' 


Table of stress-strain information for inner and outer surfaces 
at each grid point selected: 


Flag ("S") indicating stringer portion of stress and strain 
X coordinate, in. (XG) 

Beta position, in. or deg (XB) 

Depthwise position, in. 

Axial strain, dimensionless 

Circumferential strain, dimensionless 

Shearing strain, dimensionless 

Axial stress, lb/in 

Circumferential stress, 1b/in?2 

Shearing stress, 1b/in? 

Flag ("*") indicating equivalent strain has exceeded yield 
strain (elastic runs only) 

Counter indicating number of unloading and reyielding (KY) 
(elastic-plastic runs only) 


Table of reactive force information for each grid point selected: 


Normal reactive force (Vy or Vg), lb/in. 
(VRX or VRT) 


Tangential reactive force (Ny, or Ng), 1b/in. 
(ENX or ENT) 


Reactive forces at corners: 


Reactive force (R), lb . 
(omitted for panels clamped on all edges since forces are 


all zero). 


Table of displacement-pressure information at each grid point 
selected: 


cit eee aaa 


X-coordinate, in (XG) 
Beta-position, in or deg (XB) i he eects 


—— == ses = Gs 3 = > ain 
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Saree, 


TABLE XIV. (Concluded) 


Tangential displacement, in (VF) 
Radial displacement, in (WF) 
Pressure, psi (PPP) 


Summary Output 
Message indicating whether run was terminated normally or 


abnormally, and the time at which computations stopped, sec 
(TIME) 


Net CP time for response, sec (CPT) 


Number of integration points which yielded, if any 
(Elastic-plastic response only) 
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TABLE XV. ERROR MESSAGES 


CANNOT TOTALLY CORRECT FOR OVERSHOOT. XXX (PANEL or STRINGER). 


An iterative process to correct for overshoot associated with 
yielding has not convergec in five trials. This probably means 
a numerical instability is creeping into the solution. Program 
Continues until such errors occur 100 times. (SIGMA) 


DEPROSP IS ABORTED AT T, SEC = XXX. 


DEPROSP cycles back to attempt next case. This case is aborted. 
(DEPROSP) 


EPP IS OUT OF RANGE (PANEL or STRINGER). 


Numerical instability detected. A smaller At may be required. 
This case is aborted. (SIGMA) 


IMMEDIATE RELOADING, XXX (PANEL or STRINGER). 


Probable numerical instability creeping into solution. Program 
continues until such errors occur 100 times. (SIGMA) 


SINGULAR MATRIX IN S/R SOLVE. 


The relaxation process has generated a singular matrix in deter- 
mining static equilibrium or a singular matrix has been generated 
during a stringer solution of the equations of the 2nd order 
differential equations of the axial or radial displacement com- 
ponents as in Eq. (33). Program aborts this case. (SOLVE) 


SOLUTION DIVERGING IN DEPROSP 
Very large accelerations have been computed in DERV2, indicating 
a numerical instability. A smaller At may be required. This 
case is aborted. (DERV2) 

SOLUTION DIVERGING IN RELAXP. 


The iterative process to find static equilibrium has failed. 
Program aborts this case. (RELAXP) 


SOLUTION IS UNSTABLE. 


Numerical instability has been detected in elastic-plastic solu- 
tion. A smaller At is required. This case is aborted. (SIGMA) 


TABLE XV. (Concluded) 
STRINGER CALCULATION. 
This message precedes another error message signaling that the 
problem was encountered while calculating a stringer stress or 
strain. (SIGMA) 
THE VALUE OF LBAR IS INVALID. LBAR = XXX. 


An incorrect value of LBAR has been specified. This case is 
aborted. (LEGEND) 


TEE VALUE OF LBARST IS INVALID. LBARST = XXX. 


An incorrect value of LBARST has been specified. This case 
is aborted. (LEGEND) 


TOO MANY TRIALS IN STATIC SOLUTION. MTR = XXX. 


To avoid looping indefinitely in attempting a solution repre- 
senting static equilibrium, an upper limit of 10 is placed on 
the number of trials. Program may need more trials and adjust- 
ment of the variable CON in RELAXP. Program cycles to next 
case. (DEPROSP) 


**WARNING** INCONSISTENCY IN SYMMETRY 
A clamped-simple boundary condition has been specified, while 
a symmetric solution has been indicated. Program continues. 
(DSET3) 

**WARNING** - TIME EXCEEDS TABLE 
Either TSTOP should be reduced or the time-pressure table 


extended. Program uses last value in table and continues. 
(Load Option 4) 


ore AIRE EEN IRRE cM mR = 


B.S 


ayer 
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any 


GMD 


59 


ary 


7.9 


anny 
1199 


#29 


4%04 
1499 


Program Listing (DEPROSP) 


PROG24™M MEPPASE (INPUT,OCUTPUT, TAPSS=TNPUT,TAPFS=OUTPUT) 

SOMMON ZFIRST/S ITCOUNT \ 
COMMON/SSNOWYAS = CEIT(5) DELTIM, GAMMA(41),TCIMP, INOUT, KALT gK8y 

1 KIAMA KIS 9 CE 27g KOK K TYPE SNOALL gNTESE SNCHOTZNNRUGNYASS sNIRIAL» 
2 PBI47) ePDAY, OPP ,PRINT,ZRE 2, 2TPIAL(5Y, TIME, TITLE (C20), TSTOP, 

3 771449) 

SOMMOM SPLOAI F74,P50,TTICy»TIPF IME AN ANNZOTT1, OTT IgAZy 

1 JLeNTIMF eNLOSD,OT (C290 TT (20) g FEF SPH 415279 VS,y 

2 NET CANSADIeENAKeMNPVGDTIMy FT (A161 7DXPCIOD  VOCIND, 

BP TXT CATV, IVICA TV, ILI CIM 199 ,P OTT (1% ,19),0%1 023) ,0Y1 (23) 


FOOMAT(451T12) 
FOPMAT(5F 92019 
FORMAT (2NA4) 
NCAS= 
Inout 
tcom? 
Icouvt = 3777799000009 70F0003 
REAL (5,1) NEOASTS 


oot 
is 


PEAD (5,3) CFTTLE CT) gT=1520) 
NCAS= = NCAS= + 1 

KEPR = 9 

NTPTAL = 0 

KNAM = ? 


READY (5,1) KTY27- KOS sNARUG 
IF CINOJTLEN.M) SO TO 1400 
WRIT=(5, 7099) (TITLE (CT), T=1,20) 
GO TD C3BNDGOV,5SINAANN A 7COI, KTYOF 
ARIT= (46,3590) 

sO T) £950 

WOTT= (A, F409) 

60 TY? ihS 

WETT (CF, 27009 

69 79 195"/ 

WRITE C4, 3899) 

£0 T) 19690 

WETTE (5, 3999) 

60 Td 1969 

$0 TD (L109 91290 13090, KAS 
WRIT= 05,4390) 

GO. ¥3 1466 

AQPTT~ 0594499) 

GO T? 1400 

W2TT> (45,4590) 

NCALL = 2 

GREE EV 

TF (45 99,6T.9) GSO TN 1609 
NCOLuL = 4 

CALL PINITCO) 

Pact io he Be 

TF €€956EN1) 5D TA 14699 
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AAAI TELE Ea TE AA STB ABI A pt 


1570 


1499 
1705 


zanA 
3590 
36H2 
3700 
3805 
39090 
&3Nn9 
G&D 
450) 


IF (4E>2.6769) ©9 TD 1500 


NTALL 
KOK = 


CALL PINTT(1) 

ATRIAL M1 )=1.% 

= NFRIAL +. f 

CALL PRIP 

TF (NSASELLTLNTASES) GO Te 100 


NTRIAL 


STOP 


TAP MAT 
FO2MaT 
FOOMAT 
FOOMAT 
FORMAT 
FOSMAT 
FOPMAT 
FORMAT 
FORMAT 
ENO 


0 


(1H1,39¥,13H9 © BPO SF PL/LY¥,20A4) 
TABHASINGLESLAVER MITTAL PAMEL ’ 
(3NHOSINGLFeLAYFR PLASTIC PANEL ) 
{2S5HOHCNTYCOMA METAL PANEL ) 
C27HONAONEYCO“S PLASTIC PANEL ’ 
(29HOMULTI-LAYFR PLASTIC PANFL ’ 
(214HOSTATIC SOLUTION ONLY) 

(22HONVNAMIC RESPONSE ONLY) 

C37HOSTATIT SOLUTION AND OVNAMIT RESPONSE) 
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Se se Cd 6d 


qd oo 


oe ae) 


19) 


1 
2 


1 
2 


SUBROUTINE ROLT 


THES SUSROUTENE ‘SETS (UP iH 
SES CRis 


TOMMINATALKAS Neg TYMASTIVK7 SL IO AL PAPSTLYAXGLMAXST,. IOST2(6) MB, 
MIN PMASTI, Gy MGMT TD MIMD, MGMAMMUSE (13 e150 go NIV IAPZNBN(13), 


MINDNIT, INT IVINGSNONG, 
SOMMINSTALKAS FETICPRIACE 


C< 04) 9 COS86299),COS6029 9)» S982 70299) 50525 (299) 4 TARTS WIAT1, 


F324 (799) ,FO 71299) ,FPR(9 

F27(72799) pFPRC13,2) 6 

GNM (739 4KS ePIMAC23) PIN 

© IN22 (299) ,STN2G (299) ,¥ 

MLO SeYL te XL A eXL Fe XL Gy XL 
DIMENSION €010°99,092(270), 
DIMENSION COL@20), COAC20) 


DATA SOV71 oe SISEAUS7TI1L4 yg 26 GAIA7TS2GH72 15, TaSOIDILG7IGSg 4. 49999953 847, 
5 F090 999199 G5 5049999999 I15 sg 7e 5 Fed Fe ye SLL e517 e513 Sy 


14059150 9915 0 54170 5918S 

SATA 2f270.98259221656%,4 1. 
1670099146989 ,9.99999993 
LSP Pai0:7 


DATB CVI 1 6269275 33595 42 oe 2G4IIDITA 925 4 36 249999999595 402595025 ye 25y 


725g he 25970 25917025 GLI 
170255190 2591902592 .25/ 
DATA 394/71099977731192,1.0 


FAC = S9°T(?.0) 

90 190 T=194 

wee: = RAG 

It = ¢ 

GO TD (509,799,509,700,900 


CLAMIEN = CLAMBEN, GAMMA, 


DO S20 Jetty 

M = “4GMOT) 

SOLITh = COLUY} 

GOAT) = CHh2t) 

09 FAN Mets 

Xft = COL CM) 

¥2 = COACM) 

90 499 Y=i,NcT 

Ries hoe IF 

¥3 = YLFGAM(T) 

Fy™( x?) 

FyYD(=¥3) 
= SINCY) 

CL = COS (YF) 


" 

4 

_ 
“om 


MONE SHAPES FO? JOUNDARY TONNITTONS 


NGN3T, NGYT, NOLT NSTO NS YM3, NSYMG,PL 
1110 3)9,CC201 3) ,FC5013),CCO(1L 304 


FIV, FOG (13, 2945512999 ,F95 0299), 
&(239_9SINF(299) »SINGI299), 
SsX IS 2 eV S38 eK SG aX IF oe VL KLM PL ye XLO2, 


SeXL7eS TRO“ 1,STRON2 
797029) ©9094 020) 


21Ge 5.2L -5/7 
NON7 77311%9496 9999664591246 
733551.2.99900009 270,%.99993999998 B81, 


25 gL 2 PMG at Fey Lhe 259152591 G25 yg 


1000144989,1.00 9909 70269,1L7¥*1.07 


9500 49079790,990) 4 NONI 
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SE a TE ATS TIE Ma i Ala nial Nii 


onl 


FPECHE) sft # X2*SL + eS*EL. X20 EX 2 SEL eX 2) SEXY 
FP2¢TT) KEFCSL © BAFCL = oSP CLs t¥2ZVFEK2? + 65 *E16-K29 FEXT) 
FPSeIT) NUPSP RCCL — KASSEL + OP CLOEX2VPEX? + 2 5¥( Le X2) VEX) 


EF CESED LE FPGE MI) = —2.8*X 2X TFTS 


EF €PeEYoMIAM FOR 42) = KXLFFS*O-SL = K2FCE + .S*C0L,. — K2* 


1 Et] fe & KNFPSELSES 
SONTINUE 

CKUSY = Ye.7FAC 

GO 3 10990 


SIMOLY = SIMOLY, GAMMA. 


AGN M=1—™5 
= MGOW(M) 
any TH=ieNnst 
FE 4 & 
X1*GAM(T) 
x3 SIV¢%2) 
FPLt{rt) Xz 
FP2(TT) XiF0os(¥2) 
FP3(TT) “XLF¥2*xX3 
EF €£oGie 2 <PGtee Pils =KPs*s 
IF (L2EVeMBAP) FOU(M,2) = $X¥LF*Z*¥COS(X2) 
CONTINUE 
CK(S) = FAC 
59 ™7 1000 


CLAMPED - SIMPLY, GAMMA. 
D0 929 T=14M5 

M = 4GM(T) 

SOLET) = COZ) 

CDACT) = COG) 

$0 12 54 


Ebe y 
GO TI €1109913974170741106 5110 991500 6135N05150901509) » 


CLAMPSN = GLAMPEN, BETA, 


O00 1120 T=i1,4% 
NM = N8BNCT) 
SOL(ty = COIN) 
CI4¢I1) = CO2KN) 
90 1209 N=1,M8 
X21 = COLON) 

X¥2 = COA(M) 

99 1790 J=1,N9°T 
as S31 + 2 

MZ = VIF BETO S) 
€wi = Fy>(x3) 


NBNO 


i , 1 
1200 


carved 


ESCH 


1400 


Gace 


150 


1529 


“4600 


Fx? = FYDO(-Y%) 

StL = SINCX3) 

cL = COS(X3} 

FPSCLITY Ss -€l + XOFStl © SFL, 4X2) 4EX Oe) ASP (4, XK 2S FS KS. 
FPRCIID = KI*¥CSL + K2*CL = .5¥ (12X27) FFX? + 254 (10-02) FEX1) 
FR7(TI) = X9**¥2¥*(CL = X2"CL + S¥ Let XA FEV? & SIL. - X29 FEX1) 


TF (JeE201) SORT Nol? 
TF (J.E2.NSA7) FPR( Ny?) 
“EX (1. + X2)¥*EX?D) 

SONTINUE 
cK (48) 1./FAC 
G9 T2 44507 


~2 DEK OFX HHS 
KLFFZE(ASL 


= KOFCL! # 5 oF C04 s 


x2)* 


MP = SIMPLY, PETA, 
1490 N=i,¥% 
NANON) 

1400 J=1,NAT 

Pr 4 
YIFRFTRC YS) 
SIN(Y¥2) 

FPS{(ITI) x3 

FP6tITd KXLFF OSX?) 
FP7(IT) =KX1**2*Xs 
IF (JeEQe1) =PRIN,1) 
IF (JeFDeNSA2d FPBIN,2) 
CONTTNUE 

CK(45) FAC 

GO 1) i600 


-X1F¥3 
“XL *¥*3* COS (1X2) 


CLAPPED =" SISPEEs 
N9 1520 T=i, 5 

N NANCT) 

COLT) CIOZINY 
COAT) COGEN) 
GO T7 1149 


SETA. 


PETUPN 
END 
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a? 


1¢n 


SURSIUTINE D=oV? 


COMMINZSTILKLZ A, IMASTE KZ »LBARgL BALSTyLMAXsLMAXST,-ICSTRIG) a M39 
1 MID% MEST, 4MG_yMGM(1 3) gp MSMVy MGM 24 MUSE (13513) 5°B,NIARZNBN(1 3) 
2 NIN DgNRT NIE RV ING INGNSGNGNITONGTSNPLT gNSTOoSNSYM3y9 NS YMG PI 
SOMMON/SIBLK2/ = BETROAZIACCL(13) 5602013) gCC5(13),CCH(13), 

C< (5) ,CO0S31299) ,COSG(299) -69052R(299),005 261299) .I3°RTSDPRI1L, 
F241 (799) 4597(299) ,FPS(299) gy FOU (13, 29 9 F P56 299) F925 (299), 
FI7OPAGD) g=PAC1L 32,4 
GAM(23)4KTePIMA(23) gPINA(23),SINII299) ,SING (299), 4 
SIN23 0299) 4S IN25(4299) 9X Sa XS2gXS3eX Sug KX IG eg XL gp XLP eK _P19XLP2,y 
YL XL LeXL AX LIpXLG,XL5,XL7,STOCNL,STREN2 

SOMMON/SSLK3Z GX C5VySXSTPULG) gHGO(5) gHSOSTPO (14) 

SOMMIMNSSELKGS NV2gVXO(LE7IV XX CLG 7) AVV(LGT), 

1 BAS 1499499 4 NAW(49549) -B39UN (4994 82W (49), TPS (49) TPN (49) 

COMMOMNZTILKRZ ‘IU ,9(361) FP A127 ,23) 

COMMON/SALKLIOS NWBI361),D0WG0361),D0WN( 361), U(351), 

1 U203519,US C351) .V(361I,VB(361) ~VG(351) pW0351) WIC 351), 

2 ¥33(361).NG (361) -WGRC 361) »pWGGI 2751) 

SOMMIN/T3ILKA1Z TM11—9OM11ST »CM1 2, C0M22,CN33 59411 59M115T,0M1250M22, 
4 YYB3—eFML14F 4124 FM22,F M33 

SOMMON /CBLKIG/S NRUSE (23923) eNRO C1 a 2 gC Sache C5yCGyi7,y 
41 ODELX,IELT, W5GG(44) »WB38 (44) ,WG69 (44) ,WSBI (4G), 

2] VRYX(G4IV,VETIGEI,RRIGI,ENKIGG ID SOE NTOG4) NKP, KPC (KG), <9B (46) 

SOMMON/TNOVAS CPTT(5S) sDELTIMs GAMMA(G1) g ICOMP, INOUT, K<ALT 9 KBy 
1 KIOMy KNS pK E27,KOK aK TYPE pNOALL » NCASE gNCHOT sNDBUGSNYASSsNIRIAL 
2 PR(GAN) gPDAY, POP ,PRINT,ZRFR,RTAIAL(S) sTIME,TITLE (20) e TSTOP, 

3 77109) 

COMMON CN10,CN11.CN85CN9,EPB80 (1805) »sEPROST(L380) ,ETT EXT, 
EXY¥ gIXXSTIVINZ(2) gINZST2C2]) gKSUMA L 3619 gKSUMAS (233) »KY(1805), 
KYST2 (1380) ~NUSE(23,23) oSTT (1895) ,SXT(1805) pSXX(LBI5), 

SY YST2CL3R1) SAAC 3A1) eS2N0361) S34 C361) eSGASH1),SSA(361) 4 
SSACTALD gp S7A C731) SBA C230) gy UCTS LTV VVCL Te 13) HAC L313, i 
KAI ZW) gXGO7 Bd yg XKTT yg XKXT gy XKXM pV KXYSTyX140361)5X24( 351) »X3A(361),4 
XY4A (TALI SXSACTZAL) oe X60 (35199 X74 (2700, XRA(230),7A(29,ZASTRI2), 
7302) p7BST2R C7) ZF LH) ZESTRCLG) 7616) »ZGSTRI14) 9Z4(56) » ZHSTRI1G) 


WF wnanre 


NOVF WNP 


IF (NCALLSEQ.0) CALL PRESS 


T=1 

090 139 M=1_M% 

90 190 N=igMs 

TF (4U5STOUNe MP0 FI69) GN TO LAN 
UU(IN,“)d=X¥CT) 
VVCMNSMIZKXXOMSMUGTY 
WHON,M)= YX OMG24T) 
I=T+1 

CONTINUE 

K=) 

00 £79 T2i,NsT 

Pp = 7 
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PEA pect es GR hhectE A cath eRe E i Relea Nea ak a ek a ll cs ta 


EF UChr yet} EY = 4 

IF (TeFO.MIS2. NNNSYYGL FQ) IT 
D9 FAN yet,Nar 

IF CNUS=(€ Je T?.F O29) GO Fe Soin 
NBC = N3USEC),T) 

JU = 9 

BE $206 Ret) 25 54 

IF (JeEQ-NSA2. AND NSYMS.ES.1) JJ 
K=K +t 

SS Y=T. 

SS2=7.8 

SS$3=9.9 

$S4=13.59 

SS5=17.9 

SS5=).9 

SS7=9.9 

$$3=2.9 

$S9=9.0 

SSinm=0.n 

$$11=0.0 

$S$1?7=0.9 

Ss13 de 

SSt& iP 

SSi5 
SS16 Ne 

NO GIG YW=1—M5 
MM= (MeL) ¥NGT#CT 
SM=STNGCM™) 
CM=TISGOCMM) 
SM2Z2=SIN?G(M) 
CM2=39S26 (MM) 
FPL (MM) 
FP2 (MM) 
EPZCMM) 

iD] 


+e 


TF C872 .6%e8) FR = FOREN, IY) 


“4 
an r 
yon Wo 


Si& = Oe 

00 270 W=tyM? 

IF (MUSEONSM).©O.9) GSO TO 200 
NM= (N-<1) ¥NATOS 

5 = he 

TF tJ Sestah? VS | FOSIN, JIS 
SN=CTNSOCNND 

CN=C3S30NN) 
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"! 


nt are hE sR AR AS IES sa NOI et — 


400 


S57 WAPTEINRS) 


SN2=SIN220NM) 
GN2=395290NN) 

UAN= JIN M) 

VMNEVVOM,™) 
WAN=EAW ON GM) 
SL=SL+UMN*SN? 
S3=S34¢C72 (Nd *UMNECN2 
S4U=SS#V4NESN 
SAESAECT SIND EU NECN 

S7 = WMVFFOG(NND # S7 
SQ = WMN*FOR(NND # SQ 
S11 = WYN*¥EO7ENN)D & S11 
S14 = WYN*TS @ Siu 
CONTINU= 

SS2=S1¥S4SS] 

SS2=31 ¥CCS(M)FCM + SS2 
SS3=S3*5M & SS? 
SSL=SHESM24SSG 
SSS=54*FSC1 0") FCM24+SS5S 
SSH=SA¥SM2ESSS 


SS7 = S7*¥i © SS? 

SSA = S7*T2 + SSB 

SS9 = S3*T1 + SS9 

$S10 = S7*T*% + SS10 

SSii = SLL. + SStt 

Ssi2 = SO*F2? + SSi2 

TF (NSTC EQLNLCINSC.ST.100) 
SS13 = S7*T& + SS13 

S$S1& = SiG*T1 # SS14 

SSiS = Sif*T2e + SStS 


S$S16 = 39*T3 # S$S145 
CONTINU? 

U(KY=SS1 

US¢(K)I=S52 

VRACKI=HSS3 

V(K)=SS% 

VG¢K)=SS5 

VK) =SS6 

WKI=SS7 

WG(K)=S5A 

WACK) =S59 

WOG¢<)=S5S19 

WAIL =SS11 
W690<CP2S812 

IF (NPC CEN MOP NAC STLLO0D 
NPC = JASSENAS) 

IF €1Te£9299 GN TM &SN 
WGOG(INS*Y = S81? 
WGOPBIN@C) = $S15 

30 T7 5A 


SSi4% 
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39 Te &OP 


$9 TO 500. 


j 
; 


500 


COMPUTE ST2ATNS ON STRESSES 


510 
550 


Sar 


LUGF)) 


_WGRF 


WOSPENAT) = 335145 
CORVENUE 


X=0 

KSTP=0 

OO FOO EstansrT 

QO 7IN JsieNRT 

IF CNUSE (Se TPA E709) GO TD 708 

K=Kel 

IMASTPR=N 

IF «NSTZ.FQ.9) GO TI 55N 

990 Sin LSTP=t gNSTP 

IF (LOCSTRELSTR) NEW) GO TO 519 

IMacTe=1 

KSTP=KSTP eT 

SONTINUSE 

UF=U(K) 

UGF=J504%) 

UBF=U3(K) 

VF=VtK) 

VGF=NG(<) 

VBF=VB(K) 

WE=WEKS 

WGOFHYG(<) 

WRF =AB(K) 

OWGF=DWS(K) 

QOWRF=0H3 (K} 

EXX=¥LA¥ CUGFEYLL* CAGFENWGFE 0 5 * CHG ERE QEVGFEEZEUGE FE?) DD 
ETT=XS¥(VRF EXIF CWRE FOU RE 4 SFC NSE FFP YRE SEP EUQE EDD) Y 
EXTHXIFISFOY_ 1 CVG F(T. DEK SEV BED XSF CHUGE FCW SE ENW RED ENWGF EWBF EUBF® 


AC = YJ®VBF + XLI*UGF + 1.0 

WGGF AGGCK) 

WBRF NBRBCK) 

AGACK) 

XL7*WSGF AC 

KKTT ¥JP7#WAAF FAC 

XKXT XJSPWZIF FAC 

IF (IMASTP.F9. 1) GO TO 580 

EXYSTR ££ FRY = UXLEFE2Z9S  SOuce fees 

MKMXST = XKXK = YL7*WGGE ¥X JFVBF 

IF tNPLT.FQ.99 GA TO END 

EFTH=TTAWE* (1. MEX IEV IF HKD SEWED EV CY PBF OD SEV) 
EXT =EXT#XL 1* (WSF PVF A@VGF FWP) 
KKTTEYKTT eX SP VIFOXL LS UG Fear 

KKXT=XKXTeXLiFVGF 

AC = YJ#UXF & WE 

Y¥KYY = ¥YKXYK = YL7FWSGF EWE 

TF CIMASTQALNTL UD YKYKET = XKYXST = YL 7*¥HGGE EME 


xKYY 
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XKTT = ¥KTT + XJFVBEFOXSFV3IF = WF) = XJ2FWEFWBRE © 
1 (CXL*VRF = WED **? &€ ACFECAS + XJ* WRF) 

XKXT = WKYT = XJS*WGRFEYNF + XL 3*NGFEAC 

TF (NDF RAV.EC.23 GO TN HbN 

S144) CMLL¥EXY ¢ FMLZFETT © EMLLFXKXX + FML2*XKTT 
S2A<) CMOZFETT + CML2*EXX + FM22*XKTT ¢ FMLO*XKXK 
S3A(<) CM33¥*EYT + FM33*XKXT 

SHACK) QMALFXKYX & OMIZEXKTT + FMLLFEXX # FML2*ETT 
SSA(K<) OM?2¥*XKTT #€ DML12*X<XX # FM22*ETT + FM1L2*EXX 
S6A(K) QMSR*¥XKXT + FMSSFEXT 

IF (IMASTP.EIL 0) GO TO 66N 

S7ACKSTRI= CMLISTFFEXXSTR 

S8A(KSTRI= OMLLIST*XKXYST 

50 T2 659 

TALL SISMA (JIeTsKykKSTR) 

XL 404) cXy 

x20 64) FTT 

X3A (4) FYT 

X4A(K) XK XX 

X5A(<) Y¥KTT 

XHA(M) XKYT 

IF (IMAST®.EQ20) GO TO 700 

X7ATKSTRY = =FXXST2 

KRATKST29 = YKYXST 

SONTINUE 

IF (KER20GT29) SO TD 2°90 


CT | 


K = 9 

<STrR =f 

00 750 T=1,N35T 

90 750 J=1eN3T 

IF {NUS=(3,I).F 0.9) 560 TO 750 
KS Bet 

IMASTR = 0 

IF (NST2.F0.99 GO TO 740 

20 7725 LSTPR=1,NSTP 

TF (LOCSTRILSTP).NEWJ) GO TO 775 
IMAST? = 4 

KSTP = KST? + 4 

CONTINU= 

IF (NRUSE(SeTV NF OM) CALL PEIT (Ly Sak KSTO} 
CONTINUE 


IF (NCALL~EN.1) GO TO 900 

KZ=0 

TF CKDAML Te Ze ANDCKCLEN1N) KZ = 1 
IF (2PINTLEOL%.) 5D TO 800 

IF {(TIMELTDPCTY 5D TO 800 
K7=K7+2 

IF CCDAMSLTL9Y KZ = 3 
JPRT=NPIATENPI TY 


P2INT PESULTS ANO/JOR CHTEK MAYTMUMS 


880 TF «INDFAV.FQO.1) CALL LISTS 


agn 


910 


TF (NOFRVOEQ.?) CALL LIST2 
TE (KZ er Meteo I®.KZ-E0.3) KC = 0 
C = KC # 1 


I7=f 

90 919 TRP=1,43 
SRUCTRY=A0.N 
PAWI(TPV=N LN 

JO 9196 TS=1449 
AAUCI?,TS=9.0 
ARW(T2,TSIV=NA 
IN POND T9=1—MG 
MMO=( TOL ENGT 
N99 2908 TS=1.483 
IF (4USECTIS,&T929~-5%.9) GA TD 2009 
NNO=CTS<1)*N°T 
IZ=17+1 


KSTR=9 

99 1790 T=1NGT 
YM=M4Oe! 
SM=SING(MM) 
SM=CISGUMM) 
SM2=SIN2G(M™) 
SM2=359S260M™9 


Ti = FPL (MM) 
T2 = FP2 (MM) 
T3 = FP3(MM) 
SU = 0. 

SV = 0.6 

Sw = 06 
SUST2=0.0 


SWST?=0.0 

POLM = DTMACT) 

90 1599 J=teieT 

TF (NUSE CS, TPF IU-NY GO TC 1699 
RK = 4 FT 

IF (NUSE OCS, TI 297619 69 TC 1699 
IMASTQ@=f 

IF CNST2.FQ.9) GO TD 1909 

BS G2 ‘LeieVs FR 

IF LOCSTPIL I eNSe dS) SO TO 928 
IMASTR=4 


177 


KSTR=<ST 244 
920 SONTINUF 


7909 PRLN = JIMA) 
PRLNST = STRIN1*PARLN 
NNENND +9 
SN=SINAICNND 
SNESCSACNN) 
SN2=STN23 0NN) 
SN2=50S27R CNN) 
UF=U(K) 

UGE=SIG(¢) 
USFHIA(<K) 

VF=\yi4) 

VGF=HVG5(<) 
VBFoy2(<) 

WF=W(K) 

WGF=NG(0) 

WRF HAAR(K) 

WGGF AGG{K) 

WSBFE W329 (K) 

WGRF AGRCK) 
DWGF=0US (CK) 
JWBF=HDN3 CK) 

IF CNULE9.0) PP F=P(K) 
PU=S4*SN2 
PUGHSTS(T92D FSMFSN2 
PUR=ST2CTS) FSM FIN2 
Dy=sS42"*sn 
PYG=SCLCUIPRETY2FSN 
PYB=TCECTS)FSM2*0N 
PH = TIFFPSINN) 
PWH = T2¥*¥FPS (MSN) 


PWR TLFFPECNND 
PWGCG = TR*FPS(NN) 
PWBR = TI*¥FP7ENN) 


®PWGR = T2*FPHINN) 

PEXXJ=Y¥_1¥PUGE (1.0 +XLI*FUGF) 

PEXYY=XL7¥*VGEFOYS 

PEXXW=¥_7¥*PWS* (WSF ENWCS) 

PETTI=XJ2*UR= ¥PU3 

PFTTV = yy*PvVAR*¥(1.0 + XUV IF) 

PETTH = ¥JP¥FOWRFIWAF + DWP=) 
PEXTU=XI*IPUI* (1. NFXL A FUGF EYL L¥UBEFIUG) 
PEXTV = XLIF(°VG*ttoh & KXJFVIFE) & YI" VGFFPV3) 
PEXTA = YJIFTFTOWZFECWGFEDWGF) + PWG *(WIFEOWRE) ) 


a 
le 


OKYYN=Y_ 7FPWSS . 
PKTTW=VJ2?*OW9A 

PKXTW=VIS*OWSR 

BE = ¥JF VRP + YLT SUG ® 
PKYX) = YLIFXL7*FAGGFFOUG 


PKXXY = XSFYL7FAGGF ROY A 
PKYYH = OKXYA & YL7*PWGGFAD 
PKTTJ = X¥J2¥YLLFENRVFERUG 
OKTTyY = OVE FX I*XJ2FWQQE 
PKTTW = OKTTH # SWARFYJ2PFLD 
PKXTJ = XJIG*XL TENGRFEFOUS 
PKXTY = XLB*¥XJ27FWGAF EM 


OKXTA = PKXTH + ¥JSFIWGIFAS 

IF(CIMAST2.EQ.N) GO TO 1158 

PEYXIJS=7EXXU 

DEYYAS=I=5X¥W 

PKYXIS=*KXXU 

PKXX¥AS=OKXKW 

21150 Si=.% 

Seo. 9 

IF ({NPLT.EQ. 7? GO TO L2n6 
PETTY=S=TTVeOVF(VFtXJFWSF) - KYFWFEFOYR 

PETTW=CETTW = PWC LINEXIFVIF RAF) + KYFVFFPHP 
PEYTV=PZXTV + XLI* CHGS FOV NF FIG) 

PFEXTW=P>XTW # XYLI FIVE FPNG-VGE FPN) 

PKTTS = XLL*¥2UG & PKTTU 

PKTTV=XJG¥*¥PVZ + OKTTV 

PKTTN=PKTTW - OW 

PKXYTV=X_1*PVS & PKXTV 

PKXXW = PKXYN = XL7T#FCOWGGFWE & PW*WOGF) 

IF CIMASTSOSNE SD) PKXKXWS = FKXXKW 

PKTTY = PRKTTV + XJFPYRF(GLEXJFVOF - 3.FWF) + 
1 OVF(SLFYJ*FHIF + 2. FVF) 

PKTTH = OKTTA & XJZFOW SA (HWE + XYLI*UGF) + PH*I2.*WE 
1 B3e®YSFVRE = KIP*WBIFED & FPWAEKIF(GLFXJFWRE + 3. FVD 
PKXTV = OKXTVY + XL3*PV*WGE 

PKXTW = PKXTWH = XJS*(PHGB¥AF #€ PWF¥WGQF = PW3*WGF) & 
1 PHIFEYLZFCKXJ*WIE & VE) 


A 


Si = DKI(K) + WE 
S2 = VF 
L2OG FU = XLOPSPPeeo UF INGE + HWSE >} 
PY = XLOL*PPOFOVECXIFONSFE & DWAS) + S2) 
PWS SUP PS CCPPONe (St = XELSUGS = XUFVEr os Seu 


IF (NGFRVSFO.1) GO TD 178¢ 
JT = LRAR* (K-11 

KSUM = <SUMATK) 

JIST® = LRAPLT*(CKSTO=“43 
KSUMST = KSUYAS(KSTR) 

IF (“SUMOLFeLAARY GO TO 1300 


1279 Gi = SiS {K) 
G2 = $?4(KY 
G3 = STK) 
GH = SGA{KY 
G65 = S5A(K) 
GH = SAS UK) 


F4 


POYYU*RY & PE TTUFEG? # PEK TYF GZ 


Sa A Se et 


H 
F2 = PKyYUtGh + PKTTUEGS # OKXTUFGA 
| FR = OFVYYFGL # CFTTV*G2 + PEXTV*G32 
Fh = OKYVUFGY + OKTTV¥GS + 2XXTV*GS 
FR = PFYYW#51 # CFTTHWES? # PEXTH¥S3 
1 SHR = PKYYWEGK &€ OKTTHFGS & PKK TWEGSH 
FU = CNLO*FFL & DNILFF2 
Fy = C\IN*¥F3 + ONLLFFG 
; FW = TNIMFFS + SNLLEES 
GO 12 1419 
| 1799 TOTUY=".9 
j TOTV4=5.°0 
TOTW’ = 0.0 
TOTUR = 0.9 
TOTV3 = 9.0 
TOTW? = ALN 
DO 1499 <K=21,L R42 
L = JI + KK 
$1 = HS) (KK) 
S2 = CX(KK)*51 
Gi = SYY{L) 
G2 = STT(L) 
63 = SYT(L) 
TOTUM = TOTUM + S1¥*CPEXXUFSL * CETTUFG2 & PEXTU*G3) 
TOTU3 = TOTUR * S2¥*{PKXKXUFGL & PKTTUFG2 & OKXTU*G3) 
TOTV4 = TOTYY # S1*®CPEXXV#SL + PETTVEG2 + OFXTV*G3) 
TOTVR = TOTVS + S2¥*(OKXXVFSL + PKTTV¥G? + PKXTV*G3) 
TOTWY = TOTWY + SL¥*¥C(PEYXW*SG1 # PETTW*G? # PEXTW#G3) | 
TOTW3 = TOTW3 # S2*{PKYXW*351 # PKTTW*G2 # 2<XTW*G3) “44 
1400 CONTINU? E 
FU = CNAI*TOTIM + CNOFTOTYUR 
FV = CN@*TOTV™ + CNO*TOTVE 
FW = CN@*®TOTWM # CNO*TOTHA 
1410 IF (IMASTOLEQ.%) GO TO 1509 
IF (NOF2V.EO.1) SO TO 1420 
IF (<SU4STELTeLSBAPST) GO TO 1450 
1420 G7 = S7A(KSTR) 
G8 = S83(KST?) 
F7 = PEyyUS * 67 
FR = OX¥YUS * 58 
FQ = PryYYWS * 57 


FO = PKXYXWS * GS 


FUST2 = CNINFF7 + ONILFFB 
FWST> = CNiN*F9 + CNLL*FO 
GO T) 15° 
1459 TOTUMS = 1.20 
TOTW™S = 9,0 
TOTUIS = 060 
TOTH3IS = 0.0 
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1479 


? 


15949 


1A0N 


1900 
41950 


2909 


N99 1470 KKK=1,LRAQST 
LSt® = JISTR + KKK 

$3 = HSISTIA(KKK) 

S& = GXSTR(KLH) FSR 

G& = SXY¥STRILST9) 


TOTUYS = TOTIMS + S3I¥FPFYYUS*S& 
TOTURS = TOTUPS + SU*Ouyvusesy 
TOTWYS = TOTAY“S ¢+¢ S3*PFYYWSF54 
TOYWIS = TATH2AS #* SGFOKYYHS*54 
FUS TF? CN8*TOTUMS # CNO*TITUZRS 


FWST2 = CNS*TOTHMS # CNI*TOTA3S 


= SU + (CES ££ -PUPFPSEN 
= sv + {Fu + DY), FOPLN 
Sw = SW + (FN + OW) FPRL 
CIMASTPLZED. 7) GO TO LHR) 
SUSTS = SUSI & FUSTRSPREMSE 
SWST2 = SWST? + FWSTR¥PSLNST 
CONTINUE 
SURS = SURS + PR LM*SY 
SVRS = SVPS + PRLM¥FSY 
SWRS = SWPS + PRLM*SW 
IF (TMASTRLEQ. 9) GO TO 1709 
SUPSST = SUPSST + PRLNFSUSTR 
SW2SST = SWRSST & PPLM*SUSTR 
CONTINUE 
IF (ABSTSWRS#SHWISST) GT 120630) GO TI 2159 
AAUITIZ,3Z¥ = 1-29 + AAU(TZ,12Z} 
EZW = MOMB2 * FZ 
AQW(T7,TZ) = 129 # AAW(TZ,1Z9 
BRU(TZ) = =<(5U2S5 + SUOPSST) FCKU1) FCKC2) 
YY(MSMS4T7) = -SVRPS*CKIS) FCKIG) 
3BW(TZ) = =-€5SWRS # SWRSST) ¥OKC5) *CK(5) 
TF (NSTPLEG. IS GO TO 1950 
ISUP = (TR-1) #48 
90 1399 ISTI=1L1,NSTP 
ISURL = LOCSTICTSTR) + TSUB 
SIRCIN = STPUNZ © PSCISUS.) 
99 1300 IN=1,48 
SUS? = £0 +- F508 
STPAID = STOTINFEOS (TSUR2) 
AAUCT7Z TI) = AAUOTZ,I199 *& STFADH 
QAW(T7,19) = AAW(I7,TQ) + STRAOD 
50 TJ 2500 
vVV(TZ) = BRUIT7) 
YVQT7W) = RAwWUCT7Y 
SONTINUE 
TF (NSTZOFU.4) GO TO 200 
MAYDIM = IZW = MGMA2 
SALhL SOLVE COA, AAYDNI% 49,9, TFN I. F983) 
CALL SOLVE CAAH pMAVOT Me G9 5F LPM ee 1997) 


r4 


ust 


: 
| 
] 
| 
| 
| 
] 


po 2190 YT=i,4AxOIM 
YYtT) = BRUIT 
YYC(T#eMGMP2) = BRAWCT) 
30 ¥2 2200 


KFO> = 14 
ARITE (552151) 
(ZTHOSOLUTION DIVERGING IN DEPROSP ) 


FORMAY (31 
RETUON 


END 


SUB9IUTINE MSETH 


COMMIN/TALKIZ Ay TMAST 2 KZ LIAQZLRAOST,LYAY,LMAYST, _ ITSTRIU5) o M3, 
1 MINA YMBST2gMGyMGMCA ZY gMS4%B,MGMIPSYUSE( 175139 gNReVRAPENBN(1 3) g 
? NINO SNRT ANTS QVING NON NGNSTANGTSNPLT NST? gNSYM3,NSYMS oPI 
COMBON/SAELK AS AE T2 (23),501113),C00201 3) ,CC5(139,CCS(13), 

CK OKY CAST EZ?IIY,COSG1299),508S29(299),C0826(299), Y2RT~NPRT1,y 
F249 0299) 47576299) gFOR( 299), F724 (13452) 5 FP5(299) 5F725(299), 
FI7(7279Q),=-PH(13,2)y 
GAMO23),KS,°IMA(23),PTNA(239,SINI(299),SING(299), 
SEIN230299),STN260299) g¥SeXS 2% S39 VS XS59XLSKXLPYYL P14 KXLPZ, 
Kee KL De KL Oe KL 3g KL Ge XL OG XL 7eSTAICN1L STAN 

COMMIN/OBLK 3/7 3¥ (62 65XSTRI1L4) gHCGO(H) ~HCOSTIACLGY 

TOMMINZZSALKGLS NY 2,VXYOCL UTD ,YXX (AG 7VCYYULG7), 

L E4U (49,49, AAW (69,69) ,FISIGEI), SIANI49) ,TP5f4 9D. TPAL49) 
SOMMONSSPLKGS FM (RD GEPO LE K7) gH G0172513) pH*C20) sYOUT(159) » 
+ NIUT C1569) ,FHIMNGB) -UL013 513) VIL (1L3,139,nW1013,13)9 
SOMMINZITBRLK7/ RS T2A,CN1eCN12,0N13,CN2,EN25T2, CN CNG eC NS CNG 2CN7, 
| 1 FLeFPZFPO,e PP, FPASTOgH,ASTR yg IE TOST pJEI PST ZIST OFT LS gL CMAX 

| 2 LOMAXS,LOSTPO NEL »SIG0,SIGI2,TNUyTNUSQ 

SOMMIMNSTALKBZ NU e270 351) ,PN(72F, 23) 

SOMMON/STAL KOZ 3TL(8) e3XL(A8),3IXLST(9) ~CORITCSIY,SINSI(S),ETIBD, 

i FX (3) eGYT(AIVNLZ(15) sNREG INTEC OQNZP,SAC(R) gSATI3) gSMAXK, 

2 TSRIT C8 Ys THNUCB) g TMAX XK KNUC 8) 4 706468) gZCSFRE2) 

COMMIN/TBLKIVS IWR (CZ6L),0WSI361) ~OWOC351) , U(I351), 

4 UP0351),US C3519, V (361), VB 1351) ,VGI351) gWI351),N30351), 

2 WISTZHL)d SWGTIEL) WSR 361) WCET 351) 

SOMMON/SSBLKALZS TMLL,CM1LIST »OM1 290 M22,CMS7eOM1L1,0M115T,0M12,0N22, 

P 1 M4ASSsF MLL, F412¢FM22,F 433 

COMMONSIPLK13/4 IOC,ECeEPSIF sGSyHSOFNL pNNOUT22HO, THET AD 
COMMON SOBLKL&/S NAUSE (23423) aNROCeC1 eS 2eCS Cae C5 CG S77, 
1 NELK »IELT, W666 (4G) _,HB89R8 (G44) gHGGP OGG) gHSRWGRD, 
2 VRX(G4) eVPTOG4GY -R2ATKIZENXCRG) ENT OGG) WNKP, KPC (G6), XPR(4G5) 
SOMMINZSNOVAS SPITCS) pDELT U4, GAMMACEI DV, TCIM, INOUT, KALT, KS, 
1 KIAMSKNS XE 27,KNK SK TYOF NCALL, NCASF gNCHPT gNDQUG NYASS,NTIRIAL, 
2 PRIGO) gPOAM, PAP, PRINT, PFIAARTQIAL(S) »sTIME,TITLE (29) »,TSTOP, 
3 271493) 

COMMOM TN1L90,°N21,0N8,60N9,F98 911805) ,EPROSTIL3S90) ,ETT EXT, 
EXX = ¥XSTIZTYZ O79 gTNZSTEO2Z) gKSUMA C351) gXSUMAS (23)) e-KY(1805), 
KYST2 013819) sNUSE (23,23) gSTT (1995) »SXT(1805) »SXX(1805), 

SeY¥STOE (1299) S140 361) gS AREZAL) pSVAC351) oe SGA (TH10e SSACS510 

SSA CTSLV STA C2ZVV SRA 2 TO) ge UUTLIg 1 SP eV CLF e138) NACL LSD, 
XVOPTW) gXGO7 TZ) XK TT XRT KKM KA XKEKKS TA KLACI5 1) yy X200 361) ,X3N (361), 
YEA CTBL) gXSAC THAD YAO (351) X74 C230) pXAACPF3Z9V,ZAC2) ~ZASTRI2ZI, 
Z202) yZBSTR CAV GZF CSV, ZFESTICLG) » 7G, 7GSTPCLG) g74(6V,ZHSTRIL4G)D 


VF wwnvere 


NOW FE WNP 


TNPUT DATA 


C2428 


READ (65,7090) “GyMB_,4BAQ NIA, LAA 
READY (5,790) (MGMCOT), Tai, 4S) 
PEAND (5.7099) (NANUI), T=1, 499 
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AD-A053 954 FLORIDA UNIV EGLIN AFB GRADUATE ENGINEERING CENTER F/G 20/11 
STUDIES ON THE FAILURE OF STIFFENED CYLINDRICAL SHELLS SUBJECTE=-ETC(U) 
DEC 77 C A ROSSe R L SIERAKOWSKI AFOSR-77-3237 
UNCLASSIFIED AFOSR=-TR=-78-0697 


33 


lee 


2 = 2. 
a = 4.8 


ES Ee 


NATIONAL BUREAU OF STANDARDS 
MIGROCOPY RESOLUTION TEST CHART 


a ta ns sng ts 


BY 
90 


150 


159 


199 


Gi, “PRE 


1 


PECAN (5,790) NSYMG,NSYMRB 

READ (5,790) MASTP, LAAPST 

3STR = 9.0 

HSTP = 1.0 

NSTR = MESTP 

IF (NST2eFVN) 561 TO 4D 

IF (NSY¥3.EN. NY NSTR = MBASTQ2/2 & MODN(MASTR,?2) 
READ (5,71099 ASTEAYST? 

READY (65,7099) (LICSTOCT) »sT=1,NSTR) 
READ (5,79009 NALT,NSNIZND ERY 

PEAN (5,70990) NNOQUT 

IF (NNOJTLEQN.0) GO TO 79 

CO 593 [=1,NNOUT 

READ (5.7900) “IJUTCTIY,NOUT(T) 

READ (5,7990) NKP 

IF (NKPLEQ.9) GO TO 9N 

DO &Y T=1_NKP 

RFAD €5,70009 KOG(I),K9PB8(T) 

IF (<TAMVLENS 1. AND CKTYPFLEQ.1L) NDF PY=2 
IF (COAMCEDe LP ANDCKTYPELFO.3) NOEPV = 
IF (NIERVALEDe1) LEAR = 1 

READ (5,79009 NL 

IF (KTVYIELTeSY NL 3 

TF (<TYPELTC3) NL 1 

REAM (55,7100) XLOgTHFTAD SA 

IF (NPLT.EQ.7) A=1.0 

IF (NDE2VeFENe1) GO TO 150 

00 130 I=4eNL 

READ €5,7100)9 HMCI),RHOM(T) gEMCT) 
READ (54,7100) TNU,ZSIGOZEP,EPSIF 

Go TT) 190 

00 159 T=1,NL 

READ (5,71390) HUT) ,PH9OM(T) 


2 


READ (35,7100) ©XCT) ZETCID » XXNICI) pTHNUCT) »S¥TUT) 


READ (5,714903 SATCI)D,SAC(TI) 

IF (<TYPFELNE SSL. ANDKTYPE NEG) 69 TD 199 

IF (<MO4M.NE1) GO TO 190 

READ (557199) FT,GC,9° 

REAC (547109) CLFGCI, Sd, T=1943)—I=1-%G) 

RFAND (5,7100) DJELTIM,TSTOP,PRINT 

IF (INOUT.FEFQ.9%) GO TD 2100 

NT OJT THE INPUT 

WRITE (657170) 

WRIT= (597209) MGR MBAR, NBA L342 

WETT= (597210) (MGMCID »T=145MS) 

WRITE (5572299 (NENCI),T=15"3)9 

WRITS (5,7507) MRSTR 

TF (NST26GT29) WRITE €6,751029 BSTP,HSTO,LBAPST, 
C(LOSSTPIT) »T=1,NSTS) 

WEIT= (5972759 NSY4GeNSYMANPL TeNBNDYNDERV 

WRIT= (5571599 NNOUT 


184 


1169 


Lisn 


1199 


“2490 


IF (NNOUTSGT.9) WRITE(647189) CMOUTET) NOUTCID »T=1,NNOUT) 
HRTT= (447185) NKP 

IF (NK2.5T. MY WATTE (647189) CKPGITY,KPBCT) eT=1—9 NK 
WRIT? (4,11R99) NLYYLP 

IF (NPT F009) WITTE (6,723) THET AN 

IF (NPLT SEQ.) HOTTE (697259) THETAN,A 

IF (NDE2NCEO.2) FO TO 1120 

DO 1149 T=t9Nt 

WRITE (C5g1179M) TyHMOT) »PHOMET) »EXCTY SET CID» XXNUCT) » THNULID, 
1 GYT (TY, SAT(T) SAC(T) 

IF (CTYPE NEST. ANNKTYOFWN=.4) GO TO 1190 

IF (<724.NF.1) GO TO 1190 

WRITS (5911990) FC,SC,2C 

50 7) 1999 

WRIT? (557289) CHMCT) »PHOMETO SENET) g Tate NL 

WRTTS (547309) TNU,STGO,FP,=PSTIF 


APTTE (Ay 7H499Y COPGCL 9d) Tete 43) p JEL MG) 
WRITE (543200) MFLTEM,TSTO 220 NT 

T=0 

MGMR= 


DO 21553 M=1,4G 

MM = MGY(™M) 

NO 2150 N=t,¥8 

NN = “NON ON) 

MUSF ON YD HL 

IF (Le FQeNNCIT) GO TO 2130 
OO 2149 Jal INCUT 

IF CYMCE ON MOUT CID CANDSNNOFIENOUT( UD) GO TO 2120 
CONTINUE 

GO FD 2LGt 

MUSE(N,%I9=0 

IT=Te#1 

SO ¥7 2150 

MGMBR=aM54B4L 

CONTI “UE 

MGMR2=?7*MEMA 

PO 22N0 M=1,% 

MM = MGYCM) 

CC1i(4) = MM 

CCS5(4) = MM + 1 
KIJSLIC LA STHETADN 

IF (NPLT CENT) XSJ=PT/STHETAN 
DO 2320 N=i,¥2 

NN = N&ANCNY 

ec2(4) = NN 

CCA(N) = NN #& 1 

RETU2"! 


FOOMAT(ST12) 
FORMST(5°1201) 
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ree 


7150 FORMAT {QHONVOUT = [4 

7170 FOOMAT C25HLTNOPUT JATA FO® JEPROSO 4 
71RN FORMAT (2TG) 

7185 FORMAT (7HANKP = TRY 


779M FORMAT ( LOHOMG = I12719H M8 =I 
127104 VAP = [2/7104 NSAP = I2/194 LAR = [?7) 
7210 FORMAT (19HAYGM = (1075)) 
7229 FORMAT (CLSHONR' = (4f735)) 
7225 FORMAT CANHONSYYG = T2714 NSYM3 = T2/71L0H0NPLT = T2/ : 
3 10H NOND = [2/7 10H NOEQV = I?) . 
7739 FOPMAT(L7ZH THE TAN, IN = £16.48) ‘ 
7257 FORMAT(17H THE TAN, PFC = F15.A/17H A, IN => €16,8) 
7289 FOOMAY (17H0 HMg INg &X p2LHRHOM, LB-SECFHEASING FH, aX, 
1 THEM, PST Z(3517239)) 
73PN FORMAT (LZ7HOTNU = F16.8/17H SIGO, PST = £E16.8/17H E°, 
4 PST = F15.8/17H EPSIF, IN/ZIN = E16.8) 
7400 FORMATISHOFG = /{56£14.5)) 
7500 FORMAT(LOHOMIST? = [2) 
7510 FOOMAT(1 7H ASTP, IN = £15.8/17H HSTR, IN > Fl6237 
1 104 LRARST = I5/10H LOCSTR = 615) 
P2N0 FORMAT(LSHOMEL TI“, SEC = FLB-8/15H TSTOP, SEC = £15.8/15H PRINT 
1 = =1628%) 
Zi FCO FORMAT CLOHANL = [2/717H0XL, IN = £15.28) 
11799 FORMAT (6HALAVERT7/27H HM, IN = €15.537 , 
4. 27H RHOM, L3-SFC#FQ/INF*G = £15.87 i 
a i ET = £16.8/ 
3 274 FT, OSI = F1668/ j 
4q 274 X XN!) = €16.8/ 
5 274 THNU = £16.87 
6 274 GXT, PSI = F16e87 ' 7 
7 274 SLT, PSI = £16.87 | 4 
8 274 Sac, PST = £16.8) | 
11990 FORMAT (14HO5EC, PST = F16e8/711H GC, PST = E1487 
i 4i1f4 98>. IN = £16.81 
ENO 


ae steerer coemenn —-- -— 
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1 
1 


1 
2 


1 
2 


1 
2 


1 


1 


NOovrwnr 


SUQEVITTNE ASE T2 


COMMINZSTPLKIZS Ay IMAST> 4 #7 ,L 90%, L PARST yLMAXY gL MAXST,. ICST2IG) oMBy 
SVADMAST I, MG MGMED BD a MSMDg MOMAD,MUSE (13913) 9 NByN IAP NANG13) 4 
NINQe MAT, IIT PU ANG ANONS, NONI TONG ys YPLT ey NST Re NSYM39 VS YMG PI 

COMMINZSALK2Z «= BFT2C230 oe SS1013) gSC201 3) SSCOCSUL SD COAL 30» 

CK CAD yg PISTE799) COSC (299) gCOSZPI739) ¥COS2602799),I2ATyIPRT1, 
F4£299¥ 45°57 ( 299d gFO31299) y FPG (135 299 F 951299) ,FP5 (299), 
F27(29DV FPACL eA, 

GAMI2 3) KT y FTMALP3) pPINA( 23) gSTINI(299) ySING (299) 4 

SIN 2302990 4S TM2G C299) yX SoX¥ SPX I3SeK SG eX ST oXL oXLPAK_PLgXLP2, 
YIP TeYLDeVL Pe¥LIyVL Gy ML EyXL 7S TION1g STRONP 

TIMMONSS SILKS G15) pGYSTE (14) pHODIAD gHGOSTPI1G) 

SOMMINSTALKGS NY 2eVXOCLG7I XY (CLG 7IGVYCLG7D 
LAU UGI GO 4 SBWIG9G49) pPIIOG ID, WWIGID G IPSI49V TON (G9) 

COMMOMNATIL KGS EMECRD = O2(L47) gE GCL% 91 3) HME2T) pMOUTOILB ADs 
NIT O169) pPHOM (BD ULCL 3413) V1 (135139 gWLE13 51 3) 

SOMMINSSRLKIS —BSTIgTND ae CNL Ze ONL oCN2gCN2STR CNS gCVG ONS yCNGsON7y 
FL yFP,EOO,E PM, EPPSTI, Hy, 4ST TE IU Te IFIUST y ISTOFT SLO LL CMAX » 
LOMAX Sg LOST 2g NELO,SIGCO, SIGO2y TNUy TNUSO 

SOMMONZTALKRAS = NU, 7203619 498(25 42%) 

COMMONSSOLKGS «= STLEBD GS BWLE9D SXLST(BD pCCRITCAD gCINSE (3 SETEBD 
FX (Ads GXTIRD »NLZE1LA) yg NREGeNTEC D9N7PySAC(B) oe SAT(3) yp SMAXg 
TOPITURD »y THNUCAD gp TMAX yXXNUC 9), ZOCL0 eZCSTRI2) 

COMMON/TPLK1NZ OW9(3H1)e OWS0351)5 HOU 3H1), UC351), 

UPCISALD gS UTCTSAV gS VOSSAV  VICZOL) pVGC 351) pW 361) pNBC 35104 
WIB(361) AG (3610 oe NGA 351) eWGGl 61) 

COMMINSTALKILZS SMALsOM11STyOM1 29 0M224CM33 99411 »OMLIST,OM12,0M22, 
QMS, FMALyFML2,FM22,F M23 

SOMMONSSALK137 IGeECeEOSTFE Gly HBA NL yNNOUTs 2HOy THETA 

COMMON/SNOVAS = SRITESD eS DELTIM, GAMMA(41) 4p ICOM, INOUT, KALT 9 KBy 
KAM, KOS y KE P> ye KOK yKTYPE gNCALig NCASE yNCH?T ys NOSUG,NMYASS,NTRIAL, 

PRIGM) gPNAY%, P9D4PRTINT, OF MVATATAL Loe go TIMESTITLE (200, TSTOP, 
77109) 

COMMON CLD ANA oON3 9 OND, FORT(1895) pEPROST(L 320) SETT gEXT, 
EXK gE XXSTI2g TN7C2) gINZST 2029 eK SUMA (351) oe KSUMAS (239g KV 16050, 
KYSTICLISV) ~NUSEC23 423) pSTT (1895) op SXTC1L905) pSXXUL BOSD» 
SXXSTVEL AWAY -SLAC TEA) pO 2NCZHL) pSBOCSALV SGA (3H1) 4 SSAC3EL) » 

SAA CTALY yS7AO73OD yg S9N (239) 4g UCI Tot 3d op VVES F913) WALL Tol 3) 
YVC2TV eV GCP Bd g MKT Tg XKYTS XKX Xp XKXYST yXIAC3H10 9 X2NC 361) X3A (3610, 
YUP (BBL eX SACBEL) eK EA (35104 X70 (2TDVSXSAL23NV, ZA(2) yZASTR&I2), 
2302) g7BST2 02) gZFCEIg 7FSTICIG) 9 7G(6) pZGSTRUIG) pZ4C5 Dy ZHSTULG) 


IF (NDERVEFO.1) SO TO 271° 

CALL LFSFND 

HSTP2 = HST * 0.5 

Z8STOCLY = 4S TPP 

7BST202) = HSTP2 

TF (LBAP.FOLNY GN TO B00 

IF (<TYSE NEL SY GO FY 2550 

SET J> Z=QUIVALENT LAYER FOR HONFEYCOMS METAL, NDERV = 2e 


187 


] 
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i 
e 


a? 


24UNnh 


2770 


2705 


NEL? s ? 

42=h4( 300.5 
HOL=4y(%)-4M(2) 

HO?=4™(1) 

Z8(LP=-t24N GFN? 

ZA(2) = HP? = 25 *HDL 
SUMP4=9,9 

HO=0.9 

99 2499 T=1,? 

HL=HYCT) 

SUMP 4=SJMPHF?2HO4(T) * (HL-HO) 
HQ=sHi 

NL=1 

EQC=4M (3 Y—HM(27 #HMC1) 
CQH=( HMC TPFHYCAIVHHYCLID ESAT CZ. DERM LDF CHM (32 -H4( 2090 7EOC 
PHO=SUMIAHZE OA 

H=FCY 

HM(4)=H 

FEL=FEM( 1) FEQRZEQH 
EP=FI*FICZENH 
SIGO=SIGO*FQS/EQH 

EM(1) =F 

IF (<9AM.FQe?) GI TO 2705 
TCRIT (AY =FOSTE 

GO 9) 2705 

SINGLE LAYER, NOFRY = 2.6 

H = HM{1L) 

EL=FE4C1) 

RHO=?2HOM(L) 

H2=C.5*4 

ZRt(1) =-42 

7Bt2) = H2 

NELP = ? 

TF (<9 44-19 2499,2709,2705 
DAMASE 

TCOPIT (1) =SIGID 

CEPT (1iy=SIH64 

NFLP = 1 

GO T9 2705 

TCPITT (1) =EPSTF 

IF (XTYPELEQ.1) GO TO 2705 
TCPTT(1) = SIGO 

NELP = 1 

ST = JELTIM 

FI=2A4OFCLN = TNUFFAIVSEL 
F1 = H¥*®P/S (12,0 FF 3) 

F2= 3153/PHO 

CALL OTSTEP (F1L1,F2,F3%,F3,° 1) 
IF CIT .5T ILM) NELTIM = CT 
G9 T) 2759 


188 


cal a EE a a as 


NOERY = L. 
CIMPIJTE HRA, 
a1 


HON th 


HMOUT) 
‘4teFt. - XYYNUCT) FT RN UCTD) 
=v (1) #227 
ET (T)#22? 
YYNUCT) #222 
o¥T{T) 
RxL «Ts = 844 
PTL(I) = B22 
3XLST(I) = EX(T) 
Ht = HO 
HLF¥O oe HO **2 
a R71*H2 
RAVE 
322%)? 
QP2enL 
RWZ*ND 
R3*M1 
R1{?*Q? 
Q{2¥OL 


- 
- 
- 
= 


0 5¥* M1742 
oe SFQO3704 
eo S*AS/AB 
eo SFA7/08 
0 25*(4911 & 43922 + 4HBI3 # HIL2) 


wom nm WwW 


" 
| os | 
ao 


uo Oo tm me on te te a 
yvaowveagaIawvao7tIJu:,. 


2720 


SMLIST = 0. 
OMLAST = 0. 


DD 27279 T=HL)—NL 
Hi = MMOT) 
44ST = CHSTOJNL)*T 
BTYT = RTL (Tl 
XNUS = YYNUCTI*¥3TT 
ayx¥xX = 3AYL CT) 
BxyxKST = SYLSTC(T) 
SXTL = 
PHOL = WHOMI{T) 
Hii 2 4L = 1 
411S7F = HAST 

+ SATCIV*FHLS 

= 741 & VXX*H11 
GM12 = 5M12 *& XNURFHI1L 

= 3M22 + BIT*¥H11 

= 3M33 + GXTLFHLIL 
SM1iST = CML1ST & SYXSTFHLIST 
PHOR2 = PHNRS + PHOLFH1I1 
H42 = 4LF*¥2 - HDF*2 


H22D = 412 = 2. *HAARFHL1 
FM11 = =M1t & BXX¥*H190 
FM12 = ©M12 + YNURFHL2N 
PMZ2 = ©M22 + ATT*HL20 
FM33 = =M33 # GXTL¥HL29 
H12% = ALF¥Z - 4Heey 
Hid = 415 = Se*HRARFHIe + 3. FHRA CES OFTHIS 
UMIL = JMLL € BXXF HD IC 
QMi2 = JRL2 + KNUBFHLI0 
BOM22 = IM22 * BIT*FHL ID 
OSS = IMSS + GXTLFHLSD 
DOM115T = OMLLST + CAXXST*HLIST #3) 74.0 
HO = 41 

OH3 = 1./7HMENL) 

DA = 1474 

O2A2 = 1.702, ¥4**2) 

D3A3 = 16/(13.*A**3) 

SMLi = 3M11¥*94 

GM12 = SMi2*98 

GM2Ze = 2Me2*0 

oM32 = SM33*OA 

CM115T = CMLISTFIA 

FMi2 = =M11*9702 

FMi2 = =M12*9202 

PMZc = =MZErOr Re 

FM33 = -¥33*9742 

Oia Se IMELO IES 

OML2 = IML2*D3A3 
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SS aaAiieeniniraaaaetinniaampentinenate:tcanaenaaieaeei Reena 


2730 


2740 


2745 


2707 
2769 
2ann 
70nn 


= SM22873 43 
YMZZ = PMBBEVIAZ 

MMLIST = AMLISTFOTBAZ 

DHOR2 = PHIB2PEQNY 

RHO = 247NP 

H = 4(NL) 

FRZa=2AC/SIH3 

FLEP AQF A FR R/F 3 

FO=FSE3 

Fuh = © (RFO411) 

FS = aa | TERBEZSEG 
F3=F3/(4FCOM22) 

NY = NELTIM 

CALL ATSTFR (FLO P,F RFU ,E 5S) 
TF tC. Si eteb) OEETIM = OF 
NFLP = 4 

NZ =) 2 

ZCST2C19 = -HSTP4OLS 
ZCSTAC2) = FOSTRILY + HSTO 
IF (<¥VPELEN.S) GN TN 2738 


NLZ(1) = 1 
7001) = -H3AP 
NLZ(2) = 1 


ZC(2) = 7001) + A 

IF CKTYPOELT.3% GNM TO 2745 
NLZE7) = 8 

ZOCi1) = 7OC1d + S¥*HM(1) 
ZO(2) = ZOO? = «S5S*CHM(3) = HM (2)) 
GO 13 2745 

HT = -H2AP 

HS = 4M(1) 

NO 27469 T=1,NL 

NEZE2Z*E ve SP .s Ff 

NLZO2*T) = I 

ZOC2*T = 41) = HT 

IF (©eGfFet? 4S = HMOIT) = H40TH=1) 
HT = HT # HS 

ZOC?*T) = HT 

N7P = 2FNL 

IF ¢(€764%,E0.7) 393 79 2ANN 
DO 2747 TH=1,"L 

TCPTIYTCIT) = SATCT) 

SCRIT(I = SALT) 

CONTINUE 

RETU2" 

KFPP = 4 

RETURN 

END 
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SUBRIUTINE OSFET3 


COMMIN/TSLKLZ Ny IMASTR, KZ LBAR LAYS TL MAX yLMAXST,.ICSTROG) M3 ye 
1 MAA ALMOST IG 4G MGM 3) oY GN3 9 MGM AA MUSEC LIS LIV eNSeNSAPNINEIL3), 
2 NINO PVA T SND AV eNGINGNG, NGN3A Ta NGTy VOLT ae NST Re NSYM3ANS YMG, PI 

SOMMPINZT ALKA = 3ETP(23),CC01(13)2,5C20! 3) gCO5(13) pCCBOUL ID, 

1 CK O05) ,C9S20299) ,C9SG02939),5082773(299),0905256(299),)2RT,OPRT1, 

2 F721 (239) ,=92(799) ,FP31299),F 7% (13, 279,FP5(299) ,FP3(799), 

3 FI7 (> IAI9IEP (17,2) 

3 GAAP 7B) gy KS eg PIYMAZ23) gPINA( 23) y,S INI(299) gSING (299), 

4 SIN279€2799) 4S IN2G 0799) gX SgX S29 XS39X% SG e XIS—9 XL a XLOQX¥LPLyXLP2, 

5 MLO Se KEL eV re XL G—eXL Ge XLS GXLZ7eSTALNL STRONG 

SOMMIN/SALK3S «GY (63 -GXSTPH(14) ~HGO(E) FHGOSTP(14) 

SOMMIN/JTOLKGS = NY2ZAVXOCLG7D XK (147) ,YYULE7), 

1 BAI (L949) eAAWE49,49) 230149) »-RW(49),TFAVC4G9) ,T IV (49) 

COMMINZSALKSS = FM(I8V,EPR (147), F513 41 79 4HM(29) MOUTI(IL469), 

i NIUTCEO9UV,SHOMER) ,U1013513) V1 113513) ,W1013513) 

SOMMIN/TOAOLK7L = ASTRACN1yCNLI2,CNL3 9 GN? gON2AST? CNS o ONG gp CNS gL NG QLN7y 
1 Fue =P oF POZE PFPGEPPSTI GH ISTZ IF IAT gIFIRST eo IJSTAFT SLC yLOMAX, 

2 LOMAXS,LESTASNELPZSTGO, SIGI2], T NU, TNUSQ 

SOMMON/SALKRS = NUS P0761) 524(23,23) 

COMMIN/SALKGS = ATLIS9) eBXLE9),SXLSTI8) ,CCPITIAB) sCINSTUISIGETLBD, 

1 EX 3) ¢GXTO3) gNLZ(14) 5 NREGINTEC ONZE SSAC(R) ySAT(3)9 ~SMAXK, 

2 TOSTIT C8) sTHNUC SS »p TMAX,XXNUL BD, 70040) ,ZC0STPL2) 

COMMONS 3LK197 IRKR(361)—,NDKS(361)_, ONO 361), U(351)9,4 
i U3 035.1),UST751),V0361), VPC 551) pVSCI5L1) HU 361) WIC 361), 

2 WIA (351) WS 1351) WGRC351) pAGGI 351) 

SOMMINSSALKIIZ CMLL —eCM114ST O41 260%224.0M33490%11-0M115T,OM12,0M22, 
1 DAFT EMLL,EML2sEM22,F M32 

COMMIN/TBLKL3S IT SECs PSIF yG5 4 HBAP WNL gNNOUTe2HO, FHETAD 

COMMIN ZITBLKL4Z NAUSEC2A3 23) aNRC C1 se T2eC3 ashe C5e Che s7y 
1 NMFLX,YELT, WSGGUGG9, W882 (44), NGG? (46) -WERB (GG), 

2 VAY IRYD VOITCGGD, PACK, ENXI44) FENT(SG) aNKP4KPG(46) 4 <PRUGG6) 

SOMMIN/STNOVAS SOTT(5) gDELTIM, GAMMA(G1) gICOMOINOUT gKALT KR, 

1 KVL", KIS gS LEP VU KIK A KTYPE sNOALL» NCASE gNGHPT gNOSUGSNYASS,NTRIAL gs 

2 PRIGTY,PONAM, POO,PRINT,PFR,PTATALIS), TIME, TITLE120),TSTOP, 

3 27109) 

COMMON SNIMGENAIL SONS 9 SNGg FPRD (1805) yEPRIST(L380) gETT EXT, 
EYY = ¥YSTI,TYZ (2) gINZSTI02) gKSUMA (361) gKSUMAS (23005 KY (1805), 
KYST2¢1 787) ,YJISE (23,23) ,S5TT (1395), SXTC1805) ySXX(L905), 

SKYXTT ICL SIV) ,SAAC VHLD gS 2O(3H1L) pS3NCTE1) pSUA (361) ,S5A03H1), 

SAB CZALD gS7AC23TV9D SRA (27D), UUC LI LV gVVCL FAL) WAC SoL 36 
YBICPT eg ¥GO7A Ty XKT Ty XKXT yy XKX Ke XKXKST pXLACIOLI» XZ2AC F419 ~XSA( ZBL), 
XA (TOLD XSACZALI XGA (2619, X7N (777) XBAI ZST I, ZA(23,ZASTRI2), 
Z312) g7BSTOU7) pZFLSD gp ZESTPC14) 92606) pZGSTPEC14) pZ406) eZHSTRILG) 


NOPUF AWN 


co PRINTOUT NISCOISFION OF DEPROSP DATA * 


i 
i 
' 
| 


PANT WRITS (5, 9390) 
IF (NPLT.EQe9) APITEIS, OCOD 
IF (NPLT LEO) WOTTE(6,9509) 


7 


pest AN ape COPY 2s: 


sss iacaalamaiaaiaiiaa 


VINA 


tUNY 


zOony 


BO TD (7229, 7A, PBR, PARP, 2AQDN), MYYPH 
APIT= (45,9401) 

SO -%) 123s 

WETFE (5459799) 

E019. 2250 

WRIT (54,9801) 

30 T9 2359 

WRITT (35,9879) 

69 ¥7 235 

ARITS (45,9847) 

IF SRN VEN ete Ve NINDS ENS 3.02, NINNLE%.6) 
TF OMEN Ve FQ ae MQ NBN ED 0b oI 2%0 NINIDE 76 9) 
TE CNPND C7 Fe TVA NBN CED 7002. NANOLWE 2. 3) 
TF OMEN YF 1 ee DQeNINV CFI 04 oI. NINDF 25) 
TF UNAND CE Ae OF eNINN CED 3 2% 2 NSNDCED 7) 
I CNP ED Fa Me NANT F408 97%. NINN EO. 9) 
IF (NOZIVENWLY WRITE (66,1019) 

IF (NDFTPVCENF) WETTE (68,1020) 
W2TTI CA eT AIATY MG eM gM PAP, NWI L VAP 


WRITE 
WRITE 
WRITE 
WRT TE 
WRITE 
WRITE 


IF (NSTRZONESV) WRITE (6e10F10) MIST, LBARST 


WRIT= (5519829) 
90 23279 M21,%S 
MM = MGMUM) 

N00 2370 N=41,¥2 
NN = NONON) 

IF (MUSE (UN ,M).F 720) GO TO 2970 

WRITS (54919839) “MANN 

CONTINGE 

WRITE (4,1498599 Y¥LP 

IF (NPLT.LEQLOD WEITE(5,1090909) THETAN 
IF (NPLTCEQ01) WOTTH (5,110 IN) THET AD,S 
IF (NE 2VCENWPY 57 TI 3580 

WOIT> (46120509 HRARV (Ie T=a19NL) 

WPTT® (2,12200) CHM( TI ,T=4 oNLD 

WPIT> (55127797) (PHOMIT) yTatyNL) 
WOTTE (A,4779N3 CEXOITIY, I=! 


aN 
WRITS (5912G909 CETCTY, I=L ,NL 
WRITS (5922597) CYXNUUI) ,T=1,NL9 
WATTS (A,1?2610) CTHNUCTY TEAL yNL) 
WRITD (€5912765%) (GXT(TY,T=1.NLI 


TF CKTYPS NEST CANTO KIYPE.N= 23) GO TNO S304 


WRIT (5912997) (CSATCTI» T=t5"L9 
AQITS (55973079) CSACCTD,T=16NL) 
$0 79 340N 

WRIVT (5,42700) (SATCTI), T=1,“L) 
WRITS €5,172RF0) (SACL) gs THlyNL) 


TF (KTYPE NE SZ LAND KTYPE NE 4) GO TO 3599 


IF (<N44,NE.1) 30 TO 3500 
WOTTE (5,17190) ©C,60,% 
69 TT 7AAN 


APITE (3911190) HPHDEL Tg STGDZEP,LOSTE 


193 


(6,9390) 
(65,9329) 
(559349) 
(6,995) 
(45,9339) 
(45,19999) 


PTR warns ae 


WOTTE (55122999 CLCFGIN, MD, NEL, M3) MH NG) 
WRIT TOS, 113709) OFLTIY,TSTAP,PPINT 
IF (YMFRVLEQ.1) GO TO 4820 
NO GINO K=1,L3A? 
ZHOKY=5¥ (KY FHF LS 
ZF(KY=74HCK)SA 
ZG CK) = 5K CK) #¥? 
CONTINYS 
ZAC1) = ZBCLISA 
ZO2) = 7ROAISZA 
INZ(1) 1 
IN7 (2) L3A2 
IF (NST22EQ.-9) 397 TO 4120 
00 46919 K=1,L749ST 
ZHST2CK) = GXSTACK) FHSTR*¥OG,.S 
ZEST? («) 7HSTICKI SA 
ZGSTRCKY GYSTOI(K) #¥2 
CONTINU? 
ZAST201) 7AaSTAC1) 7 A 
ZAST207) ZRASTO(?Y SA 
INZST? CL) be: 
INZST®(?) LAAPST 
NNSY'G = 9 
NNSY¥S& = fC 
ITF (NBN CEN Se CAeNIND CAD ea 7eORe NINDLED29) NNSYMG = 1 
IF (NAN FQ 5A DFPONINILGE 28d NNSYMR = 1 
IF (NNSYMGLEQS Le ANDSNSYMGLE909) WRITE (6.13209) 
IF (NNSYMRZET. LeANDSNSYM3,05°Q.0)9 WRITS (6413200) 
NGT = M3AP 
= NRAP 
NGT 
Nat 
(NSYYG.FOQ.1) NG = (MGT+#1972 
(NSYY3.2EN.1) N@ = (NQT#1)972 
ZEMGMNA 
PT JFLOAT (2* (4M8A9-1)) 
PIZELOAT(2¥* (N3A°-1)) 
NSY4YGeEN1) PIM 2e.*PIM 
IF (NSYYRLEQ.1) FIN ?."*PIN 
R=A/4 
XL=XLPs/(OT*O) 
XLI=L)0eN/YL 
¥L2=XL ¥*2 
XL3=70N*¥L1L 
XLG=7.0*YL2 
XL5=¥LGFO 
XL7=1.9/7¥L2 
CMR = X¥L¥¥O2 
CN9 = CVASE2.N* 2) 
TF CNDFSV.EQ.1) GG TD &O4D 
Ci = LeMFCAFFED¥VL FEZ) 


iit pas conc 


a] 
v 
1 


(XJSAY FRIZYL 
CYP AY FEPENY 
WIP CABYL) E#? 
-tiRR OSD 
SA = HERBAL 


(? 
nS 
ow 4 


MFLX = 2TM¥YLPsPT 
MELT = PIN®TAPTANZPT 
NEL T = AFNELT 
CNL0 = 7.0*CWR 

SNL = SNQZ(3.9F9) 


37 ¥9 4950 
BUN CN1O 2. *CNS*R 
oN19 
AME LFOSYL FHT 
(CP4L2 © G6. OFOMSS) FO *X J¥F2 7XL 
™. POBNBY JFED 
(OML2 + &. OF OMSID FAY IFXL F¥2 
Ha NENMZTEYNJFAFEDSYL 
LN5C XJ2 = XJFF? 
XJB=YIFKLI 
XJG=>, OF VS 
XJ5=7. TF YI 
DPRTIL=SPRINT*OEL TIM 
XLPA=*L5 
XLP2=2.9FXL *? 
XLOBZ=1.0/XLP 1 
PRL = t4.7 CRHOSXLO** 2) 
SIMPSON S RULE. 
MBAR AND NAAR MUST BE ONN NUMZEPS FO? FULL PANFL. 
DO 419M T=1_y4NA> 
Fo = I-14 
SAM(I) = FFPTM 
XGCT) = GAMCT) *¥XLPZET 
DIMA(T) = PPL¥*4.*PIM/3. 
IF (NEYYGSEQe1) FIMACT) = .SFPIMACT) 
IP CUTS) 72.—F 0. 142) PIMAITY = 2.°P IMO UT) 
GLiAO CONTINUE 
PIMOCLY = SIMA ULES 
OTMAC(MBAR) = CIMA(MRAR)* 25 
90 &79N Tsi,Nase 
et Sp | 
QETP(T) = FFITN 
YP(T) = BETRIT) FTHETAQ SPT 
PINA(T) = 4. FPIN/3, 
IF CNSY4FeFO01L) PINATT) = .S*OINE (TI 
ITF (CCT#197208&OM. T7727) SINACT) = PINACT) *2. 
L219 SONTINUZ 
PINACL) = oSVFPTNA(1LY 
PIMNAC™IADP) = AINA(NFIAS)Y* 5 
STHCVL aS TORY ZN 
STPCNA=4STOFSTOCMLS (FT FH) | 


Q2 
Zz 

> 

os 
wou 


oO 
Ww 
Woe 8 es 


saa 
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410 


&L2N 


&G4GN 


G5n9 


DO &4PN T=1,4AA2 
NO ANT Jat gVIAA® 

NSUSS (J.T? = 0 

NUSECHe Fr} = 2 

Ir =09 

DO 4430 T=4,NGT 

00 4279 get,int 

IF (TILEDeNK>? 39 TN 4430 
ON G4tt K=1,K? 
IF (TeEVeKPGIK) CAND. J0E RK P8IK)) GO TO 4420 i an 
BONTINUE : 

30 T) G43" 
NUSFEF(CS,T) = 3 


TT = TY * 1 
CONTINUE 

M4 = NMA> - vevu4n 
K So) 

D0 4445 J=2,Ni 
KS <a 


IF (NUS= (Se 42 0 5T2 7) NRUSE(J,1) = -K 

IF ENS YYE CED e ANDONUS® (Jy 4882) .6T22) NBUSS (Je 4349) = -K-Ni41 
CONTTIUR 

Ni = “BBR = NOVMG 

K = 0 

99 &459 IT=2,N1 

K =< #1 

IF (NUSS (1,LT)-6T.2) NSBSUSF(L,I) = K 

IF CNSYYVICENLIT CANT .NUSFINBOAPGTY2GTA22) NBUS={(NBAR, TI = K#eNin-d 
CONTINUE 

NEUS" (4,14) = 191 

IF (NSYYP,FQ.1) NPUSF(IN®AQ,4) = 102 

IF (NSY4GeFO.1)9 NFUSEC1,MBA2) = 173 j 
IF (NSYMGENSYM?.FO.1) NRUSE(NZAR,MBAP) = 1094 

NQC = 1 # NSYM2 + 2*NSYMG 

I1I=C 

NO 4399 M=4,M5 

X72 .s “GM(M) + i 

NO 4590 T=1,_'GT 

IIT=IT+1 

SING(TTY=SINGY 4 FGAMCT)) 

SOSGUTIYSCISCY1 GACT) 

SINZSCITISSTNGCXL $1.9) *GAMUT)) 
COS2TCITI=ACOS COXL HL. 90 FGAM(TY) 

SONTINUZ 

IT=e 

DO 4500 N=1_4F 

Xi = NMPNOND # 1 

99 G5N9 Jai NPT 


Pisit+e i 
SINPETTPESTNCYTERETOCS)) ; 
COSPETIDECISEXLFEBETICIN) 
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5ifo 
Sann 


93Nn90 
9409 
9509 
9599 
a7rg 
aann 
2a9n 
aay” 
aurn 
2224 
2949 
aa61 
a34) 
9494 
qyi99 
Lers“ee 
yrann 

1 

? 


a 


SINZICTTVSESTNGA EMT HLA Fe BET IGS) ) 
COS2ICTTIECNTOCVL HL. MPFSET 205) ) 
SONTIN!!= 


GALL POST 


K=0 
OA=1.078 

NO 8799 T=2,4er 
D0 S799 Joi,irT 


IF (NUST¢SeTI2 = 2.20) G 


Kok +1 

DWwO¢<=0,9 
Juntey enn 
DHA CP HTLAR 

NO S199 MeLeMsS 
MM= (4-49 *NGT # T 
NO SIND NSte49 


TE CNUSTON MDC ET.99 GS 


NN=(N=1) ¥NRAF + J 
FGMN=F GCN gMD¥IA 
DWO(<) = FGMN¥FI1 (MM) 
JWG(<) = FOMNFE OO (4M) 
YWS3¢K) = FGMN¥F91 (MM) 
CONTINU 

SONTINUE 

NGNRST = K 

LMAX = LRARFRNGNAT 
LMAYST = LAAISTFNSTP#* 
NGN® = NG*N® 

RETUON 


FORMATULHL. 25% ,13HD F 
FO2QMAT (7H FIAT) 
FORMAT (AH SUeVEN) 
FOOMAT (22H METAL, S 
FOPSMAT (24H MOLASTIC, 
FOOMAT (L9H MOTAL, H 
FOOMAT (214 PLASTIS 
FORMAT (27H PLASTIC 


QO TO 5209 


9 TO 5190 


FFOOCINNY #€ OWIC4) 
*FPSINN)D # DWGE<) 
¥FPRCNM) # OWACS) 


N&3T 


PoRagd. FLUSHES ANE “SN eCY ZEON 


| 
| 
j 


INGLF LAYE) 

SINGLE LAYERS 
ONE Y* M2) i | 
o NONEYSOMAD | 
» MULTILAYE?) 


FOOMAT (27H = CLAMDER = CLAMPED, GAMMA OTRESTTOND 


FOPMAT (35H STMPLE 


FOPM\T (%6H CAPE 
EQRMAT (35H TLAMOFT 


FOPMAT (74H STMPLE 
FOQMAT (35H TL aMped 
FOPMAT > SHTOT EO ANSE A 


- SIMILS, GAY4A OTRECTIIN) 
~ STMPLE, GAMMA JIORSTIOND 

- CLAMPED, RETA IFIOETTTON) 

= SIMPLE, BETA JIRECTION) 

- SIMPLE, PETA NtPrsrroNn) 

PTTOMN - ELASTIC? 


FORMAT (RUHQOZSHONEE OPTTOM = ELASTIC-PLASTIC) 


FOPMATELZYNSTINSTUIAL MIE Ly 
U7H WUMRED AE GAMMA NOeoS EMGY = 1% 
h74 NUMPERNE BETA MONE S EMM} = Its 


> SESE AVAILABLE COPY 


ZY74X MUMAE? If GOMMA TNT=EG2Q TION POINTS (4B8ARD = T37 
& G7H NUMBFQI OF RETA INTE SPATTIOM PAINTS (N9APY = Tt 
5 4&7 WUMBEP OF ZF INTOGPATIION POINTS (LAAPY = TBP 

17210 FOOMAT (1H SF 
1 47K NIMBSFR NF STOINGERS (MASTRY = 137 
2 &7H NUMBER OF STPINGEP 2 INTESO. PTS.CL3L2STY = IB 


10 R29 FORMAT (2UHOMOTAL COMPINATIONS SEN) 
LNA FOOMAT 13K,2Th) 
17857 FORMAT (35HD LENGTH OF PANEL, IN (YLP) = £1564) 


19909 FORMAT (35H WIDTH OF PANEL» IN (THETA) = £16.98) 
11609 FORMAT (35H SUAITENDEN ANGLE, OFC (THETAO) = £16287 
2.35 PANTUS, IN CAI = €14.8) 

11479 FORMAT (75H VHICKNESS» IN = F15e8/ 
& 35H MENSTTY, L3=SECF*A7ING Fu = £15.97 
2 35H SEASTIC “MONULUSs PST = EL16.8/ 

3 35H POTSSON'S FATTO = €156.8/ 

& 35H YIEED STRESS s FST = F16.8/ 

o 1355 STPAIN HAPPENING SLOPE, OST = ELA.87 
5 35H ULTIMATE STIAIN, IN/SIN CEPSIF)= FLAP) 


{1277 FORMAT(PAHOINITIAL IMPFOFFCOTIONS, IN/(5F16.5)) 
11300 FORMAT (L7HOTTIME TNFORMATICN/ 


1 &2H TNTEGOATTON STEP SEZEs SEC (DELYIMI~ = "15537 
2 &2H SYOP TIME, SEC @TSEQP) = F15.9/ 
3 42H JOINT FREQUENCY (PRINTI = F15.8) 


12059 FORMAT (4NHOCONFXINATE SURFACE POSITION CH3ORI, IN E16.8/ 
L L£340L AVES NIVIED, 22% go KT157031X,6T15)9 

t2ihQ FORMAT (27H CUMULATIVE THISK NESS», INe 1 3¥ 5551525) 

12200 FOPMAT (T?H MASS DENSITY, LAI“SECE*2/SING#458X55E15, 5) 

17309 FORMIT (33H MYOULUS OF FLASTICITVY - ¥, PST»7X,6E&13.6) 

1?409 FORMAT (49H MODULUS: OF ELASTICITY << THETA, St 5515.6) 

12599 FORMAT (22H DYTSSON®S RATIO = XK 918% 6615.5) 

124679 FGOSMAT (26H PITSSON®S RATION = THETO,14X%.5215.5) 

172?7NC FOPMAT (31H TEMSTLE ULTIMATE STPESS, PSTs9X 6515.5) 

q{2?aNM FORMAT (5H COMPRESSIVE JILTIMATE STRESS, 9ST 45% 95-1526) 

17999 FORMAT (28H TENSTLE YIELD STRESS, 2ST et2%—. 5E1526) 

17990 FORMAT (32H PIAIMPRESSIVE YICLA STRESS, PSTe8X%_5£15.51 

12659 FORMAT (21H SHEADP MODULUS, OST, 19X,6515-5) 

23107 FORMAT (527HNTDRAN MONULUIS OF ELASTICTYY PAPALLEL YO 399E DEPTH (FC) 


te POT = F16.4/ 
2 624 S4HEN9 YOOULUS OF SOFT (500, PST = 
3 E14. 
& 624 CPF CELL STZE (NC), IN = 
5 £15.45) 

182P) FORMAT CGTHO®*® WAOMING ** INCONSTSTENCY IN SYMMFTOV) 
ENN 


REST AVAILABLE CCPY 


198 | 


—<crenecalll esnsasepeicinte ste ter mana 


nen ae enim 


SURPIUTINE OTSTHO (FZ, OP,e eK, eS) 


THIS o TF OOUTII= PaMoyree (AS SO9WYIVATS, CINSEOVATIYFE 
TIM® [YTD EN NHSEHIA99, 
NOTE: FY J* ASSUMED FHA = 92. 3 AY PCANEL THFTAD wT. KXke 


SOMMINSTOEL KIS Ng TYASTIGY7,_ BOA eb BAPTT AL MAK gL MAXST,. OCSTP (6) oMBy 
M 14>, “48ST, 459 MOMCL 3) ae 42 MIAME USE CLL FING NIARZNIN( 13) * 
NINO NAT ANT 2QUING NEN, NSU AMNGT, NOL TSNST Rg NSYMY, NSYMG, PT 


“SOMMINSTILKPZ JOTI OPTI CTLUL 3) CC2E1 3) COS EBD COHL13Dy 


WAwWIN er 


Ve 


CKEFY gOASVW7 BAD GEMEG (299) gF0S27%(729) ~C05265(299),3°2TNOOTI, 
FIT (CPF AAV HMAC AAD) gFEZ(P 99) eg FOG (LS, 29,FP5(299) ,F 95 (799), 
EMP (PFQQYP_LREOALY yr, 
GAM OPTV, KS, PIMA CPT) FINA C7T)D STAC 799) SING (799), 
SIN IEPA) SSTNAGEIIAAD gh eK SQ aK SF eX SGX Sg ML XL PgX_O1yXLPDA, 
MOPS s Kate Xb ee XL Se howe XL 7g ML STS Se SNS 
SOMPIN SZ RALKLSS ISTO FOSTF 4 SC eg HARA ENL » INOUT, QH%, FHETAN 
COMMON SINOVAS FOTTCS) DEL TIM, GAMMAC4H1) sg TOU, INIUT, KALT KB, 
KL EAS SKE PIGKNK AK TVET SNTALL gg NCASE gNOHOT INIBUGNYASS,NTIRIAL, 
PR (4°) , PRAM, 2PO,ORINT, PF 2,57 OTAL (5) ,TIME,TITLE L299, TSTOS, 
ETE C2) 


SN NAD 

IF CONS YS SCE Oe) EN = 2YNGBE 5 
THET? = OT¥T4FTANSLAN, 

QM = MOAD 

IF (NS¥Y¥G.ENM) OM = QeURD 
emsf,° 


TF (MAND CEQ SL 0 T%. NENDLFS.3 2 DR. NANDLEN6) CM=NLS) 
IF ONAN CF 925 NINOSET 7 0 O%e NINDLEDV29) CM=02135 
SN=0.0 

TF (NSN. E Get k NSNC2ECLY 09% MANOLEQSD Ch=NS) 
TF ONAN. ENeS NOMNN.GF.8) CNE9,15 


ariy 
OT2¥ 
YT 3y¥ 
OT4Y 
QTSY 


IF (% 


Nee 

1.MF4 

1.NE6 

1 .ntFh 

Tah 

Lig. & Wetthe - ST AX THETANFESQQITIEZISISN = 2M) 
TF % oF Ge be) SES AETAETIFSOSTIERISCIN = 1.9) 
IF OMPUT CE Ve FES UTX YUSS SIS CETL CEM = te} 
SPECS SLi MOS CORR INGE Pes 

PO 28S She SS 

SaAePv = “MEM(M) 

Or ai A SM 

IF C4ISE ON MICE I6%) 60 TO 

PAI = NONCN) 

AAPYY=7A IM & OM 

BANC = RASH © fPky 

PE StVeE eee be se To See 

FLAY feb 


o 
or 


Y 


Wei eben se 


2°0 


“4909 


ELMM = 
nT1 = 
31 = 


(PAPMKEOTsYL IPED + 


IF 
60 TT) 


Tax = OT1 


ran 


CUTV=S PANEL 


(JARNK*FOTS THE TAQ) ¥¥2 
Dr PSIOTCELMNECT2 & ELMNFETY)) 

IW 2568 
CIT1.LT.0TIYY 


# FLYN**#2) = O.5*SQPT((1.9 


ZLMN®*47 (FE Z* 3M) )D 


ELMMN = JARMY¥PTEASYLD 

EKMN = JARNYFOT/STHETS 

DUM = CELMN#F2 © EKMNEED) HED 

JTL = STFA FSOPTIE 309.5 %11.1 
2 ELMN**2) FO & SPF 0. TOL MND ¥¥79)) 
DTL = OF1775,9 

DPT2 = CT eA /SOQPTEEL EC UMS AFF? + 
OT? = OT2735..0 

FLMNM = SPT OFF 

97S = Pr SSORTCELMNFCE2 & EL MNE=5)) 
OTS = H1T5/25.0 

IF COTUL TOTES OTAX = OT1 

IF COT2.LT.ANTPN) OT2X = OF? 

ITF CYTSSLTLPTSY) ATSX = NTS 


CONT us 

TF CNPLTLENLN) 
TE €NPETSEQS1) 
IF (N3955.G6T.9 
PETUON 


FORMAT C3LHONESROSD TIME 


ENN 


Seer 
JELTEM 
»y WITTE 


al a Soe 
ys eee 
200 


QMINTONTLY,IT2X) 
AMINLCOTLY IT AK TTI ,OTEY IP 5K) 
5e19009 DWTAK,VT2X,OT3IX, DT aX, TSX 


STFP CALCULATTIOINSZ5SFE1L5.45) 


rh0 


SYROVITT ME MT CK MEY MOLIIM, TIME 9 VX 3_NG yg ANS) 


SPECT UL TYPOS ATT Oy METHOD £92 21 MPDER DTFFEPENTIAL FOUATIONS 
WHT T= SAVE NO JAMOTNGS TENT OAL OTEEPPFNTE SCHEME, 


COMM THIS KC RELEID OTS” 
yyMeEyYoe TM 2s oe Be Ve(L),vYr7(t) 


FE Cercle st oF TS -200 
JTSC = YELTYMer? 

¢ = 

OH one Hetensze 


YU(T) = KB(TY = VY¥DETPFOFL TIM & RSET TSQFOKS(T) 


MQ an) ¢ 
HoPome PCR GLA KICTE, © DTS ISSee OTs 
VSCTY = VRCTY 


4G CEY oa o¥ 

FINE Ss ELVE- 8 SU TIM 
SeTyo 

END 


201 


| 


SUSPYVITIME LIGCND 


THTS PUTING CONTAINS THE LESENNPE 7FO9e5 AND 
WETGATING FATTO FI® LAD AND LAAT ole. 1h. 


Psa 690) 


SIMM IN STOLKAS Ag TMASTI KZ LIA Mol LOST SLYMAK GLMAXST,_OCSTP CE) NB, 
1 MBE Ae MAS T2*"S,4GOMEL g) ah 4X MGM APEMUSE (13,513) oN BeVA4PO,NBN(13) o 
? MIN TZNRTE NOE OU SMG yeNGN®, NUNIT ING Te NPL TS NSTA ANS YMBGNSYMG,PL 

SOMMINSTILK3S 3 Y¥ (6) eGXYST2 (14) pHS9(6) gYGISTI014) 

NIMFNSTION CO¥OC7,14), CHO(7,14) 


| DATE TIYA7THGL I, V6 5773SI26513995 3,5 FNL Vy Ve 77HS9G6HIZ41 45 y 
Be Me BALL 3S 31159405, NM, F3FII1LTG 35°43 G, 5400744 906995179984593 965, 
DeFRUAAZLILO 55% SFM My Me FSPGAIFLNPT ALG» NOEHLATISRII“G526, 
Oe PPTL ALIA IDA B77 GF V0 1g Me WGALIW7AL2F4 2754 0.7415 31133559939,— 
MGT SHSASLI7 7H DV SET Ve MGA DAIIKSHGI7SG NA TIBHSAG77 GI IGS, 
QoS 7° FR2GNAN GVW, TMLLR WG 3454249555, 3* 909,49. 96 815923753 7ER, 
PeAZACZLAN73IPA AH, Te Ht 2371% 7270915947. FAU2SZG4AZLN TAL, TF0.0, 
0697 TANFAGAKGL 71749 %. SHSNL IZED 3399249 .56794N9558 79902, 

Deb SPIISIH1 2975, De LHAB74 3389315 2%, 7 F900 02978 22855314505, 
VNeA7™7O5HAGGG7TAIHIIA*NL 7ICLS2INSS 7HCS NW. SLINGISL2AI2]N5A1L, 

06 253743155952 Hg 2FT0 1, 0.9915 59% MU24H72,0.ING117255 37047, 
0.75799 PAULIN 39, V05R 772179547 855 29%, 357K 3149899913, 

MeL 25 2BBGI1BS1147 4 Ve Vy Ve AS84139 3054 710595 0.917592 39392229R, 
De TLS7ROGN7 TT 314 VAG2ILGIZIGEIIG YN. GUAUGITSINIGES, 

Ne TITUS Z®BLGIAGS 1G Me Me Me FBS2BTATAAFE 1 oD. A2RUIGIB ISIS 7, 
De R2A72AVLFLSNA IZA V.S87AIUP ING LL HAD ELG2U3IB3ARIAL5,y 

De FLPILPZHRIS7 AIAN ANANSY IGA7N7TZG/ 


NDMNFWNPRP DOAN AWE WYP 


DATA THISTEN LI Lo MA ED ol gf 0 555555 5555555F 97. BBRY388999ABB9, 
SEV Te, FETIEGIKRGLT 7S, 9, 55274 454 S485 255 ,5FNN De C3SIZGRBSNSH19, 
Me 7VAPIB7TAGIAAZ7Z VSHRARAAIHFIRRAWGFETLDANLI7L 324497237917, 
Ne TAPTALE7T3IVUG 21, VLE PAL SITES 72H IGEN Ny VAI 2AUIGIASL ABB, 
De S7IPTF VIG IPRA, FALIINIGNS Wt Py MUL7TISMAIA7 37, SEND, 
MoTVEPIRGTEZZI7AL 7, QAATALO TLE VT74N, ULS7OSAUSI77BI, 
Me TA PAB VIB IT7A Wy BEV Ve MAL 27GB ALS THE andy 1061396431695 9485, 
Ve PITA NAQRGT 793, De WIZ NZ77I G09, Ne FBO ATIZISONW1 2H, PFN, 
Ne Se ME SUK BDMA FE M4 94 IGT FHV SIG AG Ve PLAN IAA THP515 IN, 
Oe PR PFA5ZEFLGT TNT, 1, PAGS PUPILALLTELPENL DAN SSHAISELIS174EM1, 
Det P~ TAP PHOLE SAL, Me LP HAI 7FLII2773,%. 233197%76459199,~— 
Me PAPAMUGGGSL 256% 27AVP JARF77IIVI 4 Me Dee UPL7TSTVIABIS55L PEM, 
Met FS AFWBPGIAG AW NL ALANA ZS PIASKHTTG, VY, PBL ATH2577397,q 
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S54 (3619—9S7A (270) ,SK8A(230) ,UUC1L3 9130 oVV(13513) eWNI13,13), 
X3O7T) gp XGO7RZ) gXKT Ty XKXT p XKKX gp XKXKST XLA(551)5X200 361) XSAN 0361) 4 
K4OA CRG) XSAN CZELD KAA C3510 X7A C259) gXSA( 230945 78029, Z2ASTRI2), 
7302) ,78STRU2) »ZF (8), ZFSTIC1L4) 47606) »pZGSTR(14) »Z4(6V_,ZHSTRILG) 

NBC = N3ZUSECS,T) 

IF (NOEFV.EN.2) GO TI 45ND 

NDEPY = ie 

IF (NEC26T.I9N) GO TD 200 

IF «N2C.G¥.0) GO TN 10N 

NAC = =NRC 

VRX INAS) = CLEMGGGINSIC) + S2*WGBR INR?) 

ENX(NSSYSCMALEXLACKY + CML2FY2Z ATK) 

IF (LTeFQ.19 VAEXENACY = -VRYINSC) 

GO T) 12790 


we 


NOME WNP 


oF ] 


Qe 


190 VRTIN35) = C3¥WHIRPINBS) + 34*WGGA (NBD) 
ENTEN3SCY=CM22*¥ KPA (KD & CMLAFXL ACK) 
IF (JeFV21) VRTONAC) = =VPTINRC) 
30 T7 1209 


ea 


2090 NRC = NBR - 109 
RRIMN3IT) = COFVOCR(K) Flat. 0) FFINACS 2) 
60 ¥) L200 


NOFPV = 26 


cyaate 


ELASTIC-PLASTIC. 


70) 


6nn SuM = 9 


ae BEST AvwiuAdut COPY 
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SUMST?> 
SUMiST 
JI 


IF (NRC. 
N8C 
JIi 
K2 
IF (1.GT 
JI2 L3 
DELYX = 
EF €1.F2 
90 7oe « 
L JT + 
a | JEL 
Ee = -JT2 
G1! sxx 
Sii 
53 
G31 
SUM 
1 
700 SUM1 
IF 
JIST? 
K2ST2 
IF (T.GT 
SI2STP 
00 7530 
LST? 
L2sTf? 
G1iSTR 
SLiISf? 
SUMST 
1 
759 SUM1IST 
7990 VRX (NPC) 


= 13 


S¥T 
S¥ 


IF (TeEV. 1) 


EN¥Y«MECY 
GO 13 


Kl =K- 


SUM + 


TTMAST®.EQ. 1) 


JISTP + 


n.? 
rt. 


LAN2*(K=1) 
IF (N95.GT.ING) 


5A TO 100C 


GT.9) 69 TO 849 


-N8C 


RR¥ (<=2) 


PENT od 


ah <2 
AR¥(K? 


Cl-ShENAT: eo J 
1) 


QELX 


TEL XY 
1,LRA° 


ei) “DEX 


K 


KK 
+ K< 
+ K< 
tu) 


sxx (L2) 


tL) 

T(L1) 

(S511 GLIY/DELKXX + 2.9*(G3L 
GX (K<KY FHGO(KK) 


$3) /IELT)* 


SUM1 # HEV(KK)*G1 


69 TO 790 
LRAPST * (KSTR=1) 
T¥NSTP + KSTR 
ot} <2SUR (I-20 *FSFQ2 ¢ KSTO 
LBARST * (K2STP=1) 
€KSTI=1 gL AAQST 
KKST® 
JT2ST> #& KKSTP 
S¥XSTPELSTRI 
SXXSTPC(L2STR) 
SUMSTO + ((GLISTR-SISTR) SNEL YX) 
¥ HGISTRIKKSTP) 
SUMIST # HSOSTP(KKSTRY*FGISTP 
SS*SUM/26M% + Ce HSTR¥*27200) FSUMSTRI209 
VRXINRE) “VRY(N3C) 
H*¥SUM1/209 * HSTRFSUMIST/2.90 


¥ SXS TRICKS ERP 


1290 


2 


TF CJeFVe1) <3 K 

JIt LEAP KK 

DFLTT DOFLT 

IF (He EFe LF) SEULVY -DFLT 
JI? LIAPF((T-2VFNST & J-4) 
90 On" KK Lele 

JE RK 

STE @ BE 

JT? + KK 


ETT aaa eae aaa ine 


| 
7 
4 
4 
4 
4 
| 
7 


z 
9nd 


4e 


1nNQ 


1199 


L200 


| 


G2 = STT(L) 

$22 = SFT(L1) 

G3 = SXTCL) 

G32 = SYT{L2) 

SUM = SJM + (622 = G2Y7DELTT + 2.0*(G632 - 33) 70ELX)* 
GK CK) FHGOC KK) 

SUM1 = SUM4 # H59(KK)*G2 

VRTINPC) = S5*SUM/2.N 

IF (JeED01) VRTOINAC) = -VRTCNSC) 

ENT(N? SY = H¥SU41/2.0 

69 7 1709 


NBC = NAS = 100 

DO 4199 KK=4,L3A2 

L = JT & KK 

SUM = SIM # GXTKK) FHGO (KK) *S¥T CO} 
RRON3BCY = CS*SUMFE (1,0) FFONBOZ 2) 
RETUON 

END 


SUBRIUTTINE RELAXP (NEQ PES eX yERP NOK, NPRINTeNCOUNT) 
DIMENSION RESC1) KCL, FRPCL) 

TOMMON SEALEL2S PELX (447), TPLLG7) gPPTST147) PK I1L47)_ 2P2PESC147), 
1 SISYCLG7 Ig YPESCLG7 197) XXL C147) 

QATL CONS5SONN.S 

IF (NPOINT.NE. 2) SO TO 40 

WRITS (55509 (KON) SPFSTN)D SN=LZNEO?D 
FORMAT (1H 94/7%%53H TOTAL XglLOXe7HAFST WES USY FL 3.69 4K yF1 305)? 
IF (NOCINT.FO.M ¥ GO TO 10 

IF (NCTOINTeLE NED) GO TO 14% 

NCOUNT=N 

NC=0 

NNG=19 

GO 2 PSReh=9 

Ox = X(T) = OX(T) 

IF (ARSdAvy LE. 2P(TI) NQ = ND + 1 

TF (ARSCPES (TI) -ABSCPFIFS(TIY.GTe1.0? NNQ=ENNO EL 
CONTINUE 

IF (NDEQNFEI) GO TO 10C 

IF (NNO’FQeN=E0) GO TO 101 

NO & TH1l,gNF% 

XXL (LD) =¥ CT) 

RRESCT)=9FS(T) 

DELXCIVAABS (9.9901 *XCT)) 

TF (VELXCIVLTCERS(I)) DELXCIY=HERRIT) 
CONTINUE 

GO T3 6 

DO 7 MM=1,N=9 

XRESUMM,NCOUNT) = (RFSOMM) = RRESTMM)) Z NELXCNSOUNDT) 
XENCIUNT) =XXLUNCOUNT) 

IF (NCOJNT.EQ2N£9) GO TC & 
NCOUNT=NCOUNT#14 
X(NCOINTI=HXOINTOUNTIF9OELX CMCOINTY 

NOK= 

RETU2N 

90 29 T=L,N-9 

SIGX(T=-QQFS(T) 

CALL SOLVE C¥20S,NEQ,147 59 5F%, DET SIS) 
IF (NFMEN.09 GI TO 15 

P2OP=1.1 

90 17 Tei,N&9 

IF (ARS(SIGY(TIVLLTSCONFFP RARITY) GO TT 13 
X¥PPOF = CONFEPPCTISARSISTGYIT)) 

IF (xX PROOLT.P%09) PROP=xXP IOP 

CONTINUE 

N0 17? THleg"=O 

X¥CIY=¥XLCLOFS TOY CT) FPOoe 

Py(T)=x¥4(T) 

POFSTTI=RPES(T) 

NOK=1 

NOOUNTSNEM+4 


Mee se 


199 


191 
5S 
15 
i1 


RETURN 

NOK=1 

60 T3 it 
WRITS (6,55) 
FOPMAT (32HC 
NOK = 2 
NCOUNT=9 
RETURN 

ENO 


SOLUTION PIVERSING IN RELAX?) 
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a eee beet enens 


FUNCTION SEG (7U™) 
FINO ELAPSEN FO TIME. 


SALL SESOND (SEC) 
RETU2N 
END 


CPI CPRPICICILACIOCICI GC? 


4262 


te 


ee ae 


SUBSRPIUTIME STEMA (TeJeM, MST?) 


THIS S/2 NETSOMIMES THE STRESS STRAIN BPELATIONSHIPS FNP 
ELASTTC ANI/92 PLASTIC RESPONSE. 
SURBRIUTTNE COMPLETELY REVISED MARCH, 1976. 


K = TNB=X OF THE INTEGPATION POINT IN THE Z OTRECTTIN. 

I - INDEX OF FHE INTESPATIIN POINT IN THE SFTA JIQSSTION. 

J - INDEX OF THE INTEGRATIIN POINT IN THE GAMMA JIRESTION. 

KSTR - INDEX OF THE INTEGRATION POINT IN THE STAINGE2 Z NIRECTION, 


SOMMONZTALKIS = A,ZIMASTR KZ LBARLBARSTyLMAX sg LMAXST,.-ICSTRIS5) MB, 
4 MAAR, MAST24 MGe MGM(13) pMSM3_9 MGM 32—4USE (13913) ,NSeN3ARQNBNIL3D, 
2 MIND gNST ANS = AV gNGeNGNR yg NGNS Ty NGT eo NPLTy NST 2—gNSYM3yNSYMG,PL 

SOMMON/SOALKS/S = SXLAVSGKXSTE (LG) gHGO(6) gHGISTIC1O) 

COMMINJZTRLKGS NY 2eVXO(147) XX (14 70 9YVC147), 

1 ASIS(49,69), NAW (49,49) ,B3U(4 9), B5W(49) , TPULG9) , TOW(G9) 

COMMIN S/CSLKBS  ALTT€1805),4LXT(1895) ,ALXX(1895)N, 

1 ALY¥X¥ST (2390) ,361(01805) 86 201805),8F 301805) ,BEULLIBN), 

2 FIVLRNSIV,EPBST2(13899,5TT1 (19805) ,EXTL¢1895),EXX1(01805), 

3 EXXSTL(L 380) ,SIGTTLI18059e.SIGXTLL18295) gSIGXXL (1805), 

4 STGXLS (4.339), TTNUCLANS) gs TINUST (1389) 

SOMMIN/SRLK7S BST gEN1 eg CNL2eSN1L3 gCN2 g9ON2STR 9 CN3 CNY o CNS y CNG GCN7, 
1 Fug FT PeEPN, = PP,EOPSTA GH ZHSTR AIF IRST yJEIRST, ISTOFT SLE gLOMAX 
2 LIMAKSsLESTRNELO,SIGO,SIS9N 29 TNUy TNUSQ 

COMMINSSNOVAS = CRITCIS),DELTIMN» GAMMA(K1) g ICOMP, INOUT, KALT 9 KBy 
1 KIAMA KNS XE FP 2eKAK SK TYPF gNOALL g NCASE gNCH OT gNORUGSNYASS NTRIAL 
2 PR(4I»PORY%,FPP,PRINT,RF2,RTRTAL(S) »,TIMF es TITLE L299, TSTOP, 

77109) 

COMMON CN1M4CN11 eC NA CNIy FPSO (1805) ,EPROSTIL 3809 ,ETT EXT, 
FXY gi Y¥XST2, IN7(2) gINZST2(2) gKSUYA (351) 5<SUMAS (233) -KY(1805), 
KYSTICLZAOD gy NUSF (23923) gSTT (1805) gp SXT(19305) oe SXXUL805), 
S¥XST2€1 350) ,S1AC361)9 pS2NC351) S34 (3619,S4A(3561),554(361), 

SSA CPHL) gS7NC2ID) SHA L250) 4g UUL13 413) gp VVI1%413) pWACL 35139, 
YROPT) gX¥G(2T yy KKT Ts XKXT op XKK KM yp KXKKXST 9X1 4(0351)—9 X2A0 361)» XSA0361), 
YEN (ZBL XS CBGLIXBA (C351 VX 7ZACZATNYXYIA(CZINI,ZZA(29,ZASTREI2Q), 
7307) p7ASTICF)D -ZFEK) gp ZESTII1G) »7G(59 pZGSTI014) 92416) eo ZHSTRU1L4) 


WwW 


NOU LCWNP 


DATA TOLI5.NF- 7% 
DATO PANEL /1NH PANEL 4s 
1 STOTNG S194 STRPINGEPS 
IF (TFIZST.GT.9) GO TO 30C 
INITTOLTZATION SOUTINE FO? PANFL 
EF ERS 2.4 


FPOS3IG.V/EL 
Epp = nn 


BEST AVAILABLE COPY 
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SIGQN2 = STIGN¥*2 

CNL = (925 = TMUIZEL 

CNS = teNZEr 

TNUS?) = TNY*F2 

CNG (1.9 - TNUSOV#*?2 
CN1? (1.9 = TNU & TNUSTISTENG 
oN1LS (14.9 - WeM*FTNU + TNISD) SONG 
oNG N.757 001.9 + TNUD *F2) 
CNA CLste - TNUFRFQ) 

oN7 FL¥ON.57¢1.-0 * TNUD 
CNS 1.0/CN7 

Le = 0 

LOMAX = 109 

90 109 LeteLVMAyY 

IF CHELPLEQ.4? GO TO 70 
ALXY(L) 

ALTT(L) 

ALXT(L)D 

RFL(L) 

BE2tL) 

BESCL) 

TTNU(L) = 

EP8O(L)I=9. 

KY(L) = 1 


IF (NST?22E929) GO TO 300 


INTTIALIZATION ROUTINE FOR STPINSERS 


- 
v 
-”- 
. 

~ 
- 


EPpPSTP 0. 
cN2ST? 0.9 

LCSTR 

LOMAXS = inf 

NO 208 LST®=1 ,LMAYST 

IF (NFLP,FQ.1) SC TO 170 
ALXXSTC(LST2?) = 9.9 
BEG(LST?) = 9,7 
TINUSTILSTR) TNU 
FORQOTT(LST®) 9.0 
KYSTAC(LSTP) = 1 


STA2T OF MATIN ROUTINE FOR PANTL CALCULATIONS 
= 4 
L3UAR*{M-1) 


LI9> FAP CANEL CALCULATIONS 


Hi = 25K) 
DETEIMINE APPRAAITATS PEGION. 


KEY = <Y(L) 

IF (K€&FY.GT.3) GNX TN 350 

GO T) (470,600,709), KEY 

IF (CKEV 419 77.F70KF£Y/2) GO TO 600 
50 T9 710 


REGTIN 1-2 ELASTIC CURVE. 


KSUM = «SUM # 1 

IF (KS{M.GT.1) 3% TO 4590 
abe CNSe(CXX # TNUFETT) 
02 ONS* CE TT + TNUFEXX) 
3 CN7 EXT 

4 CNS OXKYX & TNUFXXTT) 
05 CNS *¥CXKTT + TNUEXKXYX) 
36 CNP EXKYXT 

S144) m4 

S2at4) 02 

S3A(4) n3 

SGA(M) ny 

SSA(M) 0S 

S6AtM) 06 

$1 pf + H1*OG 

62 22 # H1®95 

G3 = 23 + H1iFDS 

SIGR) = Gi*¥tGi - G2) & G2**2 + 3.0*538%2 
IF (NFLI6EN-2) 39 TO 470 
EPBO(L)=SIGR 

Go T) 7n00 

IF (312329.G6F.S1607) 50 TN 599 
EPBOC(L) = STGRD 

60 T) 3f"N 


LINE APLY TNTEPPOLATE ON SIGMA BAR TO CORRECT FO® OV=°SHOOT. 


KY(L) = KEY # ft 
SOSI3 = SQ@T(STIGA)) 


SPPTCEPANCLY) 

(SIGN = R2I/0SOSIG = 2) 
S¥*(L) * SL *(Gi = SKK CLD) 
STr cL) * BL¥(G2 - STT(L)) 
syrfti) * 81°(G3 - SXTC(L)) 
S41 

62 

(ot 

TNUY*G2) 

TNUYGL) 


B2 
81 
of 
G2 
G3 
SIGY¥i (lL) 
SIGTTiC.) 


T3 = CVN5*G63 

FXX1(L) = T1 

FTTZ(L) = T2 

EXTLELY = TS 

EPBN = SOOTCINL2FCTLF¥2 + T2**2) = SCNLSFTLFT2 & CNLET 3 #F2) 
FPR (L) = FPR2D 

EPRBO(L) s Fea 

IF (JFIRST.FIe 19 JFIRST = 1 

GO T3 3900 


REGIONS 2 ANY Ge. PLASTIC LOADING. 


FP89? = FPS(L) 
H2 Fy¥ # HIE¥¥<y¥Y = SELIL) 
43 CTT # HLEX<TT = BE2IL) 
H& CyT & HLFXKXT = BES(L)D 
oN2 TTNULLY 
II 
II = 
CN22 
EP30 SOPT (C0129 - CN2 & SN22)*(H2**#2 ¢€ HIFFZ) @ 
(12.0 = BeOFON? + CN22) FH2*FHZIZ(0 109 - CN22) "#2 + 
0675 FAK FF2/ (1.0 & EN2) *F2) 
DFLEP = EPRE - FESRINCIL) 
EPP = (FP¥NFLED + ELFFOROILIYPSEPAN 
IF ¢€TNUGTo0.0) CN2 = «5 = EPOFENL 
IF CABSCON2=TTNU(LI) -LT2%20095) GO TO 620 
IF (IT.5T.229) GT TN 615 
TTINU(L) = GN? 
GO TT) 619 
WRITS (5555099 IN2,TINU(L),TIME,PANEL 
GO T? &190 
CN2 = TINUCL) 
IF €£°RN,LE.EPR99) GON Ta £50 
FPS(L) = FRPRN 
CZEPOLGCT CEL ee D2 eFPPLL Te FO) SD TO 4090 
=PI/701.9 = CN2**2) 
N.5*FPP71129 + CN?) 
Si®(H? ¢ SN2PFRTD + ALYK(L) 
S4¥ (HF * CNPFH?Y + ALTTIL) 
S2¥*4UK & ALYXTOL) 


SECOND TEST F9> UNLOADING IN EITHER 2EGTON 2 OP eo 


Cdeana 


QI=YLACMPMKEXXL (LD & HIF KGL (M4) 
Q2=X2 AMIE TTLOLY + HIF XSACM) 
DF=K 7AM MRE VTL CLD & HP FYSA(M) 
EF (2%059.029950 TO 66N 

P4=CVVELIPNRSTSYVYIIL) + ALYYEL) 


Cr eatr 6? 


P2=STT(LV-STSTTLELY # ALTT(L) 
P3=SXTCLIY-SISNTLOLY # ALXT (LI 

GO TT) 479 

P1=0.7 

22=0.1 

>3=0." 

El=Q1 - H2 

E2=02 - H3 

E3=03 - H& 

GLI=PL-CNBFC ELE TNUFEQ) 
G2=P2-CVSFCF2ETNUFELD 
GI=PI-CN7*ES 

AL=S1-%1 

42=62-5? 

AZ=G%=-PR 
SIGPISSLECALHA2VEA2IFF24 3, CEASE F? 
IF (STGID.GEeSTGI2¢6ANNMELEPLGELO DN) GO TO 53N 


KY(LI=KTYVAl 
TTNUCL)=TNU 
EPBO(LI=SIGA 
SEL(LY=71 + JFLCL) 
BE2¢{_)=]22 # BEAL) 
BES(LI=33 & FECL) 
IF (2P.29.09.9G9 TO 300 
ALXX(L)=P1 
ALTT(LY=%2 
ALXTILI=P3 

GO T7 3090 


REGIIN 3-2 ELASTIC UNLOADING - ®F LOADING. 


RBofLkh)y = FX = HLEXKXX 
REP(LI = FTF = HL¥XKTT 
PEZ(LY = FXT = HL *XKXT 
ALX XCEL) 
QuTT(L) 
ALXT(L) 
Ci CNE*CEL + TNUFF2) 
er? CNOFLE2 + TNUFEL) 
uo GAT SSS 
G4 Ci 
52 C2 
33 C3 
AL* CAL |= AZ) + A2*ES? + J. NFAZE¥D 
IF (S3IG30-6T.81592169 TH BON 
FPR3C(LISESIGRS 
60 T2 3000 


LINFAS2LY INT=FPOLATE ON STSMA BAP TO CORREST FOR OVERSHOOT. 


~ 
. 


ann 


829 


835 


84 


aR 


B82 = SN2TCFPIOCL)) 
SQSTS = SORT(STSG80) 
IF (32.3T7.ST59) GO TO 940 
nC = 0 
81 = (SYSIG - SIGOIZ(SQSTG = 32) 
NC = NO + 1 
IF (NC.S57.25) GY TO 820 
QEL1 = 31i¥*(GL - S¥XEL)) 
DFL2 = 31*(G2 - STTCL)) 
= 31*(G7% - SXT(L)) 


am 
m 
pe 

‘ai 


on 
a 
he om MoM 
5) 
mo) 
‘ 


con 
~) 


Q 


DELS 

Gi - ALYY{(L) 

5? © ALTYT(L) 

*3 = ALYT CL) 

SQSTS = SOPTURLFE(OL-A2) + A2F¥2 + 3,0 ¥A3%*2) 

IF (O8S50SQOSIS-STSOV/SIGO.STATIL) GO TO 820 

GO T) 835 

WRITE (595700) “W5yKe Te Io KEV SOS TG, 81,32, TIME, PANEL 
Le = to + 2 

TIF (LO.STALEMAX) GO TO G1IM 


QR 
~ 


pb 
aN 
1 
eo] 
nA 


= 51 = S¥¥4U) 
O&U2 = s2 = STTCL) 
DELZ = 353 = SXT(L) 
TL = KLACM) + CN3*(DEL1 - TNUFNEL2) & HL*XGR UM) 
T2 = ¥PA(MY # TNZ*F(DSL2A - TNU*DEL1) & HL*X5SA(M) 
T3 = Y3SACMY  CNSFDELS & H1*X54AUM) 
30 T7 RAN 
WRITS (6952909 Ky Ty Io KEV eo TIMES 32,5 °STSy PANEL 


T2 = ©XY # HLF XKXY 

T2 = FETT # HL*FYKTT 

T3 = FXT # HL¥X<KYXT 

LG =.tce#i 

TF (.CeS5TeLCYMAX) GI TO &109 
EX¥S(L) = T! 

ETTI(LY = T2 

EXTI(LY = T3 

He =) Wi=" Bette} 


H& = T3 = SES(L) 

EPRD = SORTUINI2F(H2F#2 + 43¥¥2) | CNIZEH2E4S & CNYFHY FF2) 
EPSC(L) = FPRSN 

FPS(L) = FPPN 


STIGY¥1(L) = 3! 
SIGTTi¢L) = 6? 
STsyT4 (Ll) = 3 
KVCER = KEY ©. 4 


50 73 4908 


233 


Ce am a i a lath 


saci lt de 


LAST PAPT TF PANTL LOOP 


C302 62 


SX¥tL) 
STTtL) 
SXT(L) 
CONTENU 
KSUMICM) = KSuM 


e4 


32 
63 


‘9 


IF (IMASTP.LED219 GO TO 3993 


START OF ROITINE FO® STRINGER CALCULATIONS 


Ca od OD 2 


<SUYST = 90 
TJSTZ = LAAPSTF#(MSTO=“1} 


LOOP FOP STRINGER CALCULATIONS 


cqacaa 


00 3350 KST?=1,LBA2ST 
LSTR = TIST? + <STR 
H1iSTZ = ZFEST2CKSTP) 


KFYSTR = KYSTPEC(LSTR)Y 

IF (KFYSTR.~GT.3) GO TO 3100 

GO TV C3299, %5910,3700), KFYSTR 

IF €(KEYSTR#LI/2.EQKFEYSTR/2) GO TD 3500 
GO T3 3700 


REGION 23 ELASTIC 


sacaecaaa 


KSUMST = KSUMST + 1 

IF (<SUMST.GT.1)9 GO TN 3250 
DOF = ELPEXXST9 

D8 = FLEXKVYST 

S7Q(MST2) = 37 

SSACMST2) = VA 

G& = 07 # HiSTP*DOR 

SIGRIS = G4*5% 

‘IF (NELPLFQS?9 GO TO 27790 
EPROSTILSTRY = STS3NS 

GO Tr) 3900 

IF (SIG3O0S.65-eSIGO2) GO TC 330N 
EPROSTILSTP) = STGBNS 

69 TJ 3390 


LINFAPLY TNTFOPOLATFS ON SIGMA 849 TO COPRECT FOP IVERSHOOT 


KYSTSCLSTRY = KEYSER € ¢ 
SQSI3S = SQPTCSTGEADS) 


RT <<x«u 


4acacree 


3513 


84 = SQITLEPRASTILST2)) 

AZ = (STS - RG4IS(ISQASIGS = 349 
GH = SYXSTRILSTR) + RRFIGH ~ SXXKSTRILSTR)) 
SIGYIS(LST®) = GG 

T& = CN3*G4 

EXXSTI(LSTP) = Th 

FOROST = ARS(T4) 

EPSSTRILSTRY = FPROST 
EPROSTILSTR) = EPROST 

IF (JSTRFTLEQ. 1) JSTRFT = 1 
GO T? 3900 


ZESGIONS 2 AND Gt PLASTIC LOADING 
EP3MIS = FPASTO(LSTO} 


H5 = FXXSTQ + HISTREXKXXST = QEL(LST?) 
CN2ST2 = TINISTILSTR) 


TIsT2 = 9 

EESTI = LEST + 4 

CN2ST2 = GCN2STO¥¥2 

FPROST = ABSTHS) 

DELFPS = FPRIST - EPROSTILSTR) 


EPPSTR = (EPFRELEPS * EL YEPROSTCESTRIN/ESBDST 

IF © TNJ*GTI2 19 CN2STR = 0465 = EPPSTRFEINIY 

IF (ABSCCNZSTO=TINUSTILSTPD) LTO 299059 GO TO 3520 
IF (IISTReGT.299 GO TO 3515 

TINUSTILST2) = CN2ST2 

GO T9 3519 

WRITE (555509) CN2STR,TTNUSTILSTOV,TIME,STRING 

GO 73 4ngn 

CNZSTR = TINUSTC(LSTR) 

IF (=PBIST.LE-EPROPS) GO TD 3559 

FONSTIILST2) = FDAOST 

TF © EPISTReST cE «ORs FPPSEPSLTSFRP) CO TO) GaSN 

S3 = FP®STR 

34 = S3FH5 + ALYYST(LST2) 

59 T7 39" 


SECOND TESY FAR UNLOACING 


OL = X7ACMST2) = FYXSTLELST2) & HISTRFEYBACYSTRD 

IF (FP.EN.N.NY SO TD 3569 

Pu = SXYXSTRILST2) = SIGKIS(LST2) + ALXXSTILSTRY 

GO T9 3570 

0.9 

Qe = H5 

Ph = FLEES 

54 - PG 

STGO9S = A¥FOG 

TF (STGINE.GFLSTSO? JANN. DJELEPS.GEL1.D) GO TO 3530 
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3700 


3800 


R20 


3330 


3R3I5 


3RAGO 


zaag 


KYSTQ(LSTPY = KEYST? + 1 
TINUSTCLSTPY = THU 
FPROST(LST2) = SIGROS 
REG(LSTR) = Qe + REGILSTR) 
IF (2Pe=Q0.9-99 GO TO 3900 
ALXYXST(LSTQY = 24 

GO TJ 3399 


REGION 32 ELASTIC UNLOADING--RELOABDING 


EU = BES(LSTOD = EXXSTR = HLSTRA*XKKXXST 
Ch = ALXXST(LST2) 

S4& = C4 - ELFEG 

AL = G4 - C4 


STGRPIS = A‘*IG 

IF (STGANS.GT.SIGN2) GO TO 330N 
EPRCSTILSTR) = SISB8NS 

50 1) 3990 


LINEARLY FNTEZSOLATE ON SIGMA FAP TO COP2=CT FO? JVERSHOOT 


BY = SQRTCEPZOSTILSTPD) 

SQSI3S = SMPrT¢siIcans? 

IF (384.3T.SI60) GO TO 3849 

NCST2 = 9 

833 = ({SASTGS - SIGNIZ(SQSISS - 849 

NCST2 = NOST? + 1 

IF (NCSTRGT.10) GO TO 3830 

JELG = 33*0G4 - SXXKSTRILSTQ)9 

G4 = G& - BELG 

Q4& = G4 = ALXYXSTI(LST?2) 

SQASISS = A& 

IF C(ABSCOCSOSISS=<STGOIZSTGO GT. TOL? 39 TO 3820 
GO TT 3835 

WRITE (8e57019 NCOSTRAKSTO gly Jeo KFYSTR ye SISIGS BS 9345 TIME STRING 
LOST? = LEST? + t 

IF (LOSTOA.LGT.LCMAXS) 59 TO 4N30 

CONTINUE 

DELL = SH = SXYSTRILST®) 

TK = X7ACMST?) # DFLGZEL + HIS T?*#XBACMST2) 

30 TD 35AR 

WOTTE (555209) KSTPT, Je KEYST2,TI ME, I4,SOSTSS,ST2IING 
Th = EXXYSTP & HISTPFEXKXXST 

LOST2 = LOST? + 1 

IF (LOSTOSGTeLOMAXS) GO TO &NDN 

EXXSTLTILST9) = TH 

45 = TH = SFL(LST2) 

EPEPNST = ASSTHS) 

FPROSTI(LS TOY = FHNNST 

FPBSTIILST2) = FSANST 
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SIGXLS(_STP? = G4 
KYSTUUL3T2) = KEYSTP + 1 
60 T2 2390 


a LAST PACT JF STOITNGEP LOQ5? 
2 
WG) SKYSTII_ST29 = Th 
FASH CONTINGS 
<SUMASEFSTR) = KSUMST 
se SALLULATIANS COMPLETE |-= RETUPM TON JERY? 


7903 2ETUIN 


Ve i  ) 


EPROR PZ TURN, 


p02 02 


GON WOITE (54,5391) TPOLK Ty Jy TINE, EOIN, EPBNPEPINILD y PAVEL 
r GO 1) 4105 
&N5D HOITE (4953099 TPPSTREKSTR ely Jy TIMEy=PIIST,£ OS IPS, Ee PBOSTILSTRD, 
1 STRING 
&O90 WRIT= (45,5197) 


GL1iN0 WOTT= (45,5497) 
KERR = 1 
RETUAN 


.? 


S100 FOOMAT (7LHOSTOTNSER CALCULATION) 

SAVY FOPMAT (22H TMMENTATE PELOODIING 9&7 34 %E1525, 419) 

S390 FOOMAT (P7BHATPO TS OUT OF FANZE, FPP = £14.57 
i BT bT15 6 ,A11)? 

5407 FOPMIT (224HOSOLUTYTON ITS UNSTAYLE) 

5590 FORMAT (76H VALI= OF NU WONT TONVERG= 6 2F1505915H FIME, SES = 

. 1 £€15.5,410) 

S729 FOOMAT (3AH SAN NOT TOTALLY SO 2PECT FOQ OVE 2SHOIT/51T544F15.5 A109) 

ENO 
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SUBRIJUTINE SOLVE (A,N,NOIM,NIET, LT PDH T,S? 


A = SPIGINAL MATRIX. 
ASTUAL DIMENSIONS OF A, 
QECLARE NATMENSTON DF A IN CALLING PPOGRAM, 
NETEPMINENT CONE. 
NM = NOT GALCULATSC. 
1 = CALCIJLATEND. 
EP TNREY OF KeTH PIVOT 20W. 
DET JETEPYINENT OF A. 
B 2IGHT HAND SIDE VESTO?, 


ey Ca cr ca ca ca casa cerca 


DIMENSION A(NOIM,NDIMD IT PINIIM) 4® (NOIM) 
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TPCND=t 

DO 6 K=1,N 

IF €K.FILN) 350 TN S 
KP1=<#l 

M=K 

90 1 T=¢71,N 

IF (CAPSTCA(T 4X99 .GTLABSIA(M,KIDY Y=T 
CONTINUE 

IO(K)=™ 

IF (4.NE2K) ITOCN) =-IPCN) 
T=A(%, Ky 

AIM,<d =A (Ky, K) 

ACK,<)=T 

IF ¢F¥.F9.0.0) 69 19 5 
00 @ T=KO1,N 

A(T 9<d=- ACT <I ST 

30 & J=<P1_gN 

T=h8(4%, J} 

A(IM, J) =O (0K, J) 

A(KeJ)=T 

IF (T.FI29.09 GO TD G 
DO 3 T=<PtgN 

ACI gJV=HS(T, SDA T KD FT 
SONTINUE 

IF €4¢K,KV.EN0029) GO TO 15 
SONTINUE 

IF (NOQETSEQe9) GO TH TL 
NMET=xTO 

90 19 I=1 N 

OFTHIS TFACTsI) 

IF (N2E7219 GO TD 14 
NML=N-1 

90 12 “=i ,NM1 

KPL=<41 

M=IF(¥) 

T=2¢4) 

30M)=3 (4) 


B(K)=T 

00 12 T=KP1L.N 
B(TIHRCTIGAIT XI FT 
30 13 KR=1,h1 
KM1=N-K3 

K=KMiL #4 
BOKI=ATMLIVZA CK, HD 
==-P(K) 

90 13% Te1.K™1 
ZCTI=ATUTIEACT “IFT 
RCLI=FCLISAC1L, 19 

SO: ¥2 47 


N= } 
RETUAN 
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